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Abstract 



I modify the quasilocal energy formalism of Brown and York into a purely Hamil- 
tonian form. As part of the reformulation, I remove their restriction that the time 
evolution of the boundary of the spacetime be orthogonal to the leaves of the time 
foliation. Thus the new formulation allows an arbitrary evolution of the boundary 
which physically corresponds to allowing general motions of the set of observers 
making up that boundary. I calculate the rate of change of the quasilocal energy in 
such situations, show how it transforms with respect to boosts of the boundaries, 
and use the Lanczos-lsrael thin shell formalism to reformulate it from an opera- 
tional point of view. These steps are performed both for pure gravity and gravity 
with attendant matter fields. I then apply the formalism to characterize naked 
black holes and study their properties, investigate gravitational tidal heating, and 
combine it with the path integral formulation of quantum gravity to analyze the 
creation of pairs of charged and rotating black holes. 1 show that one must use com- 
plex instantons to study this process though the probabilities of creation remain 
real and consistent with the view that the entropy of a black hole is the logarithm 
of the number of its quantum states. 
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Chapter 1 
Introduction 



The conservation of energy is one of the most fundamental ideas in all of physics. 
As a principle, its history dates back three hundred years to Gottfried Wilhelm 
Leibniz whose philosophy of nature led him to propose that kinetic energy (which 
he called vis viva or "living force" ) is conserved in an isolated system. This notion 
gained currency as the eighteenth century progressed and gradually widened to 
include various types of potential energy so that by the end of the century, the 
notion of conservation of vis viva was essentially equivalent to the conservation of 
total mechanical energy. 

At the turn of the nineteenth century however, there was a large gap between 
the prevalent idea of energy conservation and how that concept is understood today. 
Most conspicuously, people didn't realize that heat was a form of kinetic energy 
and instead believed that it was an independently conserved immaterial fluid called 
"calor" . In fact Sadi Carnot developed his theory of heat engines based on that 
concept and it was believed that the power of steam engines originated from flows 
of calor from high to low temperature just as water wheels are powered by the flow 
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of water. It wasn't until the middle of the nineteenth century that James Prescott 
Joule performed his decisive experiment to demonstrate that heat was also a form 
of energy and so paved the way for the modern formulation of energy conservation 
as expressed in the 1850's by W. J. Macquorn Rankine's definitive statement that 
". . .the sum of the actual and potential energies in the universe is unchangeable 
. . . " . His notion of actual energy was identical to kinetic energy and in writing 
those words he was fully aware that heat is a form of kinetic energy. 

With the recognition that heat is energy, a theory was required to explain under 
what circumstances work could become heat and vice versa. Recasting Carnot's 
theory of heat engines in the hght of the new ideas, scientists such as Rudolf Clau- 
sius and William Thomson (Lord Kelvin) developed thermodynamics to meet this 
need. As part of this science and to explain why a given amount of heat cannot be 
fully transformed into work, they developed the notion of entropy and the second 
law of thermodynamics which says that the entropy of a isolated system can never 
decrease. Then, to try to explain the macroscopic and general laws of thermody- 
namics at a microscopic and mechanical level, physicists such as Clausius, James 
Clerk Maxwell, and Ludwig Boltzman created statistical mechanics in the second 
half of the nineteenth century. 

Now apart from emphasizing the central role that the notion of energy has 
played in physics, what makes the preceding bit of history relevant to this thesis is 
that the classical statistical mechanics that was developed could never fully explain 
the reality revealed by experience and experiment. For example, with its ideas 
of equipartition of energy this statistical mechanics could not successfully predict 
the low temperature heat capacities of an ideal gas, or much more dramatically, 
explain why the heat in a closed container doesn't all shift into ultra-high frequency 
radiation (the "ultraviolet catastrophe" ) . Problems such as these ultimately led to 
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the conception and birth of quantum mechanics^. 

Finally, momentarily leaving aside the thermodynamics, no discussion of energy 
is complete without a mention of Emmy Noether's celebrated theorem which states 
that the conserved quantities of a physical system are in one-to-one correspondence 
with the transformations which leave the value of its Hamiltonian (or Lagrangian) 
invariant. In particular, energy is conserved if and only if the Hamiltonian exhibits 
a time translation symmetry and angular momentum is conserved if and only if 
there is a rotation symmetry. Mathematically this result is straightforward and 
indeed almost trivial, but from a physical point of view its influence on theoretical 
physics since it was introduced early in the last century has been profound. Not 
surprisingly it will show up in my discussion of gravitational energy in the following 
chapters. 



1.1 Gravity, energy, and thermodynamics 

Today, the situation in gravitational physics is in some ways analogous to that 
of physics in general near the end of the nineteenth century. For almost thirty 
years physics has had a theory of black hole thermodynamics. It originated in the 
early 1970s with Bekenstein's recognition that if the temperature of a black hole is 
proportional to its surface gravity and its entropy is proportional to the surface area 
of its horizon, then the laws of black hole mechanics are laws of thermodjTiamics 
0, H. The first of these speculations was confirmed by Hawking's discovery that a 
black hole emits radiation as a perfect black body with temperature proportional 
to its surface gravity and the second supported by calculations which used the 
""^A more complete discussion of the development of all of these ideas can be found in any 
history of physics. See for example |Q or [Q. 
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Euclidean path integral formulation of gravity (proposed by Gibbons and Hawking 
in ) to predict that a black hole has an entropy equal to one quarter of its surface 
area. 

Thus, the classical laws of thermodynamics were extended to black holes with 
semiclassical calculations to bolster their interpretation and application. What was 
missing was a full theory of quantum gravity that could generate a statistical me- 
chanics to explain them at a microscopic level. Now, the difficulties in constructing 
such a theory need not concern us here but the important point is that any suc- 
cessful candidate must have those laws as one of its predictions. Indeed in the 
quest for a theory of quantum gravity, the laws of black hole mechanics are one of 
the few clues to its final form. Recently the two leading candidates, string theory 
and canonical quantum gravity, have passed muster and predicted the entropy / area 
relationship (see for example and respectively) but the issue is by no means 
fully resolved. 

As such, a proper formulation of the laws of black hole thermodynamics remains 
of great interest and that is one of the reasons why a good definition of energy is 
important in general relativity. To someone who is not directly involved in the field 
it would probably come as a surprise that such a definition doesn't already exist. 
Afterall, I have just finished emphasizing how central is the role of energy in physics 
and general relativity has been part of that science for over 80 years. Thus, there 
has been no shortage of time in which to investigate how energy fits into the theory. 
What is more, the energy contained in the other fields of physics is well-understood. 
In general, all aspects of the energy content of a non-gravitational field may be 
described by a four-dimensional stress-energy tensor Tq,^. Roughly speaking, at 
any point in spacetime, the time-time components of this tensor define the field's 
energy density, the time-space components describe the momentum carried by the 
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field, and the space-space components describe the stresses associated with the 
field. Indeed, such stress-energy tensors play a central role in determining the 
general relativistic curvature of spacetime according to Einstein's field equations 
which say that 

Ga(3 = SnTais, (1.1) 

where the Einstein tensor Ga/3 describes the curvature of spacetime. Of course, 
in general relativity gravity is curvature so the equations roughly say that matter 
curves spacetime and so creates gravitational fields^ 

Seeing these field equations one gets the first inklings that there might be a 
problem in defining energy in general relativity. Gravity plays the dual role of 
being a field and determining the spacetime in which it and all other fields live, so 
it seems likely that there could be problems in isolating its energy. Still and all, it 
seems possible that a stress-energy tensor could be conjured from somewhere. Such 
hopes are dashed by the equivalence principle. Recall that this states that there is 
no way for an observer making measurements entirely at a single point in spacetime 
to distinguish between her own acceleration and the effects of a gravitational field. 
Therefore there is no invariant way for a single observer to assign a "strength" to 
the gravitational field at a point and by extension no way to assign it an energy 
density. Thus there is apparently no way to define a purely local energy for gravity. 



An extended discussion of this point can be found in section 20.4 of reference [74 



How then are the laws of black hole mechanics defined if there is no way to define 
energy in a spacetime? Well, the answer is that the prohibition against a purely 



advisedly use the word roughly here since this split is not so clear as it might first appear. 
Namely, since Ga[3 defines the geometry of spacetime it defines the background in which the matter 
dwells. So, this set of equations is much more complicated than those of, say, electromagnetism 
where electric charge determines the electric field over an immutable background space. 
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local definition of energy does not extend to the total energy of a spacetime. At least 
for asymptotically fiat spacetimes, there are well known and accepted measures of 
the total energy such as the ADM |l|] or Trautman-Bondi-Sachs masses and it is 
usually one of these measures of energy that is used in the traditional formulations 
of black hole mechanics. However, this is not really a satisfactory way of proceeding 
since the thermodynamic system of interest is supposed to be the black hole itself 
rather than the entire, often infinite, spacetime of which it is a part. As an example 
consider a black hole spacetime which also contains a sprinkling of regular stars all 
situated many light years away from the hole and each other. Then, no one would 
argue that the stars should be considered as integral parts of the black hole system, 
yet the ADM energy would include the masses of those stars. Quasilocal definitions 
of gravitational energy attempt to meet this concern by defining the energy of just 
a part of the full spacetime while not attempting to fully localize the energy in a 
stress-energy tensor^. 



1.2 Quasilocal energy 

A quasilocal definition of energy is a procedure that associates an energy with each 
closed and spacelike two-surface in a spacetime. Though there are many definitions 
of quasilocal energy in the literature (see for example [^, ^ |2^, 0] and references 
contained those papers) a large subset of them can be characterized as Hamiltonian 
approaches. That is, they start with a Hamiltonian functional for finite three- 
surfaces in a spacetime which will generate the Einstein equations in the usual 

Hamiltonian way. Then the energy of the finite region is taken to be the value 
■^An alternate view on this point can be found in reference p5|| which argues that any quasilocal 
energy is equivalent to a stress-energy pseudo-tensor that fully localizes the energy. 
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of that Hamiltonian evaluated thereon. Usually all bulk terms of the functional 
are proportional to constraints and so its numerical value evaluated on-shellQ is a 
functional on the boundary two-surface only. Given this property, the energy can't 
really be said to be associated with the three-surface but instead is a property of the 
two-boundary alone. Any number of three-surfaces could be associated with that 
boundary, but which one actually is is completely irrelevant to the final evaluation. 
This property is in accord with the equivalence principle prohibition against a point- 
by-point localization of the energy. Since the energy can't even be associated with 
a particular three-volume, it certainly can't be assigned to individual points. 

In some ways these definitions can be thought of as analogous to the Gauss 
law for electric charge. Just as that rule defines the electric charge contained by 
a closed two-surface from measurements of the electric field made at the surface, 
the quasilocal energies define the energy "contained" by a two-surface based on 
measurements of the gravitational field made at the surface. 

One of the main aims of this thesis is to extend and generalize the popular 
Hamilton- Jacobi definition of quasilocal energy that was originally proposed by 
Brown and York Advantages of this definition include its appealing geometric 
form (discussed in some detail in chapter ^ and its natural interface with the 
path integral formulation of quantum gravity which allows one to do gravitational 
thermodynamics (briefly discussed in the next section and chapter ^, and in more 



detail in references pT| , |T^, 0). Further, in common with other definitions of 
quasilocal energy, it can be shown to behave in ways that one would expect an 
energy to behave. For example it is additive, negative for binding energies, and in 
the appropriate limits (and spacetimes) it reduces to such total measures of energy 



as the ADM energy p2[, the Trautman-Bondi-Sachs energy |TH], and the Abbot 



*That is for solutions to the field equations. 
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Deser energy |T^. In the small sphere limit in the presence of matter, it can be 
shown to recover intuitive notions of matter energy density . 



The first part of this thesis reformulates the Brown- York definition into a pure 
Hamiltonian form and removes the slight dependence on the spanning three-surface 
from which their energy suffered. That is, I modify the Hamiltonian they pro- 
posed, show that it does indeed generate the correct field equations in the usual 
Hamiltonian way, and further show that its numerical value depends only on the 
values of fields at the bounding two-surface in a way that doesn't care about what 
three-volume it contains. From there I show how the value of the Hamiltonian 
does depend on the motion of the observers measuring it, allow for the inclusion 
of Maxwell and dilaton fields, and show how the energy can be defined from an 
operational point of view. 

Moving away from the mathematical formalism I get my hands dirty and try 
to develop an intuitive feel for the quasilocal energy by examining the distribution 
of energy in the standard static and spherically symmetric spacetimes. I then 
investigate naked black holes and calculate the energy flows that occur during 
gravitational tidal heating. 



1.3 Path integrals, thermodynamics, and quan- 
tum tunnelling 

The last part of the thesis deals with a quantum application of the Hamiltonian 
work. As noted above, the quasilocal formalism of Brown and York naturally com- 
bines with the path integral formulation of quantum gravity and thereby gives some 
insights into gravitational thermodynamics. Recall that in general, path integral 
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versions of quantum mechanics calculate the probability that a quantum system 
passes from an initial state Xi to a final state X2 by considering all conceivable 
"paths" that the system can take between the two states (not just those that satisfy 
classical equations of motion). The action of each of those paths can be computed 
using a classical Lagrangian action functional and then, up to a normalization fac- 
tor, the probability amplitude that the system takes a specific path is e~*^ where 
/ is the action of the path. Then, the sum of all of these probability amplitudes is 
the probability amplitude that the system will pass to the final state X2. There are 
a myriad of unsolved problems involved in rigorously defining these integrals, but 
nevertheless history has shown that many physical insights can be gained through 
their judicious use. 

The problems of mathematical rigor are even more serious for path integral 
gravity than for regular quantum mechanics, but all the same its usefulness as a 
conceptual and provisional computational tool remains. In particular, as is usual 
with path integrals, one can use it to study thermodynamics by viewing the "paths" 
as members of a thermodynamic ensemble and so reinterpret the path integrals in 
terms of partition functions. 

The connection with the quasilocal formalism arises because the classical be- 
haviour of a system is not sufficient to specify the action functional that should be 
used to assign the probability amplitude to each path. However, it turns out that 
the choice of an action functional also corresponds to a choice of restrictions on the 
ensemble of paths considered. The Brown- York formalism provides a convenient 
way to see those restrictions from a thermodynamic perspective. With this insight 
one can associate each action functional with a specific thermodynamic partition 
function (for example grand canonical, canonical, or microcanonical) as was first 
discussed in reference [^. 
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That said, one can also use the path integrals in their original form to estimate 
the probability that a quantum event will occur. In this case, one must recognize 
that the action functional still places restrictions on the physical properties of the 
paths considered and so should be chosen to conserve essential physical properties 
(such as the angular momentum or electric charge of a spacetime) through the 
quantum transition. 

As a specific application, in recent years there has been a considerable inter- 
est in black hole pair production. Inspired by the well understood particle pair 
production of quantum field theory (for example 27 — > e"^ + e~), theorists have 
investigated the corresponding phenomenon for black holes and studied the possi- 
bility that a spacetime with a source of excess energy will quantum tunnel into a 
spacetime containing a pair of black holes. The earliest investigations considered 



pair creation due to background electromagnetic fields |3y, ^ but since 

then have been extended to include pair creation due to cosmological vacuum en- 



ergy [^, [Tj] , cosmic strings ^ ^ , and domain walls p3|, ^ |7^ . In all cases 
the chance of such an event happening has been found to be extremely small, but 
perhaps an equally interesting outcome of the calculations has been the evidence 
that they have provided that black hole entropy does indeed correspond to the 
number of quantum states of the hole. 

In the last part of the thesis I show how the pair creation results can be ex- 
tended to include pairs of rotating black holes, which were not considered in the 
above referenced papers. This is quite an involved process which starts with the 
identification of classical solutions to the Einstein equations that properly describe 
pairs of black holes in the appropriate context. From there instantons are con- 
structed from the classical solutions that will be used to approximate the path 
integrals and it is seen that requirements of regularity restrict the possible physical 
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parameters of the created spacetime. The Brown- York formahsm is used to choose 
the correct action for use in each situation and finally, with all of the preparation 
completed, I calculate and interpret the creation probabilities. 

1.4 Overview 

With these ideas in mind I now outline the rest of the thesis. As its name implies, 
chapter ^ establishes the background for the work that follows. Much of it is a 
review of well-known ideas and results but it will serve to refresh these ideas for 
the reader who is not especially familiar with this area of general relativity and 
establish notation and sign conventions. Since I will be working with a Hamilto- 
nian formulation of general relativity, section explains how a spacetime may be 
foliated into "instants" of time and how a vector field is set up to define the "flow 
of time" from instant to instant. I focus on a finite region of that spacetime and 
discuss its boundaries and the fields on those boundaries in some detail as well as 
give a physical interpretation of the boundaries as being defined by the history of 
a closed two-surface of observers. Extending the spacetime foliation to the timelike 
boundary, I foliate it with closed two-surfaces which define the observers' notion of 
simultaneity. The quasilocal energy will be defined for these surfaces. 

Section reviews the field equations for fields of interest to this thesis. Specifi- 
cally they are gravity, electromagnetism, and a dilaton field where a coupling exists 
between the dilaton and Maxwell fields. First examining these from a covariant 
four-dimensional perspective, I then review how they become constraint and evo- 
lution equations if they are projected into the leaves of the time foliation. I discuss 
how one-half of the Maxwell equations are implied by the assumption that a gauge 
potential exists, a simple fact that will have larger consequences later on, and dis- 
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cuss duality for these three fields. 

With this theoretical stage set, chapter |^ begins the main work of the thesis. 
Starting with a modified Einstein- Hilbert action for gravity, section ^]T| briefiy re- 
views how it variation produces the standard field equations for gravity. The action 
(proposed by Geoff Hayward in ref . |Q ) differs from the classical Einstein-Hilbert 
action in that it is formulated for a finite region of spacetime bounded by a combi- 
nation of spacelike and timelike hypersurfaces and disagrees with the one used by 
Brown and York in that it allows for those boundaries to be non-orthogonal. 



From that action, subsection |3.2.1| derives a Hamiltonian functional defined on 
the slices of the time-foliation of the spacetime. This Hamiltonian differs from 
the Brown- York Hamiltonian in that it does not restrict the time foliation to be 
orthogonal to the timelike boundary. It is noted that even though the functional 
is defined for a finite region of a spatial three-surface, its actual on-shell numerical 
value depends only on the fields at the boundary of that surface and how they 
are evolving in time. It is indifferent to what three-surface it bounds. This means 
that a quasilocal energy derived from the Hamiltonian really is a functional of the 
boundary two-surface only; a fact which is crucial for its proper definition since 
there is no natural way to uniquely associate a spanning three-surface with that 
boundary (or for that matter even guarantee that such a surface exists). Thus, 
this approach differs markedly from that taken by Hawking and Hunter in ref. [ET| 
which required reference terms to remove the dependence of their Hamiltonian on 
the intersection angle between the foliation surfaces and timelike boundary. 

With the proposed Hamiltonian functional in hand, subsection |3.2.2| confirms 
that it really is a properly defined Hamiltonian (this is the first time that this has 
been exphcitly demonstrated) and shows that, as would be expected, the calcu- 
lated variation of the Hamiltonian is in accord with the full variation of the action 
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functional as considered in such papers as [E^, [TO 



Section p.3| presents a definition of quasilocal energy in terms of the quasilocal 
Hamiltonian of the previous section. Its exact form is dependent on the time 
four-vector that determines the evolution of the boundary observers and it is seen 
that if that vector field is a Killing vector field for the induced metric on the 
boundary of the spacetime, then the quasilocal energy is a conserved quantity. A 
special case of the general quasilocal energy where the observers are stationary 
relative to the foliation hypersurfaces and measure proper time is considered and 
geometrical interpretations of that energy are discussed. Transformation laws for 
the quasilocal quantities with respect to boosts of the measuring observers are 
derived and investigated in section |3.4| . These laws are shown to be Lorentz-like 
and a comparison is made with corresponding laws from special relativity. 



Next, section |3T5| defines reference terms for the quasilocal energy of the previous 
sections. These terms are necessary because without them the quasilocal energy 
of a spherical region of flat space is non-zero and actually diverges as the radius 
of such a sphere goes to infinity. Within the quasilocal formalism there is quite a 
lot of freedom to define these reference terms and the choice of a particular one is 
essentially a choice of where to set the zero-level of the energy. I examine three 
choices of reference terms, starting with the original Brown- York term which is 
defined by embedding the instantaneous two-surface of observers into a reference 
three-space. A well recognized problem with this term is that it is not always defined 
and I point out that it also runs into problems for boosted observers in fiat space. 
From there I discuss an alternate proposal involving embedding the instantaneous 



two-surface of observers into a four-dimensional reference space (discussed in |jTO|, |TT 



and from a different perspective in [^). It is more likely to exist than the three- 
dimensional proposal but unfortunately is not uniquely defined. Finally, I briefly 
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comment on the so-called intrinsic reference terms that have recently been inspired 
by the AdS-CFT correspondence. 



The chapter ends with section which investigates the close relationship be- 
tween the quasilocal energy and the thin shell formalism of Israel |5^. I show that 
there is an exact correspondence between the mathematics of the quasilocal en- 
ergy and the thin shell formalism. This means that for a two-surface of observers 
with a specified time evolution in a given spacetime, the quasilocal energy with 
the two-into-four reference terms discussed above is defined if and only if that two- 
surface could be replaced with a thin shell of matter so that outside of the shell the 
spacetime would be unchanged while inside it would be isometric to a part of the 
reference space. The quasilocal energy measured by the observers is exactly equal 
to the total matter stress- energy of the thin shell. This equivalence means that one 
can reinterpret the (modified) Brown- York quasilocal energy from an operational 
point of view. That is, the quasilocal energy "contained" by a two-surface could be 
defined as the matter stress-energy required to reproduce that spacetime outside of 
a matter shell that is isometric to the two-surface and embedded in the reference 
space. 

Though similar in length to the previous one, chapter ^ can be summarized 
quite a bit more quickly since it covers much the same ground except that this 
time Maxwell and dilaton matter fields are included in the mix. It starts in section 



O] with a review of the Lagrangian action whose variation will generate all of 
the field equations. From there, subsection 4.2.1 derives a Hamiltonian functional 
from that Lagrangian action and it is noted that the assumption that a single 
gauge potential exists over the region being studied implies not only two of the 
Maxwell equations, but also that no magnetic charge can exist in that region. 
Subsection [4.2.2| checks that the proposed functional really is a proper Hamiltonian 
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and again compares the variation of the Hamiltonian with previously calculated 
action variations |Tl|. Section |4.3| reviews conserved charges, reference 



terms, transformation laws, and the thin shell correspondence when the matter 
fields are included along with the gravitational field. Finally, in recognition of the 



fact that the formalism as constituted cannot handle magnetic charges, section 
uses duality to define an action and Hamiltonian that can handle those charges, 
though in doing so it loses the ability to deal with electric charges. 

Chapter |^ applies the work of the previous two chapters to several spacetimes 
both to gain insight into the quasilocal energy and to demonstrate its utility. Section 



5TT| is targeted mainly towards the first goal as it examines Schwarzschild and 
Reissner-Nordstrom spacetimes. It starts with static and spherically symmetric 
sets of observers and shows that the quasilocal energies that they measure are 
physically reasonable though not entirely in accord with intuition. Interestingly it 
is seen that the definition of the quasilocal energy derived for a purely gravitational 
field appears to also include contributions from matter fields. The extra terms 
generated by matter fields are gauge dependent and are not directly related to the 
physical configuration of the system but instead seem to give the potential energy 
for the system to exist relative to an (almost) arbitrarily set gauge potential. The 
next section considers radially boosted observers which the original Brown- York 
formalism couldn't easily handle and so demonstrates the nonorthogonal formalism. 
Finally a spherical set of "z-boosted" observers is considered for Schwarzschild space 
and interesting but slightly enigmatic results are obtained. 

Section ^.2| applies the formalism to study naked black holes (first discussed 



by Horowitz and Ross in [|56| , [57[1). These are massive, near-extreme, Maxwell- 
dilaton black holes that are characterized by how different sets of observers feel 
gravitational tidal forces close to the event horizon. Specifically, static observers 
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measure relatively small transverse tidal forces while those who are infalling on ra- 
dial geodesies measure huge (though not divergent) forces. Though at first thought 
this might not seem to be especially surprising, it should be kept in mind that 
equivalent observers near to similar Reissner-Nordstrom black holes all measure 
small tidal forces irrespective of their radial motion. I calculate the quasilocal ener- 
gies measured by corresponding spherical sets of observers and find that the static 
ones measure a very large quasilocal energy while the infalling ones measure it to 
be extremely small. I explain all of these measurements in terms of the geometry 
of the naked spacetimes. 



The final classical application is found in section [5.3| where I apply the formal- 
ism to calculate energy flows during gravitational tidal heating. The prototypical 
example of this in our own solar system is found in Jupiter and its moon lo, where 
the gravitational tidal forces provide the energy that powers lo's volcanism. I suc- 
cessfully reproduce the results of energy flow calculations that in the past have been 



done with Newtonian and stress-energy pseudo-tensor ||80| methods. The calcula- 
tion is cleaner than the pseudo-tensor calculations and has the added advantage of 
providing a simple geometric interpretation of gauge ambiguities in the energy flow. 
Thus, this section can be viewed both as an application of the quasilocal formalism 
and as an additional check on its physical relevance. 

Chapter ^ contains a quantum mechanical application of the formalism as it 
applies it to the pair production of charged and rotating black hole pairs in a cos- 
mological background. It begins with a brief review of the Euclidean path integral 



formulation of quantum gravity in section |6.1| . Section |6.2| examines the classical 
spacetimes that describe pairs of charged and rotating black holes in a cosmological 
background. It starts with the generalized C-metric of Plebanski and Demianski 
which can be interpreted as describing a pair of charged and rotating black holes 
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accelerating away from each other in a cosmological background, and then shows 
that matching the acceleration of the holes to that of the rest of the universe (as 
is demanded by conservation of energy) reduces this metric to the Kerr-Newman- 
deSitter (KNdS) metric. Thus, those will be the class of spacetimes that I aim 
to create and so I examine them in some detail, working out the allowed range of 
their physical parameters and examining limiting cases. Traditionally it has been 
asserted that only spacetimes in full thermodynamic equilibrium can be created 
by quantum tunnelling processes, so I finish off by considering the various KNdS 
spacetimes from this point of view. I show that three limiting cases of the KNdS 
spacetime, the cold limit (which corresponds to a pair of extreme black holes), the 
Nariai limit (where the outer black hole and cosmological horizon become coin- 
cident), and an ultracold limit (the overlap of the cold and Nariai limits) are in 
thermal (but not full thermodynamic) equilibrium clS dbTQ db class of non-extreme 
black holes that are dubbed lukewarm (where the outer black hole and cosmolog- 
ical horizons simply have the same temperature). If the rotation parameter is set 
to zero, then these reduce to equivalent cases considered in extant non-rotating 
calculations. 

Section assembles instanton solutions to mediate the creation of each of the 
classes of KNdS spacetimes that are in thermal equilibrium. These instantons have 
complex metrics, in contrast with the usual ones used to create non-rotating pairs 
of black holes. I show that this is a necessary feature of the instanton metrics if one 
requires that these solutions match onto the corresponding Lorentzian ones along 
a spacelike hypersurface. Finally I construct the actual instantons treating the 
non-degenerate two horizon (lukewarm and Nariai), non-degenerate single horizon 
(cold, ultracold I), and zero horizon (ultracold II) cases separately. I note that only 
thermal (rather than full thermodynamic) equilibrium is required to construct a 
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smooth instanton. 



Section examines the instantons to identify their essential features and then 
uses the Brown- York formahsm to select the actions that will preserve those char- 
acteristics during the quantum tunnelling. Then with that action in hand section 



53| evaluates those actions for the instantons and so finds the probability for a pair 
creation event to occur. It is shown that the probability of pure deSitter space 
tunnelling into a spacetime containing a pair of black holes with opposite spins and 
electric/magnetic charges is proportional to e"^'^'''^ where T^Ai is the sum of the 
areas of the non-degenerate horizons in the created spacetime. This is in accord 
with the corresponding results for non-rotating holes |7l[] as well as an interpre- 
tation of black hole entropy as the logarithm of the number of quantum states of 
the hole. Section shows how the methods that I have used compare with the 
procedures that other people have used. 

Finally, chapter ^ attempts to summarize the results of the thesis, put them 
into some perspective, and discusses future work related to the topics of this thesis. 



Chapter 2 



Set-up 



This chapter sets the stage on which the rest of the thesis will play. In the first 
section I define a quasilocal region of space, show how it may be foliated, and 
define a variety of geometric quantities that will be used extensively in the rest of 
the thesis. The second section reviews the interlocking field equations for gravity, 
electromagnetism, and the dilaton field in their full four-dimensional and projected 
three-dimensional forms. It further discusses electromagnetic potentials and the 
duality inherent in the electromagnetic and dilaton fields. Much of this chapter is 
a review of well known facts but it serves to establish notation and cast some of 
these ideas in a new fight. 



2.1 The geometric background 

Let M be a compact and topologically trivial region of a four-dimensional spacetime 
Ai. It is specified to be the region bounded by two spacelike surfaces Ha and 
(each homeomorphic to M^) and a timelike surface B (homeomorphic to R x S^). 
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Figure 2.1: A three-dimensional schematic of the Lorentzian region M, assorted 
normal vector fields, and typical elements of the foliation. 



Such a region is depicted schematically in figure |2.1j^| . Let M be foliated with a 
set of three-dimensional spacelike surfaces {S^}, labelled by a time coordinate t, 
such that Sfi and St2 are leaves of the foliation. This then induces a corresponding 
foliation of B by spacelike two-surfaces {Qt = n B} each with topology 5^. 
Finally, in association with the foliations define a smooth timelike vector field 
such that T°'dat = 1 and is tangent to B. These conditions are not sufficient to 
uniquely specify T", so there is a certain arbitrariness in its definition. 

Intuitively one can think of B as defining the history of a set of observers and 
each Sf as defining an "instant" in time. Then each fit defines an "instantaneous" 



^In later chapters M and B will sometimes be taken to have more complicated topologies. The 
extensions to those cases will be straightforward so for simplicity I now consider only the simplest 
case. 
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configuration of those observers and in the regular way, I will say that Qti "happens" 
before flt2 if ti < ^2- Further T° can be thought of as the (unnormalized) four- 
velocity of the observers and so B can be viewed as the history of a set of observers 
who started out in the configuration Qti and then evolved through time with 
as their four-velocity. Because isn't normalized the time t doesn't correspond 
to proper time. The freedom in the definition of corresponds to how individual 
observers can evolve differently while leaving their evolution as a set invariant. 

Note that while a foliation surface uniquely specifies a corresponding fij, the 
converse isn't true. Any number of foliations can be defined that are compatible 
with a given foliation. In fact in spite of the way that the foliations have been 
set up in this section, a main goal of this thesis is to show that only the foliation 
of B is important. The foliation of the rest of the spacetime is irrelevant, basically 
because there are no observers in the interior of B to define it. The only observers 
are thought of as residing on the boundary B. 

Up to this point no real use has been made of a metric. Terms like spacclike 
and timelike have been used for clarity but everything could equally well have been 
formulated in terms of a manifold without metric. Now however, I'll introduce a 
signature +2 metric field ga/3 over M. With this metric one can define a (unit 
normalized) forward-pointing timelike vector field -u" normal to the surfaces as 
well as induce a spatial metric field ha/s = Qais + UaUp on those surfaces. Then, one 
can project into its components perpendicular and parallel E^. Namely 

T" ^ Nu^ + V^, (2.1) 

where N and V"' are called the lapse function and shift vector field respectively and 
V^Ua — 0. Conversely one can define the spacetime metric in terms of the spatial 
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metric, lapse, shift, and time vector field by 

gaf3 ^j^ap _ ^^j.a _ya^^j.f3 _y(3y ^2.2) 

Define unit normal vector fields for the various hyper surf aces. Already m° has 
been defined as the future-pointing timelike unit normal vector field to the 
surfaces. Similarly, define as the future-pointing timelike unit normal vector 
field to the surfaces Qt in the timelike hypersurface B. The spacelike outward- 
pointing unit normal vector field to Qt that is perpendicular to (and thus in 
the tangent bundle TEj) is n° and the corresponding normal vector field to Qt 
perpendicular to ■u'^ is n° which is also the outward-pointing unit normal vector 
field to B. 

Next define the scalar field rj = u°'na = —u°'na over B. If rj = everywhere, 
then the foliation surfaces are orthogonal to the boundary B (the case dealt with in 
refs. ||2^, ^]) and the barred vector fields are equal to their unbarred counterparts. 
If ?7 7^ then and may be written in terms of and n° (or vice versa) as, 

n° = + T]u'' and = ^u" + i^n", (2.3) 
A A 

or, 

1 1 
= - r/M" and = - r^n", (2.4) 
A A 

where A^ = V ^nd A may also be written without direct reference to the 

barred normals. To do that first define 

= {V"n^)/N, (2.5) 

which is the three-velocity in the direction of an object with four-velocity T" as 
measured by an observer with four- velocity u°'. Then, 



r] = v^/Jl-vl and X = Jl-vl. (2.6) 
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These quantities then begin to look hke the terms that appear in special relativistic 
Lorentz transforms. This correspondence will be explored in some detail in section 



3.4 



On the surface B one may write, 

= iVw" + V", (2.7) 

where 



N = \N= - [V^n^f and V = a'^pV^ = V - [1/^n^]n" (2.8) 

are respectively labelled the boundary lapse and the boundary shift. This split 
is possible because T° on B has been restricted to lie in the tangent bundle TB. 
Equivalently, B is the history of the observers VLt and T"" is their four-velocity, so 
naturally G TB on B. In any case T'^Ua = 0. 

Next consider the metrics induced on the hypersurfaces by the spacetime metric 
gaf3- Just as haf3 = Qap + WqM/j is the metric induced on the surfaces, the other 
metrics may also be written with respect to the normals. 7q,/3 = Qap — Uafip is the 
metric induced on B and aap = h^p — naUp = •jais + UaUp is the metric induced 
on Qf. Raising one index of these metrics defines projection operators into the 
corresponding surfaces. These have the expected properties: h^^u^ = 'j°'pn'^ = 
a'^^u^ = a'^pTi^ = a'^pu^ = a'^pfi^ = 0, and = 7^)37^^ = 7%, and 

a pa ^ — a ^. 

Let e^^T^ be the four-dimensional Levi-Cevita tensor defined over M. Then, 
fix the orientation of the corresponding Levi-Cevita tensors on S^, 5, and Vtt by 
setting 

e^^^ = M„e°^^^ (2.9) 
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e"/ = r^;3e"^^^ and 

Often, where it won't cause confusion, I drop the subscripts to get a shghtly tidier 
notation. 

Coordinate invariant integrals on M and the various hypersurfaces are defined 
in terms of tensor densities (relative tensors of weight one). Thus, the rest of this 
thesis should really be formulated in terms of tensor densities rather than tensors 
to maximize its aesthetics and remove any appearance of coordinate dependence 



(similar work is formulated in that way in §0[)- -^'-'^ ^^^^ reading however, it 
is more convenient to principally stick with tensors and a coordinate system over 
the region M. The final results will come out the same. 

That said, assume that one can define a coordinate system {r, 9, 0} on S^i such 
that VLt is surface of constant r. Continuously extend it to the other Sf surfaces 
so that {t, r, 6', 0} is a coordinate system over M and is a constant r surface. 
Then, if e"^"^"^, eajj-fS, e^^"^, e^^^, ef"', ef^^, , and e^^ are the Levi-Cevita symbols 
(relative tensors of weights ±1) in the spaces M,J]t,B, and Qt respectively with 
orientations chosen to match those of the corresponding tensors, the determinants 
of the coordinate representations of the metrics are the scalar functions g, h, 7, and 
a that satisfy the relations 

- g^afS-yS = e'^^'"'9aK9(iX9j,,95u (2-10) 

-7ef/37 = eB''7a^7/3A77/^' and 

Combining these relations with equations ( p.9|) it is straightforward to show —g = 
N'^h and -7 = iVV. 
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Define the following extrinsic curvatures. Taking Va as the covariant deriva- 
tive on A4 compatible with gap, the extrinsic curvature of in A4 is Ka/s = 
—K^Ji^qV^us = —\£uhai3-, where is the Lie derivative in the direction m". The 
extrinsic curvature of S in is Q^p — ~l'^al^ i^s while the extrinsic curvature 
of in is = —(T~'^(T^ijV^ns. Contracting each of these with the appropriate 
metric define K = h°'^Kaf3, © = ^°'^0a/3, and k = a°'^kap. The addition of an 
overbar to any quantity will indicate that it is defined with respect to and/or 

rather than and - for example, k = — u^^Van/j. 

Further, define the following intrinsic quantities over M. and E^. In Al, the 
Ricci tensor, Ricci scalar, and Einstein tensor are T^-q,/?, TZ, and Gai3 respectively. 
Da is the covariant derivative on compatible with hap-, and is the covariant 
derivative on VLt compatible with cTq,^. i?Q,^ and R are the Ricci tensor and scalar 
intrinsic to the E^ hypersurfaces. The sign convention for the Riemann tensor is 
such that VaV^a;-y — V fjVa^-y — T^ap-y^^s ^ covariant vector field cua- 

Finally, from the preceding it is clear that tensors defined over M will usually 
be written with Greek indices. However, in cases where these tensors can defined 
entirely in the tangent and cotangent bundles of the surfaces E^ they will often be 
written with Latin indices instead. 

2.2 Field equations 

With the stage set, I now consider the fields that are the players in this spacetime. 
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2.2.1 The 4D field equations 

Consider spacetimes with a cosmological constant A, a massless scalar field (f) (the 
dilaton), and a Maxwell field In units where c, h, and G are unity, the field 
equations are: 

^e"^^^V^F^5 = 0, (2.11) 
V/3(e-2"'^F"^) = 0, (2.12) 
V"V„0 + ^ae-^'^'^F^^F"^ = 0, and (2.13) 
Gap + Agap - kTo,^ = 0, (2.14) 

where a is the coupling constant between the scalar and Maxwell fields, k = Sir (it 
would take on a less trivial value in other systems of units) and 

Tap ^ i- (^[V„0][V;30] - ^[W^][V,<P]gap + e-2'^<^[F„,F/ - ^gapF.sF^']^ 

(2.15) 

is the stress-energy tensor associated with the matter. From the field equations it 
is clear that there are no EM or dilaton charges or currents in the region under 
consideration. I work with the sign convention that A is positive for deSitter space. 

The first equation implies that, at least locally, it is possible to define a vector 
potential A^ such that Fa/s = daAp — dpAa- Conversely if one takes the vector 
potential as a pre-existing field and Fq,^ as a quantity derived from it, then equation 
(|2.11| ) automatically holds (it simply expresses the identity d{dA) = for any 
differential form A). This is a common viewpoint, and from chapter |^ onwards, Aa 
will be taken as the primary field and so equation ( |2.11| ) will be that identity. The 



other equations of motion will then be derived from the variational principle. 
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2.2.2 The 3D field equations 

Much of this thesis works with Hamiltonians and as such it will be useful to know 
how these field equations project down into the three-dimensional spatial hypersur- 
faces Sj. First define (dilaton modified) three-dimensional electric and magnetic 
vector fields in the usual way. That is, 

Ea = e-2"'^F„0?/ and (2.16) 
= ^u5e'^f'~'Fp^. (2.17) 

Conversely F^js and may be rewritten in terms of the electric and magnetic field 
three- vectors as 

Fan = e''^t>^u^Ep-upE^) + u'esap^B\ (2.18) 

While these are the most commonly seen definitions of the electric and magnetic 
fields, for much of the following it will be more convenient to work with the related 
vector densities on the hypersurfaces, defined by 

= and B'' = ^S". (2.19) 

K K 

With respect to these vector field densities, the Maxwell equations can be pro- 
jected into components perpendicular and parallel to the hypersurfaces using 
the identities 

^^u^e'^f'^'VpF^s = DfsB^, (2.20) 

K 

^^^^a/J75v^F,5 = h>;£u[NB^]-h;u^e'^f'^'D,[Ne'^%], (2.21) 

^ii«V;3(e-2«'^F"^) = DpSf", and (2.22) 

^^^/i^V^(e-'«'^F"^) = h>'^£u[NS^] + h>'^u^e''^^^D^[Ne-^'"^Bs]. {2.23) 
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If these equations were written with respect to vector fields instead of vector field 
densities they would include unaesthetic extrinsic curvature terms. 

Then, for time derivatives defined as the Lie derivative with respect to the time 
vector T° rather than and using Latin indices to emphasize that all quantities are 
defined exclusively in the hypersurface, the three-dimensional Maxwell equations 
are 

DbB'' = 0, (2.24) 
h^p£TB^ = e^''^D,[Ne-^''^£d]+ £vB\ (2.25) 

DbS^ = 0, and (2.26) 
/iJ£t^^ = -e''"^D,[Ne^'"^Bd\+ £v£'', (2.27) 

where e^'^'^ = Ua€°'^'^'^ as was defined in the previous section. 

Next consider the dilaton equation ( p.l3| ). It takes its simplest three-dimensional 
form written in terms of the scalar density 

p = -^£u(t). 2.28 

K 



Then 



2^ ( V-V„0 + ^e-^^-^F^^F^^) (2.29) 
K \ 2 / 

or equivalently the time rate of change of p is 

£tP = 'b^DmoM + (e-^'"^B'Bb - e^'^'^S'Sh) + £vp, (2.30) 

where again I've used Latin indices to emphasize the three-dimensional nature of 
the equation. 
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Finally consider the projections of Einstein's equation ( p.l4| ). There are three: 
time-time, time-space, and space-space. Again it is most convenient to work with 
a tensor density, namely 

P°/3 ^ ^Kh'^^ _ , (2.31) 

which is contracted as P = hapP"^- Ignoring the matter terms, the equations 
project as 

\/h 

n = - — {G^^ + Ag^p)u''uf (2.32) 



Hp = ^hliG^s + AQjsW (2.33) 



K 

-2DpPf, and 



^ !:!^h''ihf'\G^s + A9,6) (2.34) 
= h''^hl£u{NP^') + ^^(R'^-\Rh''^ + AhA (2.35) 

{D'^Df'N - h-^D"'D,N) P'''P,s - -pA h"^ 

2k vh\ 2 / 

-J_4^ ( P7{"p 13) _ lpp"/3^ 

V " 2 
Next, the matter terms project as 

= ^[D^^][D,<P] + ^p' + ^{e"^^£'^£^ + e-"^^B"B^), (2.36) 
7^;, = -^fhhlT^iU^ = pDfi(l)^ ^^efi^^E^'B^ and (2.37) 



2v^ 



(2.38) 
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Nk f 



4v^ 

Again the right-hand sides of the equations are composed entirely of three-surface 
terms and so could be written with Latin indices. 

Combining the projections, the Einstein equations may be rewritten as two 
constraints and a time evolution equation. They are 

rr = H + n' = o, (2.39) 
H];; ^ nb + K = 0, and (2.40) 

h''^h\£TP^^ = -— {N^^^G"^ + Ah"^ - [D^D^N - h^^D^D'N] ) (2.41) 



(^P^'^Pcd - ]^P^]h''^ - 4[P"("P/) - ipP"^]^ 



Nk 

+£vP"'^ -n""^ 



where ^^^Gab = Rab - {l/2)Rhab. 



2.2.3 3D electromagnetic potentials 



As noted at the end of section p.2.1| , either by assumption or equation (|2.24|) there 



(locally) exists an electromagnetic vector potential Aa such that Fa/s = daAjj — 
d/3Aa- Here I examine how that four-dimensional potential breaks up into the 
regular Coulomb potential and a three-dimensional vector potential. To wit, define 

$ = -Aau" and i„ = h^Ap. (2.42) 

Then, 

E0 = -e-^'^'^ (J^DpiN<l>) + hl£^A^^ and (2.43) 

5/3 = Uae'^'^^^D^As. (2.44) 
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Strictly on the hypersurface and with respect Lie derivatives in the T° direction, 
these become a time evolution equation for the three-vector potential and a defini- 
tion of i3" in terms of the vector potential respectively. That is, 

Nk 

KI£tA^ = -^e^'''^£b + £vAa- Da[N^], and (2.45) 
B' = '^e^'^'DM (2.46) 



2.2.4 Duality 

Defining the dual *-Fq,/3 = |e~^"'^e^^'^'^F^5 of Fap, the four-dimensional field equa- 
tions ( |2.11| , |2.12| , |2.13| , and p.l4| ) may be written as 



V/3(e2'^'^*F"^) = 0, (2.47) 

—e'^f^'Vfs^F^s = 0, (2.48) 

V"V„0 - ^ae^""^ ^Fo,p ^F"^ = 0, and (2.49) 

G,/? + A(7«/3 - SvrT,^ = 0, (2.50) 

respectively where this time the stress energy is 

Tap = ^ (^[V„0][V^s</.] - l[V^my4>]9ap + e^'^'t'l-^Fa, - ^g^p ^F,s • 

(2.51) 

Thus, the four-dimensional field equations as a set are invariant under the full 
duality transformation (0 —ip^Fap — > -kFap). Note however that this time it 
is equation ( p.48| ) that is equivalent to a statement that there (locally) exists a 
potential one form. It implies that there exists an A* such that -kFap = daA^^—dpA*^. 

Next, in terms of -kFap the electric and magnetic vector fields may be written 

as 

E^ = ^u^e^'^^'^i.F^s and (2.52) 
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= -e^-l'i^F^f^u^, (2.53) 

or equivalently 

= -e-^'"('{uaBp-upB^)+u'esapjE\ (2.54) 

Thus, with respect to these fields, the duahty transform Fa/s ^Fap becomes Ea — > 
—Ba and Ba — > Ea- Combining these two with — > — and the corresponding 
p — > — p, it is a simple exercise to check that the set of three-dimensional field 
equations for gravity, electromagnetism, and the dilaton field are also unchanged. 

Decompose in the same way as A^. That is define 

= -A>° and i* - h^A*p. (2.55) 

In terms of these potentials the electric and magnetic fields may be written as 

El^ = M„e"^^^D^ij and (2.56) 
Ba = e"''t'(J^Da{N^*) + £uAiy (2.57) 

which alternately may be viewed as a definition of in terms of A* and a time 
evolution equation for A*. Namely, 

Nk 

hl£TA*^ = ^e-^"'^ I3b+ £vAl- Da[N^], a.nd (2.58) 
'^"^e'^^'DJl (2.59) 



K 



Chapter 3 

A quasilocal Hamiltonian for 
gravity 



This chapter presents a quasilocal Hamiltonian formulation of gravity. I start in 

use 



section ^]T] with the standard Einstein-Hilbert action and then in section 
a temporal foliation and Legendre transform to define a Hamiltonian functional 
for general relativity over a finite region of a spacelike surface and with respect 
to an arbitrary time evolution. This Hamiltonian and its derivation are similar 
to the well known ADM formalism though here the analysis is conducted for 
a finite and bounded region of a larger spacetime. I confirm that the proposed 
functional correctly generates the equations of motion and show how its boundary 
terms can be used to define quasilocal quantities such as mass, energy, and angular 
momentum. These boundary terms depend only on the values of the fields at the 
boundaries. 



With these quasilocal concepts defined, in section Oil consider conserved quan- 



tities and calculate the time rates of change of their non-conserved equivalents. Next 



33 



CHAPTER 3. A QUASILOCAL HAMILTONIAN FOR GRAVITY 



34 



section examines how they transform with respect to boosts of the time evolution 
vector field and shows that those transformation laws are pleasingly Lorentz-like. 



From there, in section I survey a variety of proposals about how to define the 
zero of the action and Hamiltonian and discuss the specific instances in which each 



is useful. Finally, section examines the close relationship between the quasilocal 
formalism and the thin shell work of Israel. This relationship makes it possible to 
recast the definition of quasilocal energy from an operational point of view and at 
the same time use quasilocal insights to shed light on the physics of thin shells. 

Most of this work was published in [10] and parts of and 0. This thesis 
however is the first place where the variation of the Hamiltonian has been explicitly 
calculated. 



3.1 The gravitational Lagrangian 

Given M G Ai as described in the previous chapter and allowing for the inclusion 
of a cosmological constant, the appropriate action for general relativity is 

I = — [ d^x^^{JZ-2K) + - [ d^x^K-- [ rf^xv^e (3.1) 
2k J m Jt. J b 

+ — / d'^xy/a sinh^^ 1] + I, 
1^ Jn 

where = — J^^ , = J^^ — J^^ , and, choosing a system of units where c, h, 
and G are unity, k = 8n. The sinh~^ t] term is added so that the variation of the 
action will still be well defined if the boundaries are not orthogonal to each other at 



their intersection. It was first discussed in |5^. / is any functional of the boundary 
metrics on dM. 

To see that this is indeed the correct action, take its variation with respect to 
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the metric gap- The result is |^ 



" = 2- 



p / d'x./^{G^p + Aj„^)ij"» (3,2) 

is J B 

+ y d^a; |^-sinh"^(r7)(5v/a^ +5/, 

where P"^ = ^ i^Kh'^^ — i^"^) is the same tensor density defined by equation 
( |2.31| ) in the previous chapter, and vr"^ = —^^^ (67"^ — B"^) is an equivalent 
tensor density defined by the surface B. 

For variations that leave the boundary metrics hap and 7q,/3 fixed, the boundary 
terms, and 6I_ vanish. Then 5/ = if and only if the Einstein's equations hold over 
all of M. Thus with these terms fixed, the variation of / is properly defined and 
generates general relativity as asserted. 



3.2 The gravitational Hamiltonian 
3.2.1 Form of the Hamiltonian 

With this quasilocal action in hand, it is a relatively simple matter to obtain the 
corresponding quasilocal Hamiltonian. The process is to decompose the action with 
respect to the foliation and then identify the Hamiltonian and momentum terms. 



Details of the calculation may be found in appendix [A.l| , but here I'll just present 



the results. Breaking it up with respect to the foliation the action may be written 



1-1= [dt\[ d^x{P'^'^£TKp)+ f d'^x{P^£T^/^) - Ht\ (3.3) 
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where retains its meaning from the previous chapter, = ^ sinh ^ r], and 
Ht= [ d^x[Nn + V"na]+ [ (fx^/^{N£-V"ja), (3.4) 

where 7i and Ti^ are the matter free versions of the Einstein constraint equations 
(|2.39| ) and ( p.40D . N and V'^ are the boundary lapse and shift defined following 



equation (p7|), while 



1 2 . 1- 

e = —k + ri—=P"'nanb = —k, and (3.5) 

Ja = — '^aa^P^'^ns - -a^dfsT] = --a^tCW pfi^. (3.6) 
yh 1^ 1^ 

Shortly e and ja will be identified as related to energy and angular momentum 

respectively but for now simply note that despite the initial appearance of these 

terms, their second versions show that they are really defined with respect m", n", 

and and as such are defined with respect to the foliation Vtt of the boundary 

and the normals and n" rather than the foliation and its normals and n". 

To motivate the definition of the Hamiltonian, recall that in elementary classical 
mechanics with one degree of freedom, the action / of a path q = r{t) taken by 
a particle is given by / = L{r{t))dt where L is the Lagrangian function and 
the integral is over the path. This is related to the Hamiltonian H by the relation 
L = pq—H, where q is the variable giving the configuration of the system andp = ^ 
is the momentum conjugate to q. Extending this analysis to gravitational fields 
[p2| and referring back to equation (|3.3| ), hap may be identified as a configuration 



variable on the spatial surfaces and recognized as its conjugate momenta. 
Further ^/a is seen to be a configuration variable on Qf (albeit one that is not 
independent of hajj) and P^ is its conjugate momentum. Finally perform an 
effective Legendre transform by identifying quantity Ht as the required quasilocal 
Hamiltonian. 
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The bulk term integrand of Ht is NTi. + V^Tia where Ti, and Tia are the Einstein 



constraint equations (|2.39|) and ( p. 401) . Then, in the standard Hamiltonian way, 
the lapse and shift are identified as Lagrange multipliers rather than configuration 
variables. Further since the constraints will be zero for solutions to the Einstein 
equations, the actual numerical value of Ht will be a functional of the boundary Qt 
and its normals and ra" only (and recall that these normals are fixed by Qt and 
without reference to S^). Therefore the evaluation of Hf doesn't require any 
knowledge of the surface apart from the fact that is has a boundary Qf This 
indifference to the bulk will be considered further in section |3.3| . 

By contrast, the nonorthogonal Hamiltonian proposed by Hawking and Hunter 



in reference ||5T| focused on the foliation surfaces and normals and n", which 
meant that their Hamiltonian was explicitly dependent on the intersection angle 
parameter t]. They had to resort to a clever choice of the reference term J to remove 
this dependence. 



3.2.2 Variation of the Hamiltonian 

This subsection checks that Ht really does encode the correct equations of motion 
for gravity. To do this, consider Ht as a functional of the surface St, its boundary Qt, 
the normal ra" to that boundary, the fields hab and ^/a along with their conjugate 
momenta P°'^ and P and the Lagrange multipliers and V"". In the usual 
Hamiltonian way the conjugate momenta are taken to be entirely independent of 
hab and ^/o^. Further e and ja are considered to have their first meanings from the 
definitions ( |3.5| , ^.6| ), and N, F", A, and r] are defined entirely with respect to V"', 
N, and Ua as expressed by equations ( p.5[ - p.8|) . Thus, Ht is a functional on the 



three-space rather in the four- dimensional spacetime M. 
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Now vary Ht with respect to the three-metric ha/s, the conjugate momentum 
P"^, and the lapse and shift and V". Because ^/o^ and P ^ are functions of A^, 
V"', and hai3, these two secondary quantities are automatically varied as well. This 
is an important calculation but its details are not really pertinent to the main ideas 
of the thesis, so I banish them to appendix |A.2| and go straight to the final result. 
The total variation of Ht is 

SHt = [ d'x{nSN + naSV''+[Kb]TSP'"'-[P''']TSKb) (3.7) 

+ [ Sx^ {e6N - jJV^ - {N/2)r'6aab) 
Jnt 

+ [ d^xV^{[V^^6P^-[P^]^6V^, 
Jnt 

where 7i and Tia retain their previous values, while 

AkN 1 

[hat]T = -J=-[Pab--Phab]+2D^aVb), (3.8) 

[P"'']^ = -— {Nf^G"'' + Ah"''] - [D'^D^'N - K'^'D^D^N]) (3.9) 

2 A€ 



Nk 



(^P'^^Pcd - - A[P''^''Pj'^ - ipp'^''] 



^^h 

= ^ (A;"'' -[k- n'^aala"^) - -^rja^^a'^aP"^ (3.10) 
[v^^ = -^(^N^P'^'naUb + N'^k-^dbV'^ , (3.11) 



and [-P^jj. is an undetermined function over Qf. I'll interpret s'^'' in section 
and show that it is actually independent of rj, despite initial appearances. If r/ = 
it becomes the stress tensor s"''' = {1/ K){k°''' — [k — n'^ad\<y"'') considered by Brown 
and York. 
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The Hamiltonian equations of motion can now be obtained by calculating the full 
variation of the action (|3.3| ) (treating the momenta as independent variables) and 
solving 61 = 0. Using the preceding result, only a little work using the fundamental 
theorem of calculusQ to move total time derivatives to the spacelike boundaries is 
required to show that 

51-51 = [ d^xP''^5hab+ I (fxP^5^ (3.12) 

- j j (fx{mN + na5V''} 

+ jdtj d^x{i£Thab - [habjr) SP"' - [£tP''' - [P"'] ^) 5hab] 
+ [dt [ d'x{{£TV^~ [,/^]^)5P^-{£tP^-[P^]t)SV^} 

- j dt d^x^ (^e5N - jJV'' - js''''5aab^ , 

where again = - /^^ and = J^^ - J^^. 

If metrics 7^/3 (equivalently aab, N, V'^) and hab are held constant on the timelike 
and spacelike boundaries respectively, then 5/ = and solving 51 = while allowing 
for general variations in the bulk gives the following set of equations. 

(3.13) 
(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 



n 


= 0, 




= 0, 


£tP''^ 


= 


£ rhab 




£tVo' 


= [v^y; and 


£ tP^ 





iThatis/;fd/ = /(f2)-/(ii). 
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Now ( p.l3| ), ( p.l4| ), and ( |3.15| ) are the (matter free) projected Einstein equations 
(|2.39| ), (|2.40| ), and (|2.41|) respectively, so the Hamiltonian has recovered those cor- 
rectly. At the same time, (|3.16|) and ( p.l7|) give the correct expressions for the Lie 



derivatives of hat and y/a as compared to direct geometric calculation. 

Finally, equation ( p.l8|) correctly expresses the fact that the time rate of change 
of is undetermined by any of the other quantities - a fact that is to be expected 
since in the Lagrangian formulation = ^smh~^r], where 77 = v^/ ^Jl — f ^ 
and fh = {y"na)/N (equation ( |2.5| ) and the surrounding discussion). The lapse 
and shift are Lagrange multipliers whose time evolution is not determined by the 
equation of motion. Therefore the evolution of is similarly undetermined. 
Intuitively, this is to be expected since vy- quantifies the "radial" evolution of VLt or 
equivalently the radial "shape" of B. The "shape" of B is chosen arbitrarily so one 
would certainly not expect P^ to be determined by the field equations. 

Thus, Hf is a proper quasilocal Hamiltonian as supposed. 
Comparison with the Lagrangian approach 

Before moving on to the next section I will compare the variation of the Hamiltonian 
with the variation of the action (and its decomposition with respect to the time 
foliation) as considered in refs. |TU|. Specifically I compare with [|T^ where we 
allowed for a non-orthogonal intersection of with the boundary B. Reference 



p2| deals with the special case where = 0. In those papers, the variation of the 
action (equation ( |3.2| )) was decomposed according to the foliation, the key result 
being that 

vr°^(57a/3 = -v^ {e5N - jJV'^) + ^^^r^da^p, (3.19) 
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where all quantities retain their earlier definition though with the recognition that 
P"^ = \^/{2K){Kh°'^ — K""^) (as opposed to Hamiltonian calculations which treat 
it as an independent variable). Then, equation ( p.2|) becomes 



6I-6L = ^[ d^x^{Go^p + Kgo^p)5g''^ (3.20) 

+ / (fx {P'^^dKp) + / (fx{P^5^) 
Jt. Jn 



where again is recognized as a function of the metric (7q,/3, its compatible covari- 
ant derivative Vq,, and the embedding of in A^. With this viewpoint equations 
( 3.1(: , I!T? , j.l8D are automatically satisfied and so the Lagrangian and Hamiltonian 
treatments are equivalent - as of course they should be. 



3.3 Energy and Kt 

In classical mechanics the value of the Hamiltonian is identified with the energy 
of the system under consideration and so by analogy Brown and York identified 
(the hypersurface orthogonal version of) Ht with the mass/energy contained by 
the surface fi^Q. A key point in favour of this identification is the fact that for an 
asymptotically fiat spacetime, iff is numerically equivalent to the ADM and Bondi 
masses in the appropriate limits (as shown in ||22[ and |jl6[ respectively). 



Tentatively making this association, recall that the energy of a mechanical sys- 
tem is conserved if and only if it is isolated from all outside infiuences. Now, a 
finite gravitational system can be considered to be isolated if the metric 7q,/3 of B 
has a timelike Killing vector field and there is no fiow of matter across B (that is 
^Or, more properly Ht is the energy associated with the surface fit as discussed in section |l.2|. 
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Tapn'^u^ = 0). That said consider the time rate of change of Ht. Because I am 
considering a pure gravitational field here, there will no matter flows across B. 

This is a surprisingly easy calculation to do because the Hamiltonian variation 
calculation can be easily recycled to do all of the work. Equation ( p.7|) showed how 
the Hamiltonian is changed by general first order variations of the metrics and their 
conjugate momenta. Of course, during that calculation there was no assumption 
made that the metrics and momenta satisfied the Einstein equations - the point of 
that calculation was to derive those equations. However, that said, the mechanics 
of the calculation equally well hold for variations that do satisfy the equations 
of motion. In particular consider a region M of spacetime with metric gaf3 that 
is a solution to the Einstein equations. Then, evaluate Ht over a spatial three- 
surface St with two-boundary Qf Lie-drag that surface forward by an infinitesimal 
amount of coordinate time, in which case 6hab = {£Thab)St, ^Vo' = {£ TV^)St, 
SP'^b = {£TP"'b^5t, and 5P^ = {£TP^)St. Combining these substitutions with 
the fact that the Einstein equations are satisfied on S^, the first and third lines of 
equation (|3^ ) go to zero and leave behind the time rate of change of Ht 

£THt= hm^^ = d^xV^ !^e£TN - j^£tV'' ~ ^s^''£Tao^pY (3.21) 

or alternatively using eq. ( p.l9|) 

£THt = - [ d\{'K''^£Tlap], (3.22) 

which is often the most convenient form for explicit calculation^. Note that just 
as the Hamiltonian itself depended only on the foliation of the boundary and its 
associated normals, so does its time rate of change. What is happening in the bulk 
is irrelevant. 

■^An alternate calculation of this result which also allows for matter flows can be found in 
reference ITI. 
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£THt is zero if £tN = 0, £tV'^ = 0, and £TO'ab = or equivalently if 
£T'yaf3 = 0. This is almost the definition of an isolated gravitational system (in the 
absence of matter flows) that was proposed a couple of paragraphs back except that 
there only the existence of a Killing vector was required rather than demanding that 
T" be that vector. £THt might not be zero even if the Killing vector exists. As 
an example consider Schwarzschild space with i? as a surface of constant r. Then 
the lapse function N can be chosen so that Ht is not a constant even though the 
Killing vector exists. See for a further discussion of this point. 



Viewing Ht as a mass, it is useful to think of B as the history /future of a set 
of observers as was discussed in section Then, as noted there, the foliation 
Vlt defines the "instants" of time agreed on by those observers and T" defines 
their four-velocity. Thus, the quasilocal Hamiltonian can be thought of as a kind 
of Gauss's law for mass, in the sense that it defines the mass contained in the 
bulk without making any reference to what is actually happening there, just as the 
electromagnetic Gauss's law defines the electric charge contained by a surface based 
entirely on measurements made at that surface. It then makes sense that the time 
rate of change of the Hamiltonian should also depend only on what happens at the 
surface since the only way energy can get in or out of the bulk should be through 
that surface. 

If the boundary is made up of observers, it is reasonable that their notion of the 
energy contained by the surface should not depend on the bulk foliation. There are 
no observers in the bulk and so there is no natural way for the boundary observers 
to globally extend their notion of simultaneity into that bulk. Thinking empirically 
these observers would say that the foliation of the bulk is a fiction invented by 
theorists that has no external reality. Locally the natural foliation for the observers 
to consider is Si - the foliation that is orthogonal to their four velocity T". That 
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is, they would view as the natural normal to the surface If that foliation 
could be extended throughout M then the numerical value of the Brown- York 
Hamiltonian is identical to its generalized form considered here. However, whether 
that extension exists or not is irrelevant from the point of view adopted in this 
thesis. 

A special case of this Hamiltonian definition of energy is when = u°'. That is, 
the observers are evolved by the timelike boundary unit normal to Qt and measure 
proper time. Such observers measure an energy of 



and it is this energy that is used in applications of the Brown- York energy to 
thermodynamics (see for example [jlSl , p| , |27| , Because of this identification 

y/ae is usually called the energy density. For r] = (that is the foliation Sj is 
orthogonal to -B) it reduces to the Br own- York energy density but in any case it 
will be referred to as the quasilocal energy or QLE. 

This measure of quasilocal energy has a nice geometrical interpretation and 

that is the reason for the subscript in Eceo- Specifically, ^/o'e = {^/a / K)k = 

— {l/^^)£n^/o^ and so measures how the surface area of Qt changes if it is translated 

"radially" outwards in the direction n". Similarly, the surface area measures how 

the volume of the region contained by Qt changes if one "radially" translates the 

surface outwards. Now, of course the volume in Qf depends on the curvature of the 

space contained therein so it is not unreasonable that its "second radial derivative" 

might tell one something about the gravitational energy^. 

^ Contained volume here is a very hand-wavy notion since as has been already emphasized, 
the volume contained by fit is really very dependent on the behaviour of the foliation and the 
derivatives themselves depend on how one chooses the radial normals . However this is a useful 




(3.23) 
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Other symmetries of the boundary B correspond to other conserved charges. I 
won't consider the details here, but it is not hard to show (see for example |2^ ) 
that if B admits an angular (spatial) Killing vector C,'^ then J = J^^ d'^Xy/a^'^ja 
is the charge corresponding to this angular symmetry. In the appropriate limit it 
agrees with the ADM definition of angular momentum |^ at spatial infinity, and 
so -/aja is usually called the angular momentum density. Interestingly can also 
be identified with the connection on the normal bundle to fit. A good discussion 
of this and its implications can be found in reference . 



3.4 Transformation laws 

Having defined the quasilocal energy it is natural to ask what is the relationship 
between the quasilocal quantities e, J^, s"^ as they are seen by different sets of 
observers moving with different four-velocities. 

Consider two sets of observers who instantaneously coincide on the surface VLt. 
Let the evolution of the first set of observers be guided by the Sf forward-pointing 
timelike unit normal vector while the second set is evolved by the time vector 
field T°. Henceforth I'll refer to the observers as the "unboosted" observers 
while the T° set will be the "boosted" observers. 

The evolution of the unboosted observers is orthogonal to the foliation so they 
view and as the unit normals to Vtt. Meanwhile the boosted observers regard 

and as the unit normals. The unboosted observers measure the radial velocity 
way to think about the energy intuitively, and in any case corresponds to the volume changes 
with respect to the natural local orthogonal foliation Sj. 
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of the boosted observers as 



= ^LJhl = LJhl^ (3.24) 



and their Vlt tangential velocity as 

Vr, = = . 



(3.25) 



Then, recall equation (|2.6|) which showed that r] = u°'na and A = ^ — - can be 
rewritten as 

rj = 7f|- and v = 7; 

A 

where 7 = 1/ a/1 — is the Lorentz factor for radial v\-. With this substitution, 
equations ( p.3|) can be written as 

= 7(n° + wkm°) and = 7(m" + f^n"). (3.26) 
The extrinsic curvature of Qt with respect to the timelike is defined as 



ki(^ = -crlaf^VyUs = --alafi£uO-^5 (3.27) 



which can be contracted to = cr^^kj^p- The rate of change of n" in the direction 
it points is ai = n^V/sria. The choice of the | superscript is meant to suggest an 
interchange of u" and in these quantities (as compared to the same expression 
without the superscript) and as usual the addition of a bar means that they are 
to be calculated with respect to u°' and rather than m° and n". The quasilocal 
quantities with and n° interchanged are: 

el = -fcl, (3.28) 
jI = --a^V^us, and (3.29) 
^ l-(^ki^-[kl-u'^al]a^f,). (3.30) 
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Note that jl = ~jc 



Some of these quantities were first used in in the context of defining quan- 
tities that are invariant with respect to boosts. The simplest example of such an 
invariant is e"^ —e^"^ which is analogous to m^c^ = E"^ —p'^c^, an invariant for a par- 
ticle with energy E and momentum p in special relativity. This suggests that be 
viewed as a momentum flux through the surface flf Support for this interpretation 
comes from noting that 

y^^I = ^^^(y'^P £^(y = (3.31) 

That is, is zero if and only if the observers don't see the area of the surface they 
inhabit to be changing. However, this means that a sphere of observers moving 
at constant radial speed in flat space will measure a momentum flux so this isn't 
entirely in accord with intuition. Of course without reference terms such observers 
will also measure a non-zero quasilocal energy so this is not entirely unexpected. 
A more complete discussion of the identification of with momentum may be 
found in which also develops a notion of quasilocal energy from the invariant 



— which is closely related to the one considered here. 

A series of straightforward calculations leads to expressions for the quasilocal 
quantities seen by the boosted observers in terms of quantities measured by the 
observers. These transformation laws are 

e = --(T"^V„n;3 (3.32) 

I 



7(£: + v^e^ 



]a = --alu'Vpn^ (3.33) 

Kj 

K 
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= ja cr^df^v^, and 

Sa/3 = - ikap n^M^^c^p) (3.34) 

Kj 

= -rSal3 + VSal3+ -(^afSi^uV 

T 7^ 
= 7(sa/3 + ^^hS^/s) H 

The reader will recall that a quantity has already appeared in equations ( p.lO|) 
and ( p.l9|) . Short calculations show that s^^ = —^aa^yO'^sP'^^, and £7-77 — £ yrj = 
N £^7] so these two quantities are the same. Further, the first line of equation ( |3.34| ) 
shows that is independent of t] and the bulk foliation S^. 

If the unboosted observers and their time slice Sf are static in the sense that 

pan 

a^p = and P'^'^rianp = 0, and the boosted observers have a constant radial 
velocity over VLt (ie. = constant and -Oq, = 0), then these laws greatly simplify. 
Specifically, 

e = 76, (3.35) 
Ja = ja, and (3.36) 



^afi ^ + (3.37) 



K 

So, in this case the energy density transforms as might be expected from special 
relativity. The angular momentum density is an invariant which isn't too surprising 
considering that it is perpendicular to the direction of the boost. However, the stress 
tensor has a somewhat more complicated transformation law that is dependent on 
the perpendicular component of the acceleration of the boosted observers. Breaking 
it up into pressure (ie. trace) p = s'^^aap and shear (ie. traceless) r/"^ = s"^ — 
{l/2)pa"^ parts a little simplification results. Namely r^"^ = 777°"^ and so it loses 
its acceleration dependence. However, p = •yp + / k) £ nVi- and the dependence 
remains there. 
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3.5 The reference term 

If one calculates the quasilocal energy contained by a spherical shell in Minkowski 
spacetime it is immediate that the reference term / cannot be neglected as I have 
been doing up to now. To see this, let the sphere have radius R. Then ^/ae = 
—R/IAtt) and so the QLE is —R. This is manifestly not zero, and what is more it 
actually diverges as i? — ^ oo, which are not properties that one would expect flat 
space to have! Thus, in anticipation of the upcoming calculations in chapter ^ it 
is time to consider / 7^ 0. 

It has already been seen that 5/ = for variations that leave the boundary 
metrics unchanged, and so its exact form does not affect the equations of motion as 
derived by the Lagrangian or Hamiltonian principles. However, it is equally clear 
that it does determine the zero of the numerical value of the action and therefore 
the zero of the evaluated Hamiltonian and all quantities derived from it as well. 
In this section I'll consider some specific choices of / and discuss the merits and 
problems of each. 

Setting / = 

First, consider when / = might be of some use. As pointed out this means that 
Hf will have non-zero values for finite regions of flat space and it is not hard to see 
that a similar problem arises for the action itself. However it does have the strong 
argument of simplicity in its favour, so it is worthwhile to consider circumstances 
where it might be of use. 

If one wishes to compare the energies contained by two almost identical surfaces, 
each embedded in the same space, then this may be a reasonable choice as any 
reference terms will (at least approximately) cancel each other out. In fact, if one 
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uses reference terms defined entirely with respect to the two-boundary metric aab, 
instead of the full three-boundary metric 70/3, and considers how the associated 
energy changes as the surface is smoothly deformed, then the terms do cancel 
exactly so one doesn't need to worry about them. This is essentially because such 
terms can have no dependence on the time rate of change of aab- Examples of this 
class of reference terms are the 2D into 3D embedding reference terms considered 
next and the 2D intrinsic reference terms considered after that. 

Given these facts, J = is often used when one is doing thermodynamics ||2l| , ^ 



In section ^]3| where I examine energy flows through a quasilocal surface, I'll assume 
7 = 0. My main motivation was to simplify an already complicated calculation, 
but as I have just pointed out, for a smoothly deforming surface Qt a wide range 
of reference terms reduce to exactly this case. 

Embedding ilt in a 3D reference space 



In |2^ Brown and York suggested that one should embed the two-surface (fit, aap) 
into a three-dimensional reference space such that its intrinsic geometry is un- 
changed. One can then define 



/ d^x^/aNe, (3.38) 
Jb 



where e is calculated calculated for Qt embedded in the reference space (usually 
taken as with metric 6ab = diag[l, 1,1]). I omit the term since it fundamentally 
depends on how is embedded in M rather than on the geometry of Qt in ^t- 

For this reference term closed two-surfaces in Minkowski space have QLE zero. 
What is more, for a two-sphere of constant r and t in Schwarzschild space, the QLE 
— i> m as r ^ 00. Further, it is with this reference term that the QLE was first 
shown to be equivalent to the ADM mass 
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There are still problems. Take a spherical set of observers in flat space and 
give them a radial boost. Then as shown in equation ( 3.35|) , e = 'ye and so the 
QLE is (7 — 1)/?. Again it is non-zero in fiat space and actually divergent as 
R ^ (X). That is bad enough, but there is an even more serious concern. As 
Brown and York recognized in their paper, in general it isn't possible to embed a 
two-surface in fiat M^. There are theorems that say (see for example that any 



Riemannian manifold with two-sphere topology and everywhere positive intrinsic 
curvature may be globally embedded in M^. However, most surfaces don't have 
such an intrinsic curvature and once that restriction is broken it is easy to find 
surfaces that cannot be embedded. For example, a surface of constant r and t 
(Boyer-Lindquist coordinates) in Kerr space cannot, in general, be embedded in 
M^. For small enough r (though still outside the horizon and even the ergosphere) 
the intrinsic curvature goes negative sufficiently close to the poles and it is not hard 
to show that the surface cannot be embedded in three-dimensional fiat space. For 
a further discussion of this point see [ff^ . 



Embedding Qt in a 4D reference space 

The Brown- York reference term may be naturally generalized to deal with the 
problem of moving observers |T^. Then, instead of embedding (fij, aap) in a three- 



dimensional reference space, embed it in a four-dimensional reference spacetime A4 
and define a timelike vector field over the embedded surface such that 

1) = 

2) V"' = (in the sense that their mappings into Vtt are equal), and 

3) £T0'ai3 = £to^r (in the sense that their mappings into VLt are equal). 
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Together the first two conditions imply that the boundary lapses N and N_ are 
equal as well as the boundary shifts. The third says that the time rate of change 
of the metric is the same in the two spacctimcs. Physically these conditions mean 
that an observer living in the surface fit and observing only quantities intrinsic to 
that surface (as it evolves through time) cannot tell whether she is living in the 
original spacetime or in the reference spacetime. That is, locally (in the time sense) 
B is embedded in A4. Prom a physical point of view the observers have calibrated 
their instruments so that they will always measure the quasilocal quantities to be 
zero in the reference spacetime, no matter what kind of motion they undergo. 

Then define 

/ = JdtJ d^x^/^lNe-VjJ, (3.39) 

where e and are defined in the same way as before except that this time they 
are evaluated for the surface Qt embedded in the reference spacetime. Thus, the 
net effect of including / is to change s ^ s — e and Ja ^ 3a — 3^- 

With this reference term, the transformation laws for the quasilocal quantities 
change. Consider unboosted observers evolved by it" and watching and 
observers. Then, in general 7] = u°'na will not be equal to ?7 = u^n^. Physi- 
cally this means that in order for [flf, aap) to evolve in the same way in the two 
spacetimes, that surface will have to "move" at different speeds in each. Then the 
transformation law for the quasilocal energy density with reference terms becomes 

£ - f = + r]e^^ - + ?7£l^ . (3.40) 

With this definition of / the problem of observers in flat space seeing non-zero 
energies is solved. Taking Minkowski space as the reference space it is trivial that 
fit may be embedded and defined. Simply leave fit as it is and define = T". 
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Then observers undergoing any motion in Minkowski space measure zero energy. 
Similarly the action is zero for any region of flat space. 

However, problems remain. In the first place even though it is always possible 
to (locally) embed a two-surface in Minkowski space (see for example ref. |T5|), 



that embedding will not be unique ||3^. Thus, the problem of existence has been 
replaced by a question of uniqueness. Furthermore, there is no guarantee that the 
desired vector field can even be defined so even the existence problem has not been 
fully eliminated. 

Nevertheless for the problems that are considered in this thesis this definition of 
the reference terms will suffice. A good discussion of a closely related reference term 
(that combines aspects of this approach with those reviewed in the next section) 



may be found in [34 



Intrinsic reference terms 

Recently there have been several proposals for reference terms / that are defined 
with respect to the intrinsic geometry of B, rather than its extrinsic geometry 
after it has been embedded in M. Most but not all (for example Lau ^6] has 



a different motivation) of these so-called counterterms have been inspired by the 
AdS / CFT correspondence and are intended to remove the divergences of the action 
/ — / without having to worry about the existence or uniqueness of embeddings 
or for that matter what is the proper reference space to use - a non-trivial issue 
if one is considering more complicated spacetimes such as AdS space with periodic 



identifications |70|, H or NUT black holes g|, |3, |8|. 

Typically such terms are defined with respect to the Ricci scalar of B or Qf as 
well as other intrinsic scalars - their exact form depending on the dimension of the 
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spacetime in which they are being defined. The original proposal ^ only 

worked for AdS spacetimes but later work allows for asymptotically flat spacetimes 



as well I^B], |63| . 

I'll briefly consider one such proposal here. Its advantages and problems are 
typical of the wider class of intrinsic counterterms. For asymptotically flat space 
Lau and Mann ^] suggested using 



J = - / £xN^Vw!^ (3.41) 



where R^"^^ is the Ricci scalar for VLt. Lau showed that asymptotically, for a static 
set of observers, this reference term agrees with the embedding reference term of 
Brown and York and so the quasilocal energy is equal to the ADM and Bondi 
energies. 

Unfortunately for a flnite region of flat space the quasilocal energy deflned with 
this reference term will not, in general, be zero. The reason for this is easy to 
see. Recall from elementary differential geometry that the mean curvature of a 
two-surface in flat three-space is Cm = {i^i + /«2)/2 and the Gaussian curvature 
is Cg = K,iK,2 where Hi and k,2 are the principal curvatures of the surface. Now, 
the contracted extrinsic curvature k = 2Cm and the Ricci scalar i?*^^^ = \J2Cg so 



k > V2i?(2) simply because the arithmetic mean of two quantities is always greater 
than or equal to their geometric mean. The equality only holds if the two principal 
curvatures are equal. That is, the two are only exactly equal when Qt is a sphere. 
Lau showed that if a rigid surface Qt is blown- up to its asymptotic limit, then 



k V2R^. However, flnite regions of flat space have non-zero energy with this 
reference term unless their boundary is a sphere. 

Before moving on to the next section it is as well to emphasize once again that 
any of these choices of reference terms are perfectly acceptable from the point of 
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view of the action and/or Hamiltonian generating the correct equations of motion. 
The exact form is only important to set the zero of the quasilocal quantities. 



3.6 Thin shells 



In this section I examine in some detail a correspondence between the quasilocal 
formalism and the mathematics describing thin shells in general relativity which 
was developed by Lanczos and Israel ^9\. This was noted in passing in but 
here it will be examined in more detail and used to reinterpret the quasilocal energy 
from an operational point of view. Following that, I'll briefly apply some results 
from the previous sections to explore the physics of thin shells. 



3.6.1 The thin shell/QLE mathematical equivalence 

Israel considered the conditions that two spacetimes, each with a boundary, must 
satisfy so that they may be joined along those boundaries and yet still satisfy 
Einstein's equations. He showed that as an absolute minimum the spacetimes 
must induce the same metric on the common boundary hypersurface. Further the 
Einstein equations will only be satisfied at the boundary if its extrinsic curvature 
in each of the two spacetimes is the same. If those curvatures are not the same 
then a singularity exists in the (joined) spacetime at the hypersurface. However 
the singularity is sufficiently mild that it may be accounted for by a thin shell of 
matter defined on that boundary. The change in curvature may then be interpreted 
as a manifestation of the thin shell of matter. 

Modifying Israel's notation and sign conventions to be compatible with those 
used here, the stress-energy tensor of that matter is defined as follows. Consider a 
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spacetime M. divided into two regions V"*" and V~ by a timelike hypersurface B. Let 
the metric on V*" be g'^g and the metric on V~ be g~g, and assume that they induce 
the same metric on B. Further, let and n'^ be the spacclikc unit normals 
of B on each of its sides (both oriented to point in the same direction) and define 
©^^ and to be the extrinsic curvature of B in V"^ and V~ respectively. Then, 
Einstein's equation will only be satisfied if a thin shell of matter is present at B 
with stress-energy tensor S^/s = ^ {{^ap ~ ©"^Ta/s) - (0«/3 - ©"Ta/?)}- Note that 
this is written as a tensor field in the surface B. To write it as a four-dimensional 
stress-energy tensor an appropriate Dirac delta function must be included. 

Now let Qt be a foliation of B corresponding to a timelike vector field T° = 
TV-u" + V"" (which as usual lies entirely in the tangent space to B). Then observers 
who are static with respect to the foliation will observe the thin shell to have the 
following energy, momentum, and stress densities: 

S^Sapu'^u'^ = -{k+-k-}, (3.42) 

Ja^-Sys(^Zu^ = -{alu^V^nj -alu^V^nj} , and (3.43) 

^ S^salal = -{Ck+g-{k+ -n+'as)aap) (3.44) 
-iKp - {k- - n-^as)aaf3)} , 

where = —alapVjfif and k^ — cr°'^k^p are the extrinsic curvature of the 
surface Qt in a (local) foliation of M. perpendicular to B. a" retains its earlier 
meaning. 

The correspondence between the quasilocal and thin shell formalisms is now 
obvious. Consider the surface {Bj'jais) embedded in a spacetime {A4,ga/3) and a 
reference spacetime (Ai, q ^). Further let {M.^ga(3) be isomorphic to (V+.f?^^) (or 
more properly the portion of (A^, gap) to one side of B is isomorphic to (V^, g^a))-! 
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and in the same sense let {JSd.,9^f^) be isomorphic to (V ,9^(3)- Then for observers 
hving on B and defining their notion of simultaneity according to the foliation VLt, 

E = e-e, (3.45) 
Ja = Ja-J^, and (3.46) 

Sap = Sap - lap, {3A7) 

where s^,^ is defined in the obvious way and the energy density of the matter seen 
by the observers is 

T'^Sapu" = NE- V"Ja = Nie - e) - V^iJa - jj. (3.48) 

This mathematical identity of the formalisms can be interpreted in couple of 
ways. First, following one can note that the quasilocal work formalism provides 
an alternate derivation of the thin shell junction conditions and stress-energy tensor. 
Namely consider two quasilocal surfaces on either side of the shell and consider the 
limit as the two go to the shell. In that case any reference terms will match and 
cancel leaving only the the stress-energy tensor defined above. This derivation is 
quite different from the one used by Israel. 

From a slightly different perspective the thin shell work can be seen as provid- 
ing an operational definition of the quasilocal energy with the two-surface-into-4D 
reference terms. Given a reference spacetime Ai which is assumed to have energy 
zero, then the quasilocal energy associated with a two surface Qt and time vector 

in a spacetime Ad can be defined as the energy of a shell of matter in Ai 
that has the same intrinsic geometry as Qt (including the rate of change of those 
properties) and a matter stress-energy tensor defined so that the spacetime outside 
ofQ^ is identical to that outside of fit in while inside it remains Ai. In fact, the 
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quasilocal energy with the embedding two-surface-into-4D reference terms consid- 
ered in the previous section is defined if and only if the fields outside of Qt can be 
replicated by a shell of stress-energy with the same intrinsic geometry embedded in 
M. . Provided that Qt and can be embedded in the reference spacetime Ai, the 
construction considered in this section defines the relevant stress-energy for a shell 
in M. 



3.6.2 Physics of thin shells 



Finally, there is a nice application of equation (p.21|) to thin shells. Using that 



equation, including the reference term, and assuming that the reference space is a 
solution to the Einstein equations, 

£THt = j (fx^\££TN-Ja£TV''-^S''''£T(yah^- (3.49) 

The stress tensor can be further decomposed as it was at the end of section p^ . 
For the stress tensor 

^s''^£T<yab = p£tV^+ V^T]^^£TCTab, (3.50) 

where p = cr^'^Sab is a pressure and 77"^ = s"'^ — {l/2)pa°'^ is a shear. The reference 
space stress tensor can be broken up in the same way. 

The terms of the above can be individually interpreted. The ££tN term records 
how the energy measured changes with how the observers choose to measure their 
time (remember that in a Hamiltonian approach energy is conjugate to time so if 
one measures time as going by more quickly then one also measures a larger energy). 
The Ja£Ty°' evaluates the change in the energy contribution from matter fiowing 
around the shell - if the observers change their motion then they will observe 
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different matter motions and so see a different energy. The part of the stress tensor 
corresponding to j9£tv^ — p£T_-\/a term measures energy expenditures required 
to rigidly shrink or expand the shell, while the rj"'^ £TO'ah — v'^^^TO'ab part records 
the work done to deform it. These are all terms that one would expect based on an 
intuition on how classical, non-relativistic membranes under tension should behave. 



Chapter 4 

A quasilocal Hamiltonian for 
matter 



The analysis of the previous chapter can easily be extended to include matter fields 
when those fields have a Lagrangian formulation. Such an extension was made (in 
the orthogonal case) for dilatons and general gauge fields in ref. ^ but for 
purposes of this work I just need the coupled Maxwell and dilaton fields that were 



discussed in section 2.2 



In this chapter I will consider a Lagrangian formulation of the field equations 



from section |4.1| , and then derive an equivalent Hamiltonian in section [4.2| . The 
field equations will then be seen to follow from that Hamiltonian, though it will be 
seen that the formalism itself puts restrictions on the matter field configurations 
that it can be used to study. Issues such as the transformation laws and thin shells 



will be briefiy reconsidered in the light of the new matter terms in section 4.3 



Finally in section [4.4| I'll examine all of this in the light of the duality that was 
considered in section [2.2.4. 



60 
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From a Hamilton- Jacobi perspective parts of this chapter were pubhshed in [|rT 
but in pure Hamiltonian form they appear here for the first time. 



4.1 The gravity- Maxwell- dilat on Lagrangian 



The field equations ( |2.12| , p.l3| , and p.l4| ) are generated by the variation of the 
action 



I = — 
2 



^ / c/^Xv^(7^-2A-2(V„0)(V» -e-2'^*F«^F"^) (4.1) 

+ - / d^xVhK -- f S 
1^ Jt. J b 



where is the dilaton field, Fa/B is the electromagnetic field tensor, and a is the 
coupling constant between the two fields. I assume that at any point in M, Fa/s 
is defined with respect to some gauge potential one-form Aa such that Fap = 
d^Ap — dpAa- The existence of these vector potentials means that e'^^'^^V pF^^s = 
(equivalently d{dA) = 0) so before applying the variational principle the Maxwell 
equation (|2.11|) has been assumed. 

Taking the first variation of the metric terms with respect to the metric, gauge 
potential, and dilaton, it is straightforward to obtain 

S {-2^{Va(l>) (V"0) - v^e-2-*F,^F-^) (4.2) 
= 4:V^^Dii5(i) + 4:^/^J^^^^,j5Ap - K^f^Tap5g'^'^ 

where, J"!^ = V« [g-^'^'^F'^^] , Ton = V°V«0 + (l/2)ae-2'^<^F„^F"^ and T^p was 
defined in equation ( p^.lSj ). The equations J^^j^j = and J-'du = are equivalent to 
equations ( p.l2|) and (|2.13|) respectively. 
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Then, assuming that there exists a single gauge potential Aa covering the entire 
region M (an assumption that I will have more to say about in section |4.2.1|) 
Stokes's theorem can be used to move the total divergence out to the boundary of 
M. Combining this with the vacuum result ( p.2|) the total variation of the action 
is 

SI"" = ^ [ d'xy^\{G^p + Ag^p-8nT^p)6g^'' + 4J^mM + AJ^^EMSAp} 
+ j rf=^x|p"^5/i„^ + p(50 + ^"(5i,} + J rf\{P^5(v^)} (4.3) 

where P°'^ and P ^ retain their earlier meanings and S°', p, and A^ are the den- 
sitized electric field, the densitized time rate of change of the dilaton, and the 



three-dimensional gauge potential as discussed in detail in section |2.2.2 . 

Fixing the metric, vector potential, and dilaton on the boundaries of M, and 
solving 61^ = the Einstein, Maxwell, and dilaton field equations must hold. 
Equivalently this particular action is only fully differentiable if those quantities are 
fixed on the boundary. 



4.2 The gravity- Maxwell- dilaton Hamiltonian 
4.2.1 Form of the Hamiltonian 

From this action it is a fairly straightforward calculation to derive the corresponding 
quasilocal Hamiltonian. As in the previous chapter the action has to be broken up 
with respect to the foliation and then the Hamiltonian and momentum terms picked 
out from the detritus. Details of the calculation can be found in appendix [A.3| but 
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the foliated action is 



+ jdty^ d'x {P^£tV^) - i/r} > 



(4.4) 



where 



[ d^x[NT-r + V^T-C^ + T''AM] (4.5) 
+ / d'xV^ [N{e + e^) - V"(J, + J^)] . 

T°'Aa = — + V^Aa in terms of quantities defined on the hypersurfaces, while 
Q = -DpS^ = is the free space version of Gauss's law from electrodynamics 
(equation ( p.26| )). 

^ = --^(n^^'5)(i$-r7i,n") = --^(r2^£'3)<land (4.6) 
= -^{np£P)A^ = -^{np£P)A^ (4.7) 

which can be identified with energy and angular momentum as suggested by the 
notation. The bar retains its usual meaning, so in this case $ = —Aav!^ and 
E"' = -2v^/kF"^u^. Note that i„ = a^Af^ and n/^S^^ = np£l^ and so are left 
invariant by the bar notation. 

Then, the electric field vector density E"" and the dilaton rate of change p 
are identified as momenta conjugate to Aa and (f) respectively. Exactly what is 
happening with the T°'Aa term isn't clear at this stage, but after calculating the 
variation of H]^ in the next section it will be clear that $ (the Coulomb potential) is 
a Lagrange multiplier. Finally, if™ can be identified as the Hamiltonian functional. 
As in the previous chapter, the numerical value of H^^ evaluated for a particular 
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leaf of the spacetime foliation depends only on the boundary f2t and how that 
boundary is evolving in time. 

In the next section I will show that the functional HJ^ really does generate the 
correct field equations, but before moving on there are a couple of points to consider 
regarding the electromagnetic gauge potential Aa and gauge invariance. 



No magnetic charges allowed 

In the derivation of the Hamiltonian from the action it was assumed that there is 
a single vector potential A^ defined over all of M. This assumption meant that 
total derivatives in the bulk could be removed to the boundaries under the auspices 
of Stokes's theorem. However, a corollary of this assumption is that there is no 
magnetic charge in M (or contained by any surface that is itself contained in M). 
The next few paragraphs explore this statement from three closely related points 
of view. 

As a start, let Qx be any closed spatial two-surface in M with normals 
and n". Then, the magnetic charge contained within fix is Jq^ (Px^/an"Ba■ By 
equation ( p.44| ), n'^Ba = —n'^ul^e^p^^D^As = e^^'^rfaA/j) where da is the covariant 



derivative in the surface Qx, ^cf is the Levi-Cevita tensor on that surface, and 



again Af^ = crjA^. But this is an exact differential form and so integrated over a 
closed surface it is zeroQ. Thus there is no magnetic charge contained by any surface 
in M. 

Keep in mind that this is a stronger statement than just the local statement that 



^In the more efficient differential forms notation, in the spatial slice orthogonal to m", A is a 
one form and B = dA is a two form. Then if B and A are the forms projected (or pulled-back) 
into VLx, the magnetic charge contained within fijc is B = J^^ dA = since ilx is closed. 
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F = dA =^ dF = d{dA) = ^ DaB°' = 0. When working with a gauge potential, 
the manifestation of magnetic charge in the potential is global and topological 
(resulting from a twist in the U{1) gauge bundle) rather than local as is the case 
for electric charge. If one assumes that there is a single Aa that covers M then the 
f/(l) gauge bundle is trivial by definition and so there is no magnet ic/topological 
charge. Even more strongly, as just noted, no surface contained in M can itself 
contain magnetic charge. This means, for example, that if M is the region bounded 
by two concentric spheres (multiplied by a time interval), then not only is there no 
charge in M but also there is no charge in the region inside the inner sphere. 

In fact, projecting into spatial slices Sf of M, de Rhams theorem (see for example 



[36 1) says that a single vector potential is defined over all of if and only if there 
is no magnetic charge contained within any two-surface Qx C S^. Thus to allow 
for a magnetic charge in M, one must break the region of spacetime into at least 
two regions each of which has its own vector potential. Then, the frequent uses of 
Stokes's theorem in the derivation will remove total divergences to the boundaries 
of those regions rather than just the boundary of M itself. By definition some of 
those region boundaries will actually be interior to M and so observers inhabiting 
dM will not be in a position to measure all of the boundary terms and therefore 
will not be able to fully assess what is happening in the interior of M. 



The gauge dependence of the Hamiltonian 

Note that even though the action is gauge invariant (ie. it depends only on Fa/B 
and not on the exact form of A^) the proposed Hamiltonian doesn't necessarily 
inherit that invariance. The paths by which this gauge dependence can creep in are 
quite easily found but at the same time the effect is important so I'll pause here to 
point them out in some detail. 
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First, by equation ( [4 .41 ) it is clear that while the time integrated difference 
between the Hamiltonian and "kinetic energy" terms must be gauge invariant, 
that invariance can only be inherited by the Hamiltonian itself if part of the gauge 
freedom is used to ensure that = 0. If this is the case and H^^ is independent 

of the leaf of the foliation, then if™ will be independent of the remaining gauge 
freedom. For stationary spacetimes that gauge and foliation are, of course, the 
natural ones to choose but it should be kept in mind that a partial gauge fixing is 
required to ensure that the Hamiltonian is invariant with respect to the remaining 
gauge freedom. 

In the conventional usage of this work to study black holes, there is an alternate 
route by which gauge dependence can find its way into the Hamiltonian. Namely, 
components of the gauge potential Aa may diverge on the (apparent) event horizon. 
Then, A^ has a singularity in M and so the uses of Stokes's theorem in the derivation 
of the Hamiltonian aren't valid. To avoid this problem one could cut out the 
horizon with a inner timelike boundary B' , though in that case the region between 
B' and B would be under consideration rather than the full region contained by 
B. This problem is usually ignored however and to facilitate comparisons between 
singular and non-singular spacetimes, only the outer boundary is considered. In 
section |5.1.1j , this version of gauge dependence will be demonstrated for a Reissner- 
Nordstrom spacetime, and it will also be seen that this gauge dependence arising 
from neglecting the inner boundary amounts to little more than a choice of where 
to set the zero of the electromagnetic energy. 

Thus, it can be seen that gauge independence of the action doesn't necessarily 
ensure the gauge invariance of the Hamiltonian. Specializing to the case where the 
lapse = 1 and shift V"' = on the boundary, and leaving aside the issue of how 
such a choice affects the relative foliations of inner and outer boundaries, note that 
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the QLE will not in general be gauge independent either. 
4.2.2 Variation of the Hamiltonian 

This section will show that the proposed Hamiltonian really does generate the 
correct equations of motion. To do this consider Ht as a functional of the surface 
St, its boundary Qt, the normal ra" to that boundary, the fields hab, \f^i 0, A^ along 
with their conjugate momenta P"'', -P^, p, i^", and the Lagrange multipliers A^, 
and $. In the usual Hamiltonian way, the conjugate momenta are considered 
to be entirely independent variables. Their connection with hah-, a/o", 0, and Aa is 
forgotten. Following the lead of section p.2.2| , e, J^, iV, V^", A, and r] are defined 
entirely with respect to V^", A^, and n^. Similarly $, e, and Ja can be written with 
respect to $, A^, A, r], and n^. 

Then, the variation of H'^ with respect to the quantities hah-, \f^-, 0, and Aa-, 
their conjugate momenta P°'^ ^ P^fs-, Pi and and the Lagrange multipliers A^, V^", 
and $ is 



5H\ 



m 




(4.8) 
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Details of this calculation can be found in appendix |A.4| but for now note that 



[hablx retains its meaning from the previous chapter (equation ( p^) ) while 

„A Nk 



T 



K 



^[DMD^<j>]h^' j + (4.9) 



IT 



A., 

T 

B 



Nk 



D'[ND,^] + a^[e-''''''B'Bb - e^'"^£'£k] + D,[pV'] 
e'''^£a + £vA,-D,[N^], 



2Vh 

-e''^''Db[Ne^''*Bc\ + £v£\ 

\& -'qe^'"t't'"^£'^, and 
A 



(4.10) 

(4.11) 

(4.12) 
(4.13) 
(4.14) 

(4.15) 



£"■ = o'^£^ and B"' = a^B^ are the projections of the electric and magnetic vector 
densities into the tangent bundle of the boundary Qf. 

The Hamiltonian equations of motion can now be obtained by calculating the 
full variation of the action ( |3.3D (treating the momenta as independent variables) 
and solving 61 = 0. Then an application of the fundamental theorem of calculus 
to remove the time derivatives to the spatial boundaries and a little bit of algebra 
shows that 



61-61 



j d'x (^P^'6hab + £^6Aa + p6<p) + j SxP^6^ (4.16) 
- j dtj d^x |(H™ - <^Q)6N + (7^™ + AaQ)6V'' - ArQ5$} 

+ j dtj^ d^x{{£TKb-[hab]T)SP'''-{£TP'"'-[P''']j,)SKb} 

6£'' - {£t£'' ~[£'']T)S^a 



+ I dt I d-'xU £TAa 



Aa 
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+ [dt [ d^x{{£TV^- [V^^)6P^-{£tP^-[P^]^)SV^} 

J JQ.t 

Then, if hah and Aa = h^Ap are fixed on the boundaries Si and S2, and 7^/3 
(equivalently (Tab, iV, V°') and 7^ A/3 (equivalently $ and Aa) are held constant on 
-B, the solution of 5/ = gives the correct field equations. Namely 



rr 

njm 

£tP''^ 
Q 

£TAa 

£tp 

£t<\> 



0, 
0, 

{P^Tt^ 
0, 



A 
[p]t ' 



(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 
(4.22) 
(4.23) 
(4.24) 



as well as equations (|3.16|) , (|3.17| ), and (|3.18|) from the previous chapter. 

Now (|4T^ ), (|41^ ), and ( ^A9|) are the projected Einstein equations (|2]3|), (|2]4g), 
and ( |2.41| ) respectively, so the Hamiltonian has recovered those correctly. Equations 
(|4:20| ) and (|]2l]) are the projected Maxwell equations (|]26D and ( 121271) , while ( ^1231) 
is dilaton equation ( p.30| ). Equations ( [4.22 ) and ( [4.24| ) are just definitions of the 
respective Lie derivatives while ( |3.16D , ( |3.17D , and ( |3.18| ) continue to express their 
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earlier meanings. Keep in mind that the existence of Aa imphed the remaining two 
Maxwell equations (12^2^) and ( ^^ . 



From equations ( p.43| ) and ( p.44|) it is easy to see that fixing $ and Aa on the 
timelike boundary B is equivalent to fixing the component of B°' perpendicular to 
B (that is B^Ua) and the components of E"' parallel to Qt (that is a'^E"'). Thus, 
the action is fully differentiable only if the parameter space of spacetime studied 
is restricted to those with a specified magnetic charge. This fits in well with the 
discussion at the end of the previous section that said that the magnetic charge is 
fixed (to be zero) by the existence of the single vector potential generating the EM 
fields. In contrast there is no restriction on the electric charge. This issue will be 



considered in more detail in section 4.4 



Note too that while the value of the dilaton field (f) is fixed on Qt, its "radial" 
rate of change £n(j) is left free. Therefore, the dilaton charge 

= [ (fx^£n(l), (4.25) 

is not fixed either. 



n 



Comparison with the Lagrangian approach 

Again it is reassuring to compare this Hamiltonian analysis with a Lagrangian 
analysis and in particular show that the variation of the Hamiltonian properly fits 
into that of the action. The matter fields considered above were examined from 
that viewpoint in full nonorthogonal form in [|lT| and for rj = in E7\, p8|. 



Breaking up the matter term of equation (|4.3| ) and bringing in the full variation 
of the gravitational action ( p.20 ) from the last chapter, it is straightforward to show 
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that 



61"^- 61 = ^[ d'x^\iG^p + Ago.p-8nT^p)6g''^ + 4J^mM + 4J^^E,,6AA 

+ J d^'x |P"^(5/i„^ + p5(t) + + j d^x {P^5{^)] (4.26) 

- j dt d^x^\^{^e + e^)5N-{j^+j'^)5V^-Ys'"^5a^p^ 

+ - j dt j d^x^/^N ^^-£fi(p5^+ {npE'^)5^ - e-^'''^Uanpe"'^^^B.^5A 

With this approach the momenta are functions of the metrics, normals, gauge 
potentials, and dilaton field so equations (^), ( |^ ), ( |CT| ), (^^), and (|]^ 



automatically hold. Gap + Ag^p — SttTq,/? = 0, J^du = 0, and J^^m ~ ^ ^^^^^ 
of the field equations and so again the Lagrangian and Hamiltonian treatments are 
equivalent. 

4.3 Properties of the Hamiltonian 

In this section I'll discuss some of the issues that arose in the previous chapter in the 
light of the new matter terms. As will be seen the required changes are incremental 
rather than qualitative. I will also examine the action and Hamiltonian in the 
light of the duality discussed in |2.2.4| and show how it may be used to extend the 



preceding analysis to magnetically charged spacetimes. 
4.3.1 £tHI^, conserved charges, and energy 



By the same arguments as used in section 3.3 the time rate of change of the Hamil 



tonian functional with matter fields included is 



£THr = I d^x^{[e + t'']£TN -[ja+j':]£TV'' -{N/2)r'£Tcyab) 
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+ / d^x^^^ ([S''na]£T^+ e-^^'^Sanbe'''"'£TAc) (4.27) 



,2 2iVy^ 



This is zero if the vector field T" defines a symmetry of all fields on the boundary B 
and so in that case is a conserved charge and is conventionally identified with 
the mass contained by the surface Qf Of course by Noether's theorem it is to be 
expected that a symmetry corresponds to a conserved charge, but once again note 
that it is only symmetries of the fields at the boundary that matter. The properties 
of the fields in the bulk are completely irrelevant. 

Even if HJ^ is not a conserved charge, I'll still label it to be the mass contained 
by Qf. Then, the discussion of section R]^ largely applies here as well. In particular 
one can consider the special case where N = 1 and V"' = and define the quasilocal 
energy 

Etot= [ d^xy^{e + e'^). (4.28) 
Jn 

Note that a gauge choice can be made to set e*" = on Qf in which case this reduces 
to the Eceo defined in equation p. 23 . 



4.3.2 Transformation laws 

It is easy to extend the transformation laws to the matter terms. Again considering 
reference frames associated with the normals (u",n") and (u°,n") and reusing the 



I notation of section |3.4| define 

2 



^ml ^ -(n^Efs)^^, and (4.29) 
^ l(n^Ep)A^, (4.30) 
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where = A^n". Note that j^^ = j™. Then, it is almost trivial to show that 

= ^e^ + r^e™! =7(e™ + i;K£"^) and (4.31) 
A 

= r- (4.32) 



Thus, 



e + e"" = \{e + e'^)+ri{e^ + e'^^) (4.33) 
A 



J {{e + e"^) + v^{e^ + e^^)) , and 



Ja+JT = Ja+j::-^a%v (4.34) 

K 



4.3.3 Reference terms 

Virtually no change is required in the discussion of the reference terms from section 
p3| . In principle, with no implications for the field equations, / could be allowed 
to depend on the matter field terms that are fixed on the boundary. That is J 
is a functional of h^Ap on Ei and S2, 7^^^ on 5, and over all three of those 
boundaries. In practice however, the usual use of the reference term is to calculate 
how different the action of M is from a similar M in an "empty" reference space 
and so this option is not generally taken up. 



4.3.4 Thin shells 

The analogy between the quasilocal formalism and thin shells can be extended to 
encompass the Maxwell and dilaton fields if one allows the shells to have electric and 
dilaton charges and currents embedded in them. These charge and current densities 
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are defined to account for discontinuities in the electromagnet ic/dilaton field just 
as the stress tensor is defined to account for discontinuities in the gravitational 
field/geometry of spacetime. 

If one assumes that there are no electric charges/currents inside the thin shell, 
then calculating the electric charge densities is an exercise from undergraduate 
electromagnetism |^ . Specifically for the foliation Qt of B the electric charge 
density on the shell is —y^/^/hn^Sa- At the same time, given an electromagnetic 
potential A^, observers on the surface of the shell whose evolution is guided by the 
vector field T" will define a Coulomb potential —T°'Aa = — V°'Aa. In the usual 
way the energy of the charge density in the electromagnetic field is then the charge 
times the potential. That is, -^/a/ y/hn°'£ai-T°'A^) = Ne^ - 9°"]^. As usual this 
component of the energy is gauge dependent. 

Similar reasoning gives the dilaton charge on the shell. The dilaton charge in 
a given volume is given by the integral of n"Va0 over the surface enclosing that 
volume. For black hole solutions, the value of the dilaton charge is constrained 
by demanding the spacetime has no singularities on or outside of the outermost 



horizon |^T|. In the thin shell case, n^V^^ then yields the dilaton charge density 
on the shell Qf 

The surface charges do not change the definition of the surface stress energy 
tensor which was defined entirely by the Einstein equations. As such they also don't 
change the definitions of S, and Sap- Therefore, including the stress energy with 
the energy of the shell in the gauge field, the total energy density in a thin shell 
evolving by the vector field T" is iV(e + e^') — ^"(Jq + J^) minus the corresponding 
reference terms. This of course is exactly the same as the Hamiltonian quasilocal 
energy of the region of space on and inside of the shell as measured by a set of 
observers being evolved by the same vector field, and so the correspondence between 
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thin shells and quasilocal energies remains. 



4.4 Electromagnetic duality 



In section [4.2| it was demonstrated that the formalism developed so far only prop- 
erly applies to spacetimes that do not contain magnetic charge. Spacetimes with 



magnetic charge are often of interest however and in particular in section |5.2| , I'll 
want to use the formalism to investigate to naked black holes which are magnet- 
ically charged. Thus, it is of interest to extend the formalism to allow for such 
spacetimes. 

The obvious way to do this is to make use of the duality reviewed in section 
2.2.4| . Applying this duality, the action becomes 



I"^' = ^ [ d^a;v^(7^-2A-2(V„0)(V» -e2"'^^F„;3*i^"^) (4.35) 
2k Jm 

H — / d^x\fhK / ci^xv^0 + - / d?x^/o smhr^{'q) + L 

1^ Jy. J b ^ Jn 

Note that FapF""^ = — -kFap -kF""^ so this action is not numerically equal to J™. 

Breaking up this action with respect to the foliation, one must assume that 
there is a single vector potential A* generating -kFap so that total divergences can 
be removed to the boundary. Then a corollary to this assumption is that d-kF = Q 
or equivalently V/3(e~^"''^F"^) = 0. From section |2.2.2| this relation can be rewritten 



in terms of fields in Et as equations ( |2.26| ) and (|2.27|) which are 



= and 

Recall that assuming a single Aa implied the other pair of Maxwell equations. The 
arguments of section |4.2.1| can then trivially be extended to show that there are no 
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electric charges in a spacetime where the Maxwell field can be described by such a 
single A*. However, magnetic charge is allowed. 

That said, the action may be broken up with respect to the foliation to become, 
/'"*-/ = jdt^j rf^x (p"^£T/ia/3 + pi^T0 + i3"i:ri*)| (4.36) 

and 



I d^xlNH"" + V''H^ + T''A*^Q*] (4.37) 

+ [ d'xV^{N{e + en-V''{j.+f:i}. 
Jfit 

T'^Al = -N^* + VAl and defined in section |]2]|) and Q* = -DpB^ is 
the free space magnetic version of Gauss's law from electrodynamics (derived from 
its 4D form at equation ( |2.24| )). Further 



= — ^(n^i3^)(^$*-r/i>") = — ^(n^i3^)$* and (4.38) 
= -^{n^B^)Al = -^{n^B^)Al (4.39) 

are the new matter energy and angular momentum terms. Note that they are 
different from and which depended on the regular vector potential A^ and 
electric field density E^- The bar retains its usual meaning. 

Sticking to the Hamiltonian perspective that momenta are independent of their 
corresponding configuration quantities, the total variation of this action is 

5/™* - 5/ = y d^x (p^^bKb + B'^bAl + + y d^xP^b^ (4. 

- j dt j c^^x|(7^'"-$*Q*)(5A^+ (7^^ + i^Q'')5V"- A^QW} 
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+ I dt I d^x\ {£tA 



A' 



+ j dtj^ d''x{{£T^-[^]T)6p* -{£tP* -[p*]t)S^} 
+ j dt {{£t^/^-[V^^)5P^-{£tP^-[P^]t)5^] 
- j dt j^d^x^l^ {e - £™*)5iV - {ja - jTW'' - ^r^aab 

where [P"'']™, [habjr, [p]t, [4>]t and [0]^ retain their earher meanings while 

Nk 



Al\^ ^ —^e-'^^Ba + £vA:-D^[N^^ 

[B\ = e''^'Db[Ne-^''^£,] + £vB'', and 

B = + r/e-2'^V"4. 

A 



(4.41) 
(4.42) 
(4.43) 



Applying the duality relations it is easy to see that these equations of motion are 
equivalent to the earlier ones. Note however, that the electromagnetic quantities 
that must be kept constant on the boundaries have changed. Specifically on Si and 
S2, A* must be kept constant while on the B boundary 7^ (or equivalently $* 
and cr^Al) must be held constant. This corresponds to holding £an°' and B"" = a'^B^ 
constant which means that now the electric charge must be held constant (at zero 
by the earlier comments) while the magnetic charge is not fixed. 

Thus, there are now well defined formalisms which can be used to study space- 
times containing either electric or magnetic charges. What is missing is a formalism 
that easily handles dyonic spacetimesQ. A duality rotation could be used to study 



^That is those with both electric and magnetic charges. 
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a spacetime with a particular dyonic charge but even then there would still be 
problems in spacetimes containing multiple dyons with varying ratios of electric 
and magnetic charges. Furthermore, there is something fundamentally unsatisfying 
about having the form of the action depend on the charges contained in the space- 
time. As it stands, I don't have a solution for this problem and so will not consider 
dyonic spacetimes in this thesis. 



Chapter 5 



Classical applications 



In this chapter I will apply the quasilocal energy derived from the Hamiltonians 
of chapters three and four to investigate a variety of spacetimes. First, in section 
|5.1| , I'll provide some orientation for the reader by examining the quasilocal ener- 
gies seen by static and moving observers in Schwarzschild and Reissner-Nordstrom 
spacetimes. Section |5]^ will then apply the work to study naked black hole space- 
times. These spacetimes are characterized by the fact that static and infalling ob- 
servers experience very different tidal forces from each other, and that section will 
demonstrate that they also measure very different quasilocal energies and explain 
the connection between the two results. Finally, section ^]3| applies the formalism 
to calculate energy transfers during gravitational tidal heating such as that seen in 
the Jupiter-Io system. This last section can then be seen from two different points 
of view. First of all it can be seen as an alternate way to calculate the magnitude 
of these effects from the usual Newtonian or pseudo-tensor methods or secondly it 
can be viewed as a test of the formalism to see if it produces the standard answers. 



Note that section is based on equivalent sections in and |TT| . The work 
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found in section |5.2| was published in |TT| , while section |5.3| formed part of 



5.1 Reissner-Nordstrom spacetimes 

This section examines the quasilocal quantities measured by observers undergoing 
various motions in Reissner-Nordstrom (RN) spacetimes (and Schwarzschild as a 
special case). In standard form the RN metric is 

ds^ = -F{r)de + + r\de^ + sin^^V), (5.1) 

where F{r) = 1 - iHi + fia^, m is the mass, and Eq and Go are respectively the 
electric and magnetic charges of the hole. The accompanying electromagnetic field 
is described by the two form 

E 

F = --^dtAdr + Gosinede Adip, (5.2) 

while a local vector potential generating this field is 

A = dt — Go cos Odip + dx, (5.3) 

r 

where x = x(^) ^) V') is any function defined over M. For x = 0, note that A is 
not defined for all of M since dp is not defined on 6 = 0, vr. This is in accord with 
the discussion of section |4.4| where I showed that the Lagrangian and Hamiltonian 
formalisms as constituted are not suitable for discussing dyonic spacetimes. In 
the {Faf3, Aa) form considered here, a single Aa cannot describe the field due to a 
magnetic charge and so this property of Aa is not just an annoyance that can be 
removed with a clever gauge transformation. Given this difficulty I set Gq = and 
focus on electric black holes in the following subsections. The results for magnetic 
black holes are identical if one switches Eq and Go and adds in the appropriate -k^s 
and minus signs. 
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I begin the study of these spacetimes by calculating the quasilocal quantities 
measured by a static, spherically symmetric set of observers. Their observations 
will be the subject of the next subsection, while the two that follow will compare 
their measurements with those of a boosted set who instantaneously coincide with 
them on a surface fit, which is a surface of constant r and t. The observers are 
evolved by the vector field T" = N{r)u'^, where N{r) is the lapse function while the 
shift V"" = 0. For any choice of N{r) the observers will be static in the restricted 
sense that they don't observe any changes in the Qt metric cTq/j, and the lapse just 
determines how they choose to measure their time on the surface B. In particular, 
choosing N(r) = a/ F{r) they measure time according to the coordinate t and T" 
corresponds to the timelike Killing vector for the full spacetime metric Qais, while 
choosing N{r) = 1 the observers measure proper time (that is T^T^ = —1). 

Then, with u° = --j==df and = a/ F{r)d", where and are the coor- 
dinate forms of the vector fields ^ and ^ respectively, a series of straightforward 
calculations yields 

1 



\W, (5.4) 



47rr2 ' ^ ' 



1 



(5.6) 
(5.7) 



As usual I am working in geometric units where G = c = /?-=l, soK = 87r. The 

reference terms are defined by embedding the sphere f2t statically in the obvious way 

in Minkowski spaceQ. Since this is a spherical set of observers in a static spacetime 

the ja angular momentum terms vanish. 

^ Clearly this embedding isn't unique. Still, it is convenient and given the many issues involved 
in choosing reference terms (section ^!^) this is all that I will ask for! 
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Strictly speaking, to rigorously apply the quasilocal formalism to black holes I 
should include an inner boundary B' as well as the outer boundary B, in which case 
M would be homeomorphic to x 5*^ rather than M^. Without such a boundary 
to remove the collapsing matter /singularity at the centre of the black hole from 
consideration, the quasilocal formalism is not properly constituted since these were 
not accounted for in its setup. Even for an eternal black hole foliated with Einstein- 
Rosen-bridge hypersurfaces that avoid the singularity, there is a difficulty in that 
the leaves of the foliation all intersect on the horizon. For a good discussion of 



this problem see ref. |37], but in the following I will generally consider only the 
outer boundary so as to facilitate comparisons with non-singular spacetimes (such 
as stars). When studying the quasilocal energy of black holes it is conventional to 
proceed in this way. In some sense it is equivalent to using Gauss's law to calculate 
the electric charge of a point particle without worrying about the divergence of the 
fields at the particle itself. 



5.1.1 Static observers 
Static geometric energy 

First I calculate the part of the quasilocal energy associated with the density e. 
Following section |3.3|, I label it the geometric energy since it depends only on the 



extrinsic curvatures. It can be thought of as the full QLE with a gauge choice made 
so that = 0. Then, 



Geo 



/ d^xV^e = -Jr^ - 2mr + El (5.^ 
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In the large r limit this becomes Eceo ~ — r + m + ^(m^ — Eq). The e_ reference 



2r* 

term is 



so 

1 



(5.9) 



Eceo -E = r- ^Jr^ - 2mr + E^ ^ m + ^(m^ - E^), (5.10) 

in the large r limit. Note that Eceo — E_ monotonically decreases as r increases, 
starting at 2m on the horizon and reaching a minimum of m at infinity. Thus by 
this measure the energy contained in the fields is negative which is to be expected 
for a binding energy such as gravity. 

For large r these are the results that would be expected from an application of 
Newtonian intuition to the (equivalent) thin shell situation. From this viewpoint 
consider how much energy it would take to construct a shell of radius r with mass m 
and charge Eq. First it would cost m units of energy to create the required mass at 
spatial infinity where there would be no interactions and so no deviation from the 
rest energy. Then, using Newton's and Coulomb's laws, and assuming that mass 
and charge are equally distributed throughout the matter, it is straightforward to 
show that —^{m? — Eq) units of work are required to assemble the shell out of the 
created material out at infinity. It is then very natural to say that this energy is 
"stored" in the field, outside radius r. Thus, assuming that conservation of energy 
holds once the matter is created, the energy contained on and/or inside the shell 
with radius r is 

(total energy) — (energy in fields outside the shell) = m + ^("^^ ~ -^o)' (^-H) 
as was calculated above. This limiting case was first considered in the original 



Brown and York paper |22 
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Note that for an extreme black hole where \Eo\ = m, EGeo — E_ = m is a constant. 
From the Newtonian shell point of view this makes sense. During the construction of 
the shell out of particles which also have equal mass and charge, equal but opposite 
electric and gravitational forces would act on the particles. Thus, no work must be 
done to build the shell and so no energy is stored in the fields. Alternatively equal 
amounts of positive and negative energy are stored in the electric and gravitational 
fields and cancel each other out. The only energy is then that stored in the mass 
and so the energy contained by Vtt is m for all radii greater than m. 

Note that even though this is the "geometric energy" with the matter terms 
omitted, it certainly seems to include the energy contributions from the electric 
as well as gravitational fields. Thus one may think of this geometric energy as a 
"configuration energy" that arises from the spatial relationships of different parts 
of the spacetime to each other. By contrast in the next subsection where the gauge 
dependent terms are included, the energy also includes "position" terms that arise 
due to the position of the different parts of the spacetime in the gauge potential (a 
point of view also explored in thin shell section |4.3.4| ). Of course, the form of the 
gauge potential is determined up to a gauge transformation by the matter so this 
view of the terms as being configurational versus positional is at best a rough way 
to think of them. 

Static total energy 

Next consider e + e™, the full energy density that was derived from the variational 
calculations (as opposed to the geometric energy which is gauge fixed so that the 
$ = on fi). Then 




(5.12) 
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As has been emphasized before this expression is manifestly gauge dependent. Even 
worse however is the fact that this energy will in general diverge at the outer horizon 
of a black hole. Before I deal with that worry however, consider the usual r — > cxo 
limit. 

Demanding that has the same spherical and time translation symmetries as 
the spacetime, x = ~^oot + fir) where $00 = lim^^oo $ is a constant and /(r) is 
an arbitrary function of r. Then, 



Since the total energy is the sum of the geometric energy and the gauge dependent 
term, it isn't surprising that this Newtonian limit is the sum of the Newtonian limit 
of the geometric energy and the "positional" potential term. One can think of $00 
as the zero level of the potential throughout space (it remains even if Eq 0) and 
so by the thin shell analogy -Eq^oo is the energy cost to create matter with charge 
Eq at infinity (apart from the energy costs associated with the mass). For extreme 
black holes recall that Eqco — = m so Etot — E_ — nT' + Jq^ ^x^e^ and the only 
energy is the mass m plus the energy of the charge with respect to the potential. 

In most situations the exact choice of gauge is just a matter of convenience. 
For black hole spacetimes however, Ei^t will diverge on the horizon with most 
gauge choices. This divergence can be directly traced to the fact that the Coulomb 
potential $ = —u°'Aa — ^{^ — dtx) also diverges at the horizon. To remove both 
divergences choose x such that dtX — ^ as r — ^ r_|_, where r+ is the outer black 
hole horizon. That is, set the Coulomb potential to zero on the black hole horizon. 
Then, assuming that A^ has the symmetries discussed above $00 — Making 
that choice, a little algebra leads to 



Etot-E^r 



— 2mr — E^rdtX 



(m+Eo^oo) + ^ im''+E^+2mEo^oo) ■ (5.13) 



i/r^ — 2mr + E^ 




(5.14) 
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where r± = m ± a/ — Eq are the radial positions of inner and outer horizons. 
This gauge will also be used for the naked black holes. For extreme black holes 
= = \Eo\ = m and so Etot = —r — E_ = everywhere. Physically the gauge 
has been chosen so that the electromagnetic potential energy is a constant and 
everywhere equal —m. The (in this case negative) electric potential energy cancels 
the mass-energy while at the same time the positive energy of the electric field 
cancels the negative binding energy of the gravitational field. 

So, as suggested at the end of the last section, the total energy may be split into 
two parts. The geometric part depends only on the configuration of the spacetime 
and examining the Newtonian limit one can see that it appears to include not 
only the gravitational but also the electromagnetic "configurational" energies. By 
contrast this section showed that the gauge dependent part exclusively deals with 
the potential of the matter relative to the gauge field. As has been seen, for a 
given solution to the Einstein-Maxwell equations the total QLE for a given surface 
may take any value (including zero) depending on the exact gauge choice. As such 
it is clear that this gauge dependent part of the energy should not be reflected 
in the geometry of the spacetime as indeed it isn't, since the stress-energy tensor 
Ta/3 doesn't depend on the gauge potential. On the other hand it should not be 
concluded that this gauge dependent part is meaningless. It certainly plays a role 
equal to the geometric energy in both thermodynamics I^T], ^ and black hole 
pair creation (chapter P). 

Value of the Hamiltonian 

Next consider the value of the Hamiltonian as calculated by the same static sets of 
observers who now measure time according to the (Killing) time coordinate t (that 
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Hceo = 


NEceo = 


= -rF, 


Hceo — K = 


NiEceo 


-m = 


Htot = 


NEtot = 


-r + 2 


Htot — K = 


N{Et,t - 


-E) = 


r limit, Hceo - 


H ^ m — 


2r 



is lapse = V^). Then 

(5.15) 

^(1 - Vf), (5.16) 
^ Eodtx, and (5.17) 
+ EodtX + V^-r (5.18) 

[ Htot -H^{m + Eo<l>oo) - 
Thus as is usual for asymptotically flat spacetimes the Hamiltonian corresponds to 
the QLE in the r ^ oo limit. Note too though that even in the large r limit, away 
from infinity these Hamiltonians don't agree with the Newtonian limits discussed 
earlier. In particular the contribution from the gravitational field doesn't have the 
right sign in either case and the electric contribution is also wrong for the total 
Hamiltonian. 

Finally consider the earlier comments on the gauge dependence of the Hamil- 
tonian (clearly seen above by the x dependence) from section |4.2.1| . To avoid the 



complications of singularities in the gauge potential redefine M as the region of A4 
contained by the two timelike hypersurfaces r = ri and r = r2 where r+ < ri < r2. 
Again foliate that region according to the standard time coordinate t. Since I 
am considering the gauge dependence of the Hamiltonian the reference terms are 
ignored since they are gauge invariant. 

Then the total H]^ for a spacelike slice Ht with boundary Qt is 

nri^t] = ^Htot = in - r2) + Yl9tx]:l (5.19) 

where ~ 9tx\r=r2 ~ dtx\r=ri and the sum is over the two boundary compo- 

nents. Thus, at first glance the Hamiltonian appears to be gauge dependent. Recall 
however that section 4.2.1 showed that it could only be expected to be (partially) 



CHAPTER 5. CLASSICAL APPLICATIONS 



88 



gauge independent if M was a region containing no singularities and £TAa = 0. 
Well, there are no singularities in M and a quick calculation shows that £TAa = 
implies that dtX is constant over S^. If this is true then [i9ix]ri ~ ^^"^ Hamil- 
tonian is (partially) gauge independent as expected. 

5.1.2 Radially boosted observers 

This section considers the energies measured by spherically symmetric sets of ob- 
servers who are moving radially towards or away from the gravitational source in 
the RN spacetime. As before, Qt is a surface of constant r and t but this time set 
the time vector T" = iV-u" where = ju°' + rjn°' = 7(u" + v\-n°'). As in section 
v\- is the speed of the observers in the = y/Fd^ direction as measured 
by the static set of observers that I have been working with up until now. 

Then, a straightforward calculation shows = so 



Unfortunately, from the point of view of simplicity, within the gauge freedom e*"! 
is not necessarily zero. Even if one considers only gauge choices that give Aa the 
same symmetries as the spacetime x = ~^oot + f{r), where as noted earlier <l>oo is 
a constant and / is any function of r. Then 



In the interests of simplicity however, I make the standard gauge choice for elec- 
trostatics and let drf = 0. Then the Lorentz-type transformation laws apply and 




(5.20) 



47r 



VFdrf. 



(5.21) 




(5.22) 
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As before, I choose $00 = so that this quantity doesn't diverge at the horizon. 

To include the reference terms, it is necessary to find a time vector for the 
reference spacetime such that T^T^, = T"Tq and £f^p = ^f^ap (the conditions 
from page ^l]). Such a vector field is given by 



(5.23) 



where m° = , n° = and t and r are the usual time and radial coordinates for 
Minkowski space. 

Then 

^ -1^ = ""^^^^^ and 7 = 7a/i - (1 - F)vl (5.24) 



fX yi - (1 - F)t;;^ 
which implies 



and thence 



^ = 7^=-r7Wl-(l-F)t;2, (5.25) 



Eoeo -E = r-f (^^1-(1-F)v^ - Vf^ , (5.26) 



and 







r\fF 



r 



(5.27) 



As they stand these expressions are quite complicated and their physical inter- 
pretation isn't at all obvious. To simplify things a little consider the large-r/small- 
Vy- limit. Then, to first order in ^ and first order in vf_ 

Eceo -E^m- ]-mvl + -L (m^ - El) , (5.28) 
2 2r ^ ' 

and 

Etot -E^{m + Eo$oo) - ^{m + ^o^oo)^' + ^(^0 + + 2mEo$oo). (5.29) 
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These results are interesting but unfortunately confound Newtonian intuition. Ra- 
dial motion of the observers serves to decrease the quasilocal energy measured. 
Specifically the boosted quasilocal energies are equal to their unboosted counter- 
parts minus a kinetic term equal to the | (Total Energy of Fields)f^. The thin shell 
equivalence and Newtonian intuition would lead one to expect the opposite sign for 
this kinetic energy term so this is a bit disturbing. By contrast the no-reference- 
term quantities increase with motion in the expected way. Some discussion of why 
this happens may be found in section where the equivalent effect is considered 
for naked black holes, but briefly the decrease can be thought of as occurring be- 
cause the relativistic effects of the boost compete with those of gravity. Thus, e and 
e begin to converge even as they are both boosted to larger values by the motion. 



Infalling observers 

Next consider the special case where the radially moving observers are falling along 
timelike geodesies towards the gravitational source (be it a black hole or any other 
spherically symmetric matter distribution). Let these observers have started with 
velocity zero "close to infinity" and then have fallen along radial timelike geodesies 
inwards. Rigorously, the geodesic is the one that, with respect to the standard 
time foliation, has radial velocity zero at infinity and —1 (ie. the speed of light) at 
the outer horizon (if the source is a black hole). Now, a test particle starting with 
velocity zero at radial coordinate tq and then allowed to fall towards a black hole 
on a radial geodesic will have coordinate velocity 

±. = -^F(ro)-F(r), (5.30) 

as a function of r, where r is the proper time. Thus an observer infalling on a 
geodesic that was static at infinity will have coordinate velocity ^ = — a/1 — F{r). 
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Let these observers measure time in the natural way (that is = 1), so T" 



-^u" — Y — p^-n". Then the instantaneous radial velocity of the T" observers as 
measured in the static frame is 



V, ^ = -v^r^, (5.31) 
and so the Lorentz factor is 7 = -4^. 

V F 

Substituting this value for 7 into equations ( |5.2CI|J5.2^ , ^.25| ) and making the 



gauge choice $00 = — — so that Et^t doesn't diverge at the horizon, 



Eceo = -r, (5.32) 
Etot = _ ^ , and (5.33) 



E = -rV2 - F. (5.34) 

Note that as r r+ all of these take non-zero values. By contrast Eceo and Etot 
both are zero at r+. Also, keep in mind that for a near extreme black hole, r+ ^ r_. 
Therefore for a black hole that is very close to being extreme, the observers measure 
Efot to have a very large negative value as they approach the horizon. 

Including the reference terms. 



EGeo-E = r{V2-F-l), (5.35) 

and 



Etot -E = r \^V2^ - ^—^j (5.36) 

So near the horizon the infalling gravitational energy (including the reference term) 
goes to {\pl — l)r+ compared to r+ for the static gravitational energy. By contrast, 
the infalling total energy (including reference term) attains arbitrarily large negative 
values as the observers approach the horizon for black holes that are arbitrarily close 
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to being extreme. Static observers however, will measure Etot — E_ = as they 
hover around the horizon. The difference is essentially due to the hugely boosted 
matter terms. The boosting of the geometric terms has a comparatively minor 
effect. 



5.1.3 Z-boosted observers 

Finally consider the slightly more complicated example of a spherical set of ob- 
servers in Schwarzschild space who are boosted to travel "in the z-direction" with 
"constant" velocity v^- In this case, "constant" means with respect to the usual set 
of static and spherically symmetric observers whose four-velocity is u°'. 

Then the four-velocity of the boosted observers is T"' = Nu°' + V"' in the usual 
way where 



1 — f ^ cos^ 9 
1 — vt 
1 



N = ,/ (5.37) 



1 — f ? cos^ 9 



(^-^d^ + v,VEcos9d^^ , and (5.38) 



V" = ^yS£=a?. (5,39) 



■A 



V 



where 9° and dg are the coordinate forms of the vector fields ^ and ^ respectively. 
Note that has been normalized so that the boosted observers measure proper 
time. Then, the static observers see the boosted observers as having velocity v\- = 
Vz COS 9 in the radial direction and f g = sin 9 in the 9 direction. 

Now, taking Qt as a spacelike surface of constant r and t, equation (|5.4| ) along 
with the transformation laws ( ^.321 ) and ( |3.33| ) gives 
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Sill 

Ja = —cT%v^ = ^ ^[d^U (5-41) 

K k[1 — vl cos^ 6) 

where [d9]a is dO in coordinate form. 

Next, I calculate flat-space-embedding reference terms. By embedding condi- 



tions (1-3) on page ^ and taking Vt^ as an r, t constant surface in the reference 
space (same coordinate values as in the Schwarzschild space) the time vector in the 
reference space is = Nv^ + where the lapse is the same as above while 

r = , \ ^ id^ + v,cos9d':),s.nd (5.42) 

a/ 1 — f ^ COS"' d ~ 

r = -^^,&S- (5.43) 

The underlined coordinates are of course in the reference space. Then 

Vz\fF cos6 . 
= . 5.44 

;i+t;2(F-l)cos2^ 

I — COS^ d 



and equation ( [5.6| ) along with the transformation laws ( |3.32| ) and ( |3.33D gives 



2 / l + t;^(F-l)cos^g 

76 = W ^ — and (5.46) 

— nr \ 1 — cos^ a 



ry^ 77 sill / F 

h = ^^%:^h = — 71 2 — I 2fi7 — r\ — rst^^]" (5-4'') 

K K{1 — cos^ 0) V 1 + v;{F — 1) cos'^ 

Not unexpectedly these results are quite a bit messier than the corresponding 
ones for radially boosted observers, and in particular they don't integrate over Qt 
into nice tidy forms. To clear things up a little, I consider a limiting case. For, 
r ^ oo, 

m / — , / \ mv-, sin 6 , 



/ \ 7777? Sin r? 

l-t;2cos2^ and =^ (5.48) 
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Then, integrating over the two-surface of constant r and t the resuh is, 

Eoo = / <f'x^{e — e) 
Jn 

m ( I arcsinf^A , , . 

= - [yi - vl + ' J , and (5.49) 

= v^l - v^m. (5.50) 

Thus the quasilocal Hamiltonian decreases in the same way that the quasilocal 
geometric energy did in the radial boost case. The decrease can again be thought 
of as occurring because of a competition between the relativistic effects of the boost 
versus that of the gravity. Thus e and e converge even as they are boosted. Another 
interesting interpretation of this result can be found in the non-orthogonal paper 
by Hawking and Hunter who considered this case using their Hamiltonian 



method. They interpret the decrease as occurring because some of the energy has 
been transformed into a non-zero gravitational momentum by the boost. 



5.2 Naked black holes 

An interesting application of the quasilocal energy formalism is found in the study 
of the so-called naked black holes. These are low-energy- limit solutions to string 
theory and are characterized by the fact that static observers hovering close to 
their horizons feel only very small transverse tidal forces while infalling observers 
are crushed by arbitrarily large tidal forces. Thus they are naked in the sense that 
even though they are not Planck scale themselves, Planck scale curvatures may still 
be experienced outside their horizons by those infalling observers. Several classes 
of these holes were studied in a couple of papers by Horowitz and Ross ETi but 
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here I will consider only those satisfying the equations of motion ( |2.11| - pTT^ . The 
naked black holes are then a subset of the following class of Maxwell-dilaton black 
hole solutions. The metric is given by 

ds^ = -F(r)dt^ + + RirfidO^ + sin^ Od^^), (5.51) 
F[r) 

where 

F{r) = ^ ^ and i?(r) =r(l- — ) . (5.52) 

\ r J 

In the above, r+ is the radial coordinate of the black hole horizon and r_ is that 
of its central singularity. The accompanying dilaton and electromagnetic fields are 
defined by 

(5,53) 

and 

G 

i,F = -^dtAdr. (5.54) 

These solutions are all magnetic black holes so as discussed earlier the dual form 
of the quasilocal Hamiltonian must be used. The ADM mass and magnetic charge 
are 

G, ^ [{f^^f. (5.56) 

Solving this pair of equations in terms of and r_ one finds that r± = ± 
v/M2 - (1 - a'^)Gl) for a ^ 1 or r+ = 2M and r_ = Gl/M for a = 1. Note that 
for a = these spacetimes reduce to magnetically charged RN black holes. 

Massive near-extreme members of this class of solutions are dubbed "naked". 
To see the reason for this nomenclature note that in terms of the orthonormal 



CHAPTER 5. CLASSICAL APPLICATIONS 



96 



tetrad {li", n'*, 0"} where u"" = l/y/Fd^, n"" = v^9°, = l/Rd^ and c^'* = 
1/ {Rsm.6)d'^ the non-zero components of the Riemann tensor are 

(5.57) 
(5.58) 





F 


TZ 


2"' 




1-FR^ 










^ TZ - ' = 



FR 

— , and (5.59) 
FR FR 



(5.60) 



2R R 

In this section ovcrdots indicate partial derivatives with respect to r (as opposed 
to the time derivatives that they represent elsewhere in this thesis). 

In the alternate inf ailing tetrad {w", n°',9°',(p°'}, where as usual u°' — {l/\)u°'-\- 
rjrf and n° = (1/A)n" + r/w", the non-zero components of the Riemann tensor are 
(in terms of the non-moving components) 

F 

T^unufi T^unun (5.61) 

2 FR 2 

'^u(pu(p T^wfUif ~l~ V ^J^wf>u(p ~l~ '^nif>n((>) '2R ^ R ' (5.62) 

IZfitpnif: — 'Zn(pn(fi ~l~ ^ \'Zu(pu0 ~l~ 'Zn(pn(p) — 2^ ^ ~~R^ (5.6oj 

^<^e<^e is unchanged, 7e^^^^ = 7^ji0s<^>, and 7^^g-^^ = 7^s0n<^■ Clearly if a = then 
R{r) — r and all of the components are the same as for the unboosted frame. 

If a 7^ and 5 = {1 — r_/r+)^^*-^'''" \ then the naked black holes are the subset 
of the above solutions whose parameters satisfy the conditions ^ ^ ^ ^ where 

= R{r+)- That is | » R+ which in turn is much larger than the Planck length. 
Note that if 5 < 1 then r_ ^ r+ and if i?+ » 1 then M, Go > 1- Thus naked holes 
are near-extreme as well as being very large (relative to the Planck length). 
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In the static frame as r — > r+, 



\nunun\ ^ ( 1 - , ) < 1, (5.64) 



1 / 2r 

^"^"^ - ^ (i - (TT^) « 1' ^^-66) 

Thus, all of the curvature components (and consequently the curvature invariants 
calculated from them) are small compared to the Planck scale. 

By contrast, choosing the tetrad to be that carried by the infaUing observers, 
— ^2^2 _ j-F as r — > r+, 

1 

Thus these infaUing observers see Planck scale curvatures. Interpreting these com- 
ponents in terms of the relative acceleration of neighbouring geodesies it is easily 
seen that these observers are laterally crushed by huge tidal forces. 



5.2.1 QLE of naked black holes 

Now, consider a spherical shell of observers falling into a naked black hole. It is to 
be expected that the huge transverse tidal forces will cause the area of the shell to 
shrink at a very rapid rate. Such rates of change of area arc an important factor in 
evaluating the quasilocal energy defined in this thesis. In particular is (up to a 
normalization factor) exactly the (local) rate of change of the area of an infaUing 
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surface of observers. As such it is of interest to calculate the quasilocal energies 
measured by static versus infalling observers and to see how they compare to the 
observed curvatures. As the first step in calculating these energies one finds that 

R 



e — 



AttR"^ 



\/ {r — r_|_)(r — r_), (5.70) 



el = £^"1 = 0, and (5.72) 
£ = -T^- (5.73) 

The gauge choice for the matter term is the same one that was used in the previous 
section. That is, I choose the gauge so that ^4* || Ua, as well as being static, 
spherically symmetric, and non-diverging on the black hole horizon. Though this is 
a long list of requirements, as noted earlier they amount to little more than deciding 
to make the standard gauge choice of electrostatics (or in this case magnetostatics) . 

For the type of infalling observers that were considered in the last section. 



= -^u'^ - \J^n'^ which implies that = -Vl - F and 7 = l/v^F. By 
contrast the joint requirements that T^T_^ = T°'T°' and £f^i3 = £f'^a/3 imply 
that 



\J 1 + R'^il - F)2f - RVl - Frf, (5.74) 
^^^^ and (5.75) 



^1 + R^{1-F) 



Then, 



7 = yi + i?2(l-F). (5.76) 

Egco = -V{r-r+){r-r.)R, (5.77) 
Eoeo = -RR, (5.78) 
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Etot = -R\l-^, (5.79) 
y r — r_ 

Etot = (5.80) 

E = -R, and (5.81) 



E = -i?yi + i?2(l-F). (5.82) 

Evaluating these expressions at r = r+ is straightforward with the only complication 
being 

If <^ 6 then the square of the coupling constant is extremely small even relative 
to 6, and R+ ^1. In fact even if ~ 5 then is of the same order as 1. 
By contrast for ^ 5, R+ ~ T+^T ^ ^- Thus, it is simplest to calculate the 
quasilocal energies for the cases <~ 5 (which includes the magnetic Reissner- 
Nordstrom case for a = 0) and 3> 5 separately. The results along with those for 



oo are displayed in table Note that if a <^ 5 < 1 then i?+ ^ r+ 



From table static observers outside a naked black hole measure Eceo, Etot 
near to the horizon while the infalling observers measure those same quantities to 
be very large. This effect occurs for both 6 <^ and the <^ 6 (which include 
the RN holes) and so cannot be attributed to the "nakedness" of the holes. Of 
course since the reference terms have been omitted, both of these expressions blow 
up if the quasilocal surface is taken out to infinity. 

Including the reference terms, Etot — is very large for static observers near 
to the horizon, where it is R+. It is even larger in the absolute sense for infalling 
observers who measure it as —R+/5. Again however, those effects are seen by 
observers surrounding both naked and near-extreme RN holes and so cannot really 
be attributed to the extreme curvatures. As r — ^ oo the two expressions agree 
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5 >Ri 


5 < 




Quantity 


r 


r ^ r_|_ 


r — >• oo 


—Eoeo 








r 


—Eoeo 


R+ > 1 


^+ 1 
l+a2 S ^ 


r 


Eoeo — E 


R+ > 1 


R+ > 1 


M 


Eceo — E 


Cii?+ > 1 


« 1 


M 


—Etot 








r 


—Etot 


^ » 1 


^ » 1 


r 


Etot — E 


R+ > 1 


> 1 


< < 1 


-(Etot-E) 


^ » 1 


^ » 1 




-1 < < 



Table 5.1: Asymptotic and near horizon values of the quasilocal energies for near- 
extreme dilaton- Maxwell black holes. 5 — (1 — r_/r+)^/*^^"'""^^ ^ 1, — ^ 1 
and <^ 1, where R+ — R{r+). Ci is a constant of the same order as 1. 

which is not surprising since as r — > oo the velocity of the infalling observers goes 
to zero. Note however that this is not the ADM mass. 

More interesting are the measurements of Eqco — E_- If a ~ 1 and the holes 
are large (i?^ » 1) then while static observers near to the horizon measure large 
values, sets of observers faUing into naked black holes actually measure very small 
values for this quasilocal energy. By contrast observers falling into an RN hole 
will measure large values. In fact one can see that if a ~ 1 and R\ 3> 1 then 
these infalling observers will measure Eqco — 1, if and only if the black hole is 
naked. Thus this is an alternate characterizing feature of naked black holes when 
the couphng constant is of a reasonable size. The equivalence is broken if < fa 5 in 
which case the static and infalling observers both measure large energies. Consider 
for example the case where = 6. Then the black hole can still be naked if 5 (and 
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therefore a^) is small enough that S'^R'^ <^ 1. 



5.2.2 Why do naked holes behave this way? 

At the beginning of the previous subsection it was suggested that the curvature 
results could be understood in terms of the rates of change of the surface area of 
shells of infalling observers. In this section the idea is explored in more detail and 
used to provide an explanation of the Eceo — E result. 

First I quantify the expectation that the surface area of a shell of infalling 
observers will be changing extremely quickly as they cross the horizon of a naked 
black hole. Recall that naked black holes are near extreme and so the singularity 
sits "just behind" the horizon (r_ ^ r+). More rigorously, Horowitz and Ross 



p6| noted that an observer passing through the horizon after falling from tq (the 



situation described by equation ( |5.30| )), will hit the singularity at r_ after a proper 



time of At <^ '"t ^ = f"'"'^ . Thus a set of observers infalling on geodesies that 
were stationary at infinity (-F(ro) = 1) will only have a very short time before they 
reach r_. At r_, R{r) — > and so the area of the shell goes to zero. However, by 
assumption R^ ^ 1 and so at the horizon itself, that same area is very large. For 
the area to go from very large to zero in such a small time, one would naively expect 
it to be decreasing very quickly as the observers pass the horizon. This expectation 
can be quantified by using ( p. 31 ) to show that the fractional rate of change of the 



area of the surface Qt as measured by the observers who inhabit that surface is 
^^ 8.J„,J.VF.-l ^_2j,^ 2_( . l\ 

where for the rest of this section primes indicate proper time derivatives. U ^ 6 
(that is, it isn't pathologically small), ^ ~ ^ 1 as expected. By contrast for 
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the RN case (o- = 0), ^ ~ <^ 1. However, the expectation is confounded if 
(5 > 7^ in which case the hole remains naked even while the rate of change is 
more along the lines of the RN values. In that case the extremely small value of 
a suppresses the rapid decrease in area until the observers get even closer to the 
singularity (basically r — r_ ^ a^). 

These rates of change of the area also nicely explain why Eceo — E_is small while 
the observed curvature components are large. Recall that to define the reference 
term E_, Qf had to be embedded into fiat space along with a vector field defined so 
that if Qt was evolved by that vector field and only intrinsic observations were made 
in the resulting timelike three- surface, those observations are identical whether 
they were in the original or reference spacetimes. In particular, the area of Qt 
should change at the same rate. Thus, if the area decreases extremely rapidly, the 
embedded shell of observers in the reference spacetime would have to be moving 
at a correspondingly fast speed. Equation (|5.75|) quantifies this saying that Vy = 



R/ y \ + B? at the horizon. Then for a^~l,-R^l^t;|_~l and the observers 
would have to move at close to the speed of light in the reference time to match 
the rate of change of the area. By contrast, for ^ 0, -R ~ 1 ^ ~ |. The area 
is changing at a relatively leisurely rate so the observers would not need to move 
so fast in the reference time. 

For observers moving at extremely rapid velocities there is a sense in which the 
relativistic effects of their speed become more important than those due to gravity. 
To see this recall equation ( 4.3.2[ ) where it was shown that — e^"^ is a constant 



independent of the speed of the observers. Now, by construction is the same 
in both the reference and original spacetime and so the geometric QLE can be 
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rewritten as, 

EGeo-E = j rf2xv^(vi^Til2- A/^2^il2j . (5.85) 

If el is much larger than e and e, and so in a sense the relativistic effect of speed 
dominates over that of curvature, then at the horizon 

and so as becomes larger and larger the observed quasilocal energy becomes 
smaller and smaller. Physically, though e and e_ are boosted to be very large, the 
difference between them simultaneously becomes very small. In particular for naked 
black holes 

^Geo - I - (I) = ^ = (5-87) 

and it can be seen that in this case the geometric quasilocal energy is actually the 
inverse of the (normalized) rate of change of the area. As noted in section |5.1.2 
these general ideas explain the much less dramatic decrease of the quasilocal energy 
for boosted observers in the Reissner-Nordstrom spacetime as well. 

By contrast Etot — E_ includes matter terms which are also boosted to be very 
large. There is no corresponding term in the reference spacetime to cancel these 
large terms out. The result is that the matter terms dominate over the geometrical 
terms in Etot — K and so this total quasilocal energy is very large. 

5.3 Tidal heating 

As a final classical calculation I use the quasilocal formalism to calculate the amount 
of work done by an external gravitational field when it deforms a self-gravitating 
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body. The canonical example of this effect in our own solar system is found in the 
gravitational interactions between Jupiter and its moon lo. In that instance, the 
gradient of Jupiter's gravitational field distorts the shape of lo away from being a 
perfect sphere and then tidally locks it in its orbit so that it always presents the 
same face to Jupiter. That orbit is strongly perturbed by the other Galilean moons 
and so its radial distance from Jupiter varies with time. With this variation comes 
a corresponding one in the gradient of the field and so lo is gradually stretched and 
then allowed to relax. The energy transferred by this pumping is largely dispersed 
as heat and it is this heat that produces the volcanic activity on lo. The same type 
of process occurs in principle for any two bodies in non-circular orbits about each 
other. 

To calculate the gravitational energy transferred to lo during this process using 
the quasilocal formalism, I'll need a metric describing the situation. To this end, 
first consider the situation from a Newtonian perspective. Assume that the self- 
gravitating body is far enough away from the source of the external field that 
that field is nearly uniform close to the body. Then in a rectangular coordinate 
system that orbits with the body (with origin fixed at the center of mass), the 
Newtonian potential of the external field may be written as ^^xt = ^abx"'x^ where 
Sab is the (time-dependent but symmetric and trace-free) quadrupole moment of 
the field and is the position vector based at the body's centre of mass. At 
the same time, to quadrupolar order the Newtonian potential of the body is $0 = 
-M/r - (3/2)r~3JafenV, where M is the mass of the body, r is the radial distance 
from the centre of mass. Tab is its (time-dependent but symmetric and trace-free) 
quadrupole moment, and n°' = x°'/r is the unit normal radial vector. 

From this description one can use the techniques of Thorne and Hartle 
to construct a metric that describes these situations in the slow moving, nearly 
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Newtonian limit. First, define an annulus surrounding the self-gravitating body 
whose inner boundary is chosen so that its gravitational field is weak throughout 
and whose outer boundary is chosen close enough so that the external field is nearly 
uniform. This region is termed the buffer zone. The rectangular coordinate system 
from the Newtonian limit is replaced with one that is chosen so that the metric 
is as close to Minkowskian as possible over the buffer zone. Then to first order in 
perturbations from Minkowski and first order in time derivatives the metric can be 
written as 



rfs2 ^ -{I + 2^)dt^ + 2{Ab + dt^b)dx''dt 

+ [(1 - 2^)Sab + + db Qdx'^dx'' (5.88) 

where the indices run from one to three and 5ab = diag[l, 1, 1] is the Cartesian metric 
on a spacelike slice. The Newtonian potential is still $ = —M/r — {l/2){3r^^Tab ~ 
r'^£ab)n"'n^ while 

is a vector potential that must be included so that the metric is a solution to the 
first order Einstein equations. Here, is the radial normal with respect to the 
fiat spatial metric 5ab and r"^ = + + . The diffeomorphism generating vector 
field ^b represents the gauge ambiguity in setting up a nearly Minkowski coordinate 
system. In order that the metric be slowly evolving and nearly Minkowski, ^b must 
be of the form 



ib = -^Ti^cn" + Pr^Sbcu' + ir^Ecdn^n^nb, (5.90) 



a 
r 

where a, P, and 7 are free constants of order one. 



To measure the flow of quasilocal energy, I define i? as a surface of constant 
r surface in the buffer zone, foliate it with constant t spacelike two-surface Qt, 
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and define the time vector as d/dt. Then I can calculate Ht from equations 
(|3.21| ) and (|3.22|) . As I mentioned in section p.5| , I will neglect reference terms here 



because for a wide range of choices of how to define them, they don't contribute in 
a situation such as this where I am calculating rates of change. Of course this also 
serves to simplify the already messy calculations. 

In calculating the time rate of change it is most convenient to switch to spherical 
coordinates. Making the standard transformation x° = r[sin 6 cos 0, sin 6 sin 0, cos 6], 
the metric becomes 

ds'^ = -(1 - 2^)dt^ + 2Ardrdt + 2rAgdedt 
+2r sin OA^dcpdt + (1 + 2$ + Hrr)dr'^ 
+r\l + 2$ + Hee)de^ + sin^ ^(1 + 2$ + H^^)d(f)^ 
+rHrgdrd9 + r sin 9Hj.(j,drd(f) + sin OHg^dOdcl), (5.91) 



where 



4a 



Hrr = -^J„ + 6(/5 + 7)rXr, (5.92) 
2a. 



Hge = ^l0e + 2f3r^See + 2-fr^Srr, (5.93) 



2a, 



H^^ = —I^^ + 2f3r'S^^ + 2^r'Srr, (5.94) 



Hre = -^Ire + {4:P + 2^y£re, (5.95) 
Hr^ = --X,^ + (4/3 + 27)r2£,^,and (5.96) 

He^ = '^Ie^ + 2pr^ge<t>. (5.97) 



In these expressions £„■ = i^afeC^eJ, Sro = ^ab^-r^-g, etc., with = Cg = d^e'^ and 
= (1/ sm9)d^e'^. Also, = {Ab + dt^bje^ri ^tc, but their expanded forms are 
not needed since only time derivatives of them show up in later calculations and 
the calculation is only been done up to first order in time derivatives. 
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As might be expected the subsequent calculations are quite involved and I did 
them with a lot of help from the GRTensor |]75| package for Maple. Ultimately 
though after a huge amount of algebra, equation (|3.22|) works out to become 

H = -\ I d'^X^Tl'^HTlab (5.98) 



2 



- ^-E i 

+^ |^(-3 -2/3- 2/3^ + 47 + 47^ + 8/3^)£ab£ab 
+ ^ 1^(3 - 2a + 6/3 - 127 + 8aj)SabIab 

Note that repeated indices continue to indicate summation. Since Sab and lab are 
Cartesian tensors, the index position doesn't matter. These calculations used the 
identities 



d6 d(f) sin 9 Arr B rr = (Stt/ 15) AabBab and (5.99) 
J ded<j>sm9{2Ae^Be^ - AeeB^^ - A^^Bee) = (47r/3)A,bS,fe, (5.100) 

where the integrations are over the unit sphere. 

This result requires some interpretation. As the external field changes with time 
and thereby forces the self-gravitating body to change configuration, the work done 
by the external field can be split into time reversible and irreversible parts (as seen 
in equation |5.98| ). The reversible part represents work being done to increase the 



potential energy of the system and is recoverable. On the other hand the irreversible 
part represents work being done to deform and/or heat up the system. This is the 
tidal work that 1 am interested in and by the above it is {l/2)SabIab, which is 
the same leading term obtained when one does the corresponding calculation in 
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Newtonian gravity or with pseudo-tensors |^ . It is completely independent of the 
diffeomorphisms generated by which correspond to fluctuations of the quasilocal 
surface. 

Note that a gauge ambiguity similar in form to (though not identical with) 
that found in the time reversible term is also found in the corresponding results 
obtained by the Newtonian and pseudo-tensor methods. What is much more clear 
in this calculation however, is that the ambiguity is a result of fluctuations of the 
quasilocal surface through the flelds as generated by the diffeomorphisms. Keep 
in mind that those other methods also give answers with time reversible and time 
irreversible parts so that is not unique to the quasilocal procedure but instead is a 
physical property of the system as I argued in the previous paragraph. 

Finally for completeness consider how the energy flow splits up into its compo- 
nent parts as considered in equation ( p. 211 ). In the approximation in which I am 



working, the angular momentum term is zero and what is left are the two terms 

taab- It can be shown that 

Hn = ^Sabiab + -^^ablab — ^^ablab ^^ab^ab (5.101) 

d / 47 + /3-2 ^ 1 2a-3 

+ — S — T" C-abOab — T^C-ab^ab ^ab-lab 



dt { 30 10 20r 

The second term is a bit more complicated. It is 

Ha = — — Sablab + ^Sab^ab + -^Sablab (5.102) 

d ( 

+ ^|4(3-« + 3/9-67 + Aa-f)SabIab 
dt y 15 
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Thus part of the work done is measured by deformations of the surface and part 
is measured by changes in how observers choose to measure the rate of passage 
of time. Note that individually the time irreversible sections of the two parts 
are gauge dependent but when one combines them equation (|5.98|) returns and the 
gauge dependence vanishes back into the reversible part where it would be expected. 

There are two ways to look at this calculation of tidal heating. The first is to see 
it as an astrophysical application of the quasilocal energy and so an alternate way to 
calculate the tidal heating effects. As I have argued above, it has an advantage over 
previous methods of calculating the magnitude of the effects in that the source of 
the gauge ambiguity in the final result can be clearly identified. It is also somewhat 
tidier than the corresponding pseudo-tensor methods since the integrals are defined 
in terms of tensor quantities and so are covariant. On the other hand, the second 
way to look at the result is as a check on the physical relevance of the Brown- York 
energy. That it can reproduce the results produced by other methods is a good 
argument for its physicality. 

On the down side, I haven't shown that this result is independent of the exact 
choice of the form of the reference term. For example, it would be good to show 
that the final results would be the same with the two-surface embedded in 4D 
reference term. Further, from the work of section |5]^, one is led to think that it is 
the geometric quasilocal energy that is the physically relevant quantity. Here I have 
calculated the Hamiltonian based on a physically arbitrary coordinate time vector. 
However, to resolve either of these questions would require extensive calculations 
so for now I let the result rest in its computationally simplest form that I have 
considered here. 



Chapter 6 

Quantum creation of black hole 
pairs 



While the previous chapter considered apphcations of the quasilocal Hamiltonian 
in classical general relativity the current chapter will consider its application to 
semi-classical quantum gravity. Specifically, I combine it with the path integral 
formulation of quantum gravity to calculate the probability that a pure deSitter 
spacetime will transform itself into a pair of charged and rotating black holes in a 
deSitter background via a quantum tunneling process. This work was published in 



As a short outline, the section |6.1| reviews path integrals as applied to quantum 
gravity and then the following sections fiesh out that introduction as applied to 
the case of black hole pair creation in a deSitter background. Section |ti.2| examines 



the classical description of spacetimes containing pairs of black holes. Section ^ 
constructs the instantons used to mediate the creation of such spacetimes while 
section uses the Brown- York formalism to decide which is the correct action to 
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use in the path integral calculations. Finally section |6.5| evaluates those integrals 



to lowest order and section looks back on some questions that arose during the 
calculations. 



6.1 The idea 

A standard problem of quantum mechanics is to calculate the probability that a 
system passes from an initial state Xi to a final state X2. If the classical equations 
of motion for that system can be derived from a Lagrangian action I then the path 
integral formulation of quantum mechanics provides a prescription for calculating 
the probability amplitude that that transition occurs. Basically it says that one 
should consider all conceivable "paths" F that the system could follow to evolve 
between Xi and X2 (and not just those that satisfy the classical equations of mo- 
tion). If one calculates the action /[F] for each of those paths then the probability 
amplitude that the system will move from state Xi to state X2 is hypothesized to 
be given by the path integral 

^12= f d[r]e-''^^\ (6.1) 



where the integral is over all possible paths. Note that I use the word "hypothe- 
sized" above because in general, this integral is not well defined and so the path 
integral methods are sometimes more of a way thinking about these problems rather 
than actually calculating exact amplitudes. A more complete description of the ap- 



proach can be found in [35 



Despite problems of definition, the procedure was generalized to a formulation 
of quantum gravity in the 1970's (see for example |^). The philosophy behind the 



approach remains the same but the details change quite a bit. 
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In the first place, it is no longer a trivial matter to define an instantaneous 
configuration of a system if that system is a general relativistic one. For a system 
with gravitational and electromagnetic fields (the case in which I'll be interested 
in this chapter) an "instant" will be defined as it was in chapter 0. Namely it will 
consist of a three-manifold S with Riemannian metric hab, conjugate momentum 
density P"'' (or equivalently extrinsic curvature Kab) describing how the system is 
evolving at that "instant" , and vector field densities S"' and defining the electric 
and magnetic fields on E. These four fields must satisfy the constraint equations 
(12:241 ), ( ^^261) , i^M) , and (|^ and if they do, the "instant" can be embedded 
in a larger four-dimensional solution to the Einstein-Maxwell equations. In fact if 
S is a Cauchy surface then it uniquely determines that solution via the evolution 
equations ( ^1251) , ( ^1271) , and (|^. 

Then, using the path integral approach one must consider all possible interpo- 
lations (or "paths") between the states (not just those that would be allowed by 
the classical evolution of the system). This means considering four- manifolds (with 
boundaries) Mu, along with metric fields Qajs and electromagnetic field tensors F^/? 
on those manifolds such that the surfaces Si and S2 and their accompanying fields, 
may be embedded in M12 and its accompanying fields Q. I reiterate that the space- 
time paths {M12, QajS, Fafs) are not, in general, solutions to the Einstein-Maxwell 
equations. 



^In this context, a three manifold E and its accompanying fields {hab, P'^^,£"',B'^} is said to 
be embeddable in the spacetimes (Afi2, Fafj) if there exists an embedding (in the differential 
topology sense), $ : E ^ M12 such that $*(/iafc) = h^f^l^, $*(P«'') = P'^^]^, $*(£,) = £^\^ = 
—2\fh/ nFapuf^ , and$*(Sa) = S^lj^ = ~2Vh/ K^e^^^^ F^su^ ■ In the preceding $* represents 
the appropriate mapping as derived from $ for each quantity. 
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Next, the action 



/[Mi2, gap, Fap] = / d^xy/^{n - 2A - FapF"'^) + (boundary terms) , (6.2) 

for each path must be calculated, where the integration is over all of M12 between 
the two embedded surfaces Si and S2, and the boundary terms are calculated on 
the boundaries of M12 that are consistent with the boundaries of Si and S2. How 



an appropriate action functional can be picked will be discussed in section pA . 

Finally, the value of the action for each path is used to assign a probability 
amplitude for that path. The amplitudes are summed over all of the possible paths 
to give a net probabihty amplitude that the system passes from Xi to X2. This 
summation is represented as a functional integral over all of the possible manifold 
topologies, metrics, and vector potentials Aa (generating the field strength Fa/3) 
interpolating between the two surfaces. That is, 

^12 = y"rf[Mi2]rf[(?]rf[A]e-^^[*^"'3'^l (6.3) 

Thus at least in principle, the probability that a spacetime initially in a state 
(S,/^afe,P''^£^i3")l passes to a state (S, f i3")2 is proportional to |^i2p 
(the wave function hasn't been normalized). Unfortunately the integral ( |6.3| ) cannot 
be directly calculated. In the first place, there is no known way to define a measure 
for the integral. Second, even if such a measure were known, it seems quite likely 
that calculation of the integral would be impractical, considering that the parameter 
space of paths from Xi to X2 has an uncountably infinite number of dimensions. 

Fortunately there is a well-motivated simplifying assumption available. In anal- 
ogy with flat-space calculations, it is argued that to lowest order in h, the 
probability amplitude may be approximated (up to a normalization factor) by 



'i2^e-'% (6.4) 
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where Ic is the real action of a (not necessarily real) Riemannian solution to the 
Einstein-Maxwell equations that interpolates between the given initial and final 
conditions. Essentially, it is assumed that such a solution is a saddle point of 
the path integral. This solution (if it exists) is referred to as an instanton. The 
probability that such a tunnelling occurs is then proportional to |\E'i2p ~ e~^^'=. 
Note that this interpretation requires that the action be real and positive, and 
ideally that all of the fields on its boundary match those in the Lorentzian solution 
"instants" so that it can smoothly match onto that solution. As will be seen in 
section |6.3| this is sometimes a bit much to ask for, but if one only requires a match, 
rather than a smooth match, it can be done. 

As an alternative to paths and instantons interpolating between two spatial 
slices El and S2, one can consider those with a single spacelike boundary that 
match onto a single slice labelled S2. In that case one can interpret the resultant 
path integral as calculating the probability for the creation of the three-space S2 
from nothing and the initial boundary condition is the no-boundary condition of 
cosmology |^5|. One can then compare the relative creation rates for different space- 
times (eliminating the need to calculate a normalization factor) and even interpret 
those probabilities as giving the chance that the different spacetimes tunnel into 
each other |]13[. This is the approach that will be taken here. 

Finally, before passing on to consider the classical solutions that describe the 
spacetimes that I want to create, note that path integrals (especially in the single 
boundary case) can be interpreted as sums over all the possible histories of the 
system being considered [^] and in particular this interpretation is often carried 



over into gravity [^. Then the path integral can be interpreted as a thermody- 
namic partition function and so this formalism naturally lends itself to the study of 
gravitational thermodynamics. As was discussed in and I will consider to some 
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extent in section |6.4| , the choice of the action I will determine the exact partition 
function being considered - that is the canonical, microcanonical, or grand canonical 
partition functions. Given this correspondence the terminology of thermodynamics 
will sometimes be used in the following. Ultimately I will also use the connection 
to extract some conclusions about black hole entropy from my calculations. 



6.2 Accelerating and rotating pairs of black holes 

Since I am interested in calculating the creation rate for a pair of black holes 
accelerating away from each other in a cosmological background, the first step in 
the path integral calculation discussed above is to find a solution to the Einstein- 
Maxwell equations that describes such a physical situation. Such solutions are the 
subject of this section. 



6.2.1 The generalized C-metric and KNdS spacetime 

The well-known C-metric solution to the Einstein equations (first interpreted in 



||62|| ) describes a pair of uncharged and non- rotating black holes that are uniformly 



accelerating away from each other. In this metric was generalized to allow the 



holes to be charged and rotating, as well as to allow the inclusion of a cosmological 
constant and NUT parameter. 

In general, spacetimes of this type contain conical singularities. Physically these 
arise if the rate of acceleration of the black holes does not match the energy source 
available to accelerate them. Thus, in the cosmological case, if the black holes are 
accelerating faster or more slowly than the rest of the universe, conical singularities 
will exist. Physically, these may be interpreted as cosmic strings or "rods" that are 
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pulling or pushing the black holes apart (or together) to make them accelerate faster 
(or slower) than the rate of expansion of the universe as a whole. The singularities 
arc eliminated if the acceleration of the holes is matched to the amount of energy 
that is available to accelerate them. In that case no extra acceleration is required 
and so the cosmic strings or rods aren't needed to provide the extra energy. 

The generalized C-metric takes the form 

ds^ = -J— i '^'''^ + TW (rf- I'drf \ 

with accompanying electromagnetic field defined by the vector potential 

_ eoq{dT +p'^da) gop{da - q^dr) 
l+p^q^ 1+p^q^ 

where p, q, r, and a are coordinates, 

P{P) = {— - 5o + 7) + 2np - ep2 + 2mp3 + {— - - 7)/, (6.7) 

and Q{q) = P{q) + 1^(1 + ^"^). A is the cosmological constant, 7 and e are constants 
connected in a non-trivial way with rotation and acceleration, eg and Qq are linear 
multiples of electric and magnetic charge, and m and n are the respectively mass 
and the NUT parameter (up to a hnear factor). This solution can be analytically 
extended across the coordinate singularity at p = g, so that on the other side 
oi p = q there is a mirror image of the initial solution (though with opposite 
electric/magnetic charge and direction of spin). Thus, if one views it as describing a 
pair of black holes, the two holes will be on opposite sides of that p — q hypersurface. 

In general this metric has a conical singularity in the (p, a) hypersurface which 
corresponds to the above mentioned string or rod. There are a few limiting processes 
that can be used to remove this singularity, but on setting the NUT charge to zero. 
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at least one of them reduces the metric to the Kerr-Newman-deSitter metric. Details 
of that process can be found in appendix [B.l| . In Boyer-Lindquist type coordinates, 
the KNdS metric takes the form iTSl 



O G 

ds^ = --^ (dt-asm'^ed(f)Y + ^dr'^ (6.8) 

yx Q 

G ,„2 Hsin^ ^ / ,, r 2 21 , a2 

where ^ = + cos^ 9,H = l + (A/3)a2 cos^ e,x^ = ^ + (A/3)a^ and 

Q ^ ^a"^ - 2Mr + {a' + El + Gl) . (6.9) 

The individual solutions are defined by the values of the parameters A, a, M, Eq, 
and Go which are respectively the cosmological constant (since I'm interested in 
deSitter type spacetimes, assume that it is positive), the rotation parameter, the 
mass, and the effective electric and magnetic charge of the solution. Along with 
the electromagnetic field 

F = {Xdr A{dt-a sin^ 9d(f)) + Y sin 9d9 A {adt - (r^ + a^)d(p) } , (6.10) 

where X = EqT + 2aGor cos 6, Y = CqT — 2aEor cos 6, and F = — cos^ 6, this 
metric is a solution to the Einstein-Maxwell equations. For reference note that a 
vector potential generating this field is 

A = ^{dt-a sin^ ed(j)) + {adt - (r^ + a^) d(f)) . (6.11) 

^ A, ^ A. 

Keep in mind however the restrictions against dyonic spacetimes that were discussed 
in previous chapters. Thus, even though this is a dyonic solution I'll only be able 
consider the creation of spacetimes where either Eq = Q oi Gq = 0. 

The roots of the polynomial Q correspond to horizons of the metric. As a 
quartic with real coefficients, Q may have zero, two, or four real roots. I will be 
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interested in cases where there are four real roots, three of which are positive. In 
ascending order, these positive roots correspond to the inner black hole horizon, 
the outer black hole horizon, and the cosmological horizon. 

If all of the roots of Q are distinct, then by the standard Kruskal techniques the 
metric may be analytically continued through the horizons to obtain the maximal 
extension of the spacetime [Q. Though this maximal extension is infinite in extent, 
a variety of other global structures are possible if periodic identifications are made. 
In particular, demanding that there be no closed timelike curves in the spacetime 
and also that there are two black holes in spatial cross-sections of constant time 
coordinate t, the global structure is uniquely determined and is shown in figure 
|6ni (for a two-dimensional constant (j), 9 = ^ cross section). As indicated the 
figure is repeated vertically and periodically identified horizontally, r = is the 
cosmological horizon, r = is the outer black hole horizon, and r = is the 
inner black hole horizon. The wavy lines at r = represent the ring singularity 
found there for a 7^ 0. If a = then this singularity may not be avoided and the 
spacetime cuts off at r = 0. Otherwise the singularity may be bypassed and one 
may proceed to negative values of r. r = r_ is the (negative) fourth root of Q. 
The constant t spatial hypersurfaces are closed and span the two black hole regions, 
cutting through the intersections of both the r = and r = Tq lines. The matching 
conditions are such that, in the spatial hypersurfaces, the two holes have opposite 
spins as well as opposite charges. Thus, the net charge and net spin of the system 
are both zero. Note that it is not possible to periodically identify the spacetime so 
that the spatial sections contain only a single black hole. 
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Repeat 




Repeat 



Figure 6.1: The global structure of the KNdS solutions with periodic identifications 
so that t = constant hypersurfaces contain only two black holes. 

6.2.2 The allowed range of the KNdS solutions 

The allowed ranges of the parameters so that Q has three non-negative roots are 
calculated in appendix |B.2| and shown in figure |6.2| . The parameter space is the 



region bounded by the two solid sheets plus the = 0, M = 0, and Eq + = 
sheets. The darkest sheet corresponds to the extreme black hole for which the inner 
and outer black hole event horizons are degenerate and the lighter gray sheet is the 
case where the outer black hole horizon is degenerate with the cosmological horizon. 



Taking nomenclature from the non-rotating instantons discussed in [71 1 I'll call the 



CHAPTER 6. QUANTUM CREATION OF BLACK HOLE PAIRS 



120 



extreme black hole case the cold KNdS spacetime while the second will be the Nariai 
KNdS spacetime. I'll denote their intersection the ultracold KNdS spacetime. The 
transparent sheet represents a special case of solutions corresponding to lukewarm 



spacetimes, which will be discussed in subsection 3.2.3 



Note that the extreme cases, though limits of the KNdS metric, have differ- 
ent global topological structures. In fact the Nariai and ultracold spacetimes do 
not even contain black holes. Their metrics in coordinate form may be found in 
appendix [B.2| , but here I'll just comment briefly on some of their properties. 

In the cold case the double horizon of the black hole recedes to an infinite proper 
distance from all other parts of the spacetime (as measured in a spacelike surface 
of constant t). Thus, the global structure of the spacetime changes. In particular, 
the region inside the black hole is cut off from the rest of the spacetime. Making 
appropriate periodic identifications of the global structure so that the spacetime 
contains two (in this case extreme) black holes, the structure is shown in figure 



0| . In that figure opposite sides of the rectangle are identified, r = rc is the 
cosmological horizon and r = Tg^i is the double black hole horizon. If a = 0, then 
the spacetime cuts off at the singularity at r = 0. Otherwise, one may pass through 
the ring singularity to the negative values of r, including r_, the fourth root of Q. 
Note that in this case, the t = constant hypersurfaces contain two extreme black 
holes, and so are not closed as they are in the regular KNdS spacetime. The metric 
for this case is given in appendix p. 2.1 . 



As noted the Nariai solution shown in |6.4| is no longer a black hole solution, 
and there is no longer a singularity at finite distance beyond either of the horizons 
at p = ±1. In fact, the diagram is the same as that for two-dimensional deSitter 
space. If there were no rotation (a = 0), then this spacetime would just be the 
direct product of two-dimensional deSitter space, and a two-sphere of fixed radius. 
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Figure 6.2: The allowed range of the KNdS parameters. The range is bounded by 
the planes M = 0, = 0, Eq + = 0, the cold solutions (the darkest sheet) 
and the rotating Nariai solutions (the lighter gray sheet) . Also shown as a meshed 
sheet are the lukewarm solutions. 
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t Repeat t 




j Repeat j 



Figure 6.3: The Penrose-Carter diagram for a two hole cold KNdS spacetime. 



With rotation, of course the situation is not so simple but if a = 0, it reduces to the 
non-rotating charged Nariai solution considered in [^]. The metric may be found 
in appendix |B.2.2. 



Even though the Nariai solution is not a black hole solution itself, it was shown 



m 



43 1 that an uncharged, non- rotating Nariai solution is unstable with respect to 



quantum tunnelling into an almost-Nariai Schwarzschild-deSitter spacetime. It is 



usually argued that this tunnelling carries over analogously with the inclusion 
of charge and rotation, in which case Nariai solutions decay into near Nariai KNdS 
spacetimes. Thus, in the future sections where I study black hole pair creation this 
solution will remain of interest, as a route to black hole pair creation will be to 
create a Nariai spacetime and then let it decay into a black hole pair. 

A similar argument can be made [|7I|] for the ultracold spacetimes found at the 
intersection of the parameter spaces of the cold and Nariai solutions. There are 
two possible spacetimes (appendix p.2.3| ), one with one horizon and the structure 
of Rindler space and the other which is conformally Minkowski and has no horizons. 
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Figure 6.4: The Penrose- Carter diagram for the Nariai hmit spacetime. 



Neither contain black holes. 



6.2.3 Issues of equilibrium 

Before passing on to the next section where instantons to create the above space- 
times will be constructed, I'll pause to examine whether these solutions to the 
Einstein-Maxwell equations are stable with respect to semi-classical effects. This 
is relevant because traditionally one only considered the quantum creation of black 
hole pairs in thermodynamic equilibrium, as it was thought that these were the 
only cases where regular instantons could be constructed. It is not so clear to- 
day (see for example ^) that that requirement must be enforced, but since 



I'll be more-or-less using the traditional methods here and also try to draw some 
conclusions about the thermodynamics of the spacetimes, it is an issue that must 
be considered. To check for this equilibrium, one must consider three phenomena: 
thermodynamically driven particle exchange between the horizons, electromagnetic 
discharge of the holes (due to emission of charged particles), and spin-down of the 
holes (due to emission of spinning particles and super-radiance). 

It is well known that a black hole emits particles in a black body thermal spec- 
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trum and thus may be viewed as having a definite temperature |^8|. In the same 
way, it has been shown that deSitter horizons may also be viewed as black bodies 



and have a definite temperature [^]. For a spacetime with non-degenerate hori- 
zons, these temperatures may be most easily calculated by the conical singularity 
procedure (which will show up again in the next section during the instan- 
ton construction). First, corotate the coordinate system with the horizon whose 
temperature is being calculated. Second, analytically continue the time coordinate 
to imaginary values. For definiteness label the imaginary time coordinate T, the 
radial coordinate TZ, and let the horizon be located at TZ = TZh- Next, consider a 
curve in the T — TZ plane with constant radial coordinate TZ = TZq. Periodically 
identify the imaginary time coordinate with some period Pq so that this curve be- 
comes a coordinate "circle" and may be assigned a radius Rq and circumference Co 



according to the integrals 



no 



Ro = \ I Va-R.nd'JZ 



, and Co 

r=o 



. (6.12) 

7^=7^o 



Finally, calculate limT^^^T^^ Pick the value of Pq so that the limit has value 27i. 
Then, the horizon has temperature Th = 1/Po, and surface gravity Kh = 2tx / Pq. 

If there is a degenerate horizon as is the case for a cold black hole, then that 
horizon is an infinite proper distance from all non-horizon points of the spacetime. 
In such a situation there is no restriction on the period with which the degenerate 
horizon can be identified, and it has been argued that the black hole can 
therefore be in equilibrium with thermal radiation of any temperature. 

Now consider which of the spacetimes are in thermodynamic equilibrium. First, 
consider the general non-extreme KNdS solutions. The temperature of the outer 
black hole horizon and the cosmological horizon are respectively. 



dO 



dO 



CHAPTER 6. QUANTUM CREATION OF BLACK HOLE PAIRS 



125 



There are two ways that these two temperatures may be equaL The first is if 
f'bh = fch which actually is the Nariai spacetime. However there is also a non- 
extreme solution labelled the lukewarm case. Its parameterization is considered in 
appendix B.2.4 . 

The cold limit is in thermodynamic equilibrium at the temperature of the cosmo- 
logical horizon, for as has been noted an extreme black hole may be in equilibrium 
with thermal radiation of any temperature. As noted, the Nariai limit too is in 
thermodynamic equilibrium, with both horizons having the same temperature 

TNar = (6.14) 

47r 

where e and 6 are defined in appendix |B.2.2| . The first ultracold case has only one 
horizon with temperature 

Tuci = ^, (6.15) 

and so with no other horizon to balance this one off, it is not in thermal equilibrium. 
The second ultracold case has no horizons, and so is trivially in equilibrium. 

Next consider discharge of the black holes. Even if the black hole and cosmolog- 
ical horizon are in equilibrium with respect to thermal particle exchange between 
them, there can still be a net exchange of charge between the horizons. The mecha- 
nism is that even though both may create the same number and masses of particles, 
an excess of charged particles will be created at the black hole horizon, and so it 
will discharge |^0|. This effect can be quite rapid and so in most charged 
black hole cannot be said to be truly in equilibrium. However, there are a couple of 
ways to avoid the discharge. If there are no particles of the appropriate charge that 
are also lighter than the black hole then discharge cannot occur. Thus, if magnetic 
monopoles do not exist then the magnetic holes will be stable with respect to dis- 
charge. Further, even if the appropriate light charged particles exist, the discharge 
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effects will be small if the temperature of the black hole is small relative to the 
mass of those particles. That is, the more massive the black hole, the slower the 
discharge. 

Finally consider the spin-down of the black holes. If the black holes and cos- 
mological horizons are at the same temperature, then there will be no net energy 
exchange between the horizons, but the particles created at the black hole hori- 
zons may still have an excess of angular momentum relative to those created at 
the cosmological horizons. Unfortunately this effect is not as well studied for cos- 
mological spacetimes as is the equivalent discharge case. Still, from the extensive 
calculations in asymptotically flat space |]77| , [76| , [2^ one can say the following. In 



fiat space, the direct spin-down by particle creation is a relatively slow process but 
it is greatly amplified by super-radiance. In combination the two processes cause 
angular momentum to be radiated relatively more quickly than mass is radiated 



unless there are a truly ridiculous number of scalar fields in the spacetime [P4 



Preliminary indications ||8^ are that spin-down occurs at least as quickly and pos- 
sibly more quickly in asymptotically deSitter spacetimes which are not in thermal 
equilibrium. There aren't any corresponding calculations for black holes which are 
in thermal equilibrium with cosmological horizons, but that said, in the other cases 
spin-down is a relatively quick effect which means that a rotating black hole in 
deSitter space probably cannot be thought of as being in full equilibrium. Possibly 
the presence of thermal equilibrium might cause something miraculous to happen, 
but that is unlikely and in any case a matter to be resolved by future calculations. 
However, even in the absence of such a miracle, the physically intuitive notion that 
a black hole that is rotating slowly relative to its mass will discharge slowly is sup- 
ported by the existing results, and so it seems likely that at least a class of these 
holes may be considered quasi-static in a thermodynamic sense. 
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That said, in the following section I'll show that only thermal rather than full 
thermodynamic equilibrium appears to be required for the construction of smooth 
instantons. Of course if the created system is not in full equilibrium one can- 
not really draw conclusions about its thermodynamics. Thus, the reader who is 
uncomfortable with the quantum creation of spacetimes that are not in full ther- 
modynamic equilibrium can consider all of the following to apply only to the subset 
of spacetimes that are at least quasi-static. 

6.3 Instanton assembly 

In this section I construct the instantons that will be used to study the creation of 
the spacetimes considered in the previous section. As discussed in the review of the 
path integral formalism, these instantons must both be solutions to the Einstein- 
Maxwell equations and also should match as smoothly as possible along a spacelike 
hypersurface onto the spacetime that they create. The instantons constructed here 
will satisfy the cosmological no boundary condition, and so I will not need to worry 
about matching to initial conditions. 

6.3.1 Analytic continuation 

For static spacetimes, the first step of instanton construction is usually to analyti- 
cally continue t it. For a static spacetime expressed in appropriate coordinates, 
this gives a real Euclidean solution to the equations of motion but for a spacetime 
that is only stationary it will usually produce a complex solution to the equations 
of motion. For now 1 accept this complex solution but at the end of this section 
I'll consider its relative merits compared to the more standard approach where 
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other metric parameters are also analytically continued in order to obtain a real 
Euclidean metric. That said, I proceed in the following manner (which is equivalent 
to continuing t it). 

If a spacetime is foliated by a set of space-like hypersurfaces labelled by a 
time coordinate t, the most general Lorentzian metric can be written as 

ds'^ = -N^dt^ + hab{dx'' + V''dt){dx'' + V''dt) (6.16) 
= {-N^ + habV''V^)dt^ + 2habV^dx''dt + habdx''dx\ 

where as usual hab is the induced metric on the hypersurfaces, is the lapse 
function, and V"" is the shift vector field. Using the prescription of [T^, the ana- 
lytic continuation can be made by making all of the Lagrange multipliers from the 
Hamiltonian purely imaginary. To wit, I start by changing the lapse and shift so 
that — >■ iN and V — > iV"'. The spacetime metric for the proto-instanton then 
becomes 

ds^ = {N^ - habV''V^)de + 2iKbV^dx''dt + habdx^dxK (6.17) 

If V^' = then this metric has a Euclidean signature, whereas if V^* 7^ then the 
metric is complex and its signature is not so easily defined. There is a sense however 
in which it is still Euclidean. At any point Xq one can make a complex coordinate 
transformation = — it V°'\^^ (or equivalently add a complex constant to the 
shift), to obtain the metric 

ds^ , = N'^dt'^ + hijdx'dx^, (6.18) 

at Xq. Thus the signature is Euclidean at any point modulo a complex coordinate 
transformation. Following the Lagrange multiplier prescription, the electromag- 
netic field is made complex by rotating the Coulomb potential $ — > z$ which 
changes the Maxwell field tensor as 

Fta iFta, Fat iFat, and Fab Fab, (6.19) 



CHAPTER 6. QUANTUM CREATION OF BLACK HOLE PAIRS 



129 



where as usual the Latin indices indicate a restriction to the spatial slices. If the 
original Lorentzian metric and electromagnetic field were solutions to the Einstein- 
Maxwell equations, then so are this complex metric and electromagnetic field. 

I now show that this complex solution can be matched onto the real solution 
from which it was derived. 



6.3.2 Matching the complex to the Lorentzian 

The obvious hypersurface along which to match the Lorentzian solution to its com- 
plex "Euclidean" counterpart described above, is a hypersurface of constant t. I 
specialize the general metric ( |6.16D to the stationary, axisymmetric case where 



X 



: 0, = 9, and = r. Then, V = [V'f'{r,9),0,0], N = N{r,9), and 
hab = diag[h^if,{r,9),hQg{r,9),hrr{r,9)]. This specialization will remain general 
enough to cover the cases of interest in this thesis. 

Now, consider how the complexified solution does or does not match onto the 
Lorentzian solutions across a surface of constant t. First, the unit normal to is 
Ua = ±N[dt]a where [dt]a is the coordinate version of dt. Choosing it to be forward 
pointing on the Lorentzian side and consistently oriented on the "Euclidean" side 
Ua = —N[dt]a in each case. Then, on the Lorentzian side the induced metric is 
hab = dab + UaUb whilc ou the "Euclidean" side it is hat = gab — UaUb which are both 
equal to diag[h^if,, h^g, hrr]- Thus, the induced hypersurface metrics match and so 
a geometrical matching is possible. In the same way, the same vector potential 
Aa is induced from both sides, so from a purely Hamiltonian point of view, the 
configuration variables match. 

Of course for the matching to be smooth, both sides should also induce the 
same extrinsic curvature on the surface (as discussed by Israel in and already 
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discussed for timelike surfaces in section |3.6| of this thesis). Unfortunately with 
Ua as defined above, the extrinsic curvatures Kab = e^e^VaUp are not the same0. 
Namely on the Lorentzian side, 



Kab = e^efM„;/3 







2N 



2N 



2N 






27V 






(6.20) 



while on the "Euclidean" side the extrinsic curvature is iKab- In a similar way, the 
induced electric field on the "Euclidean" side is iEa where Ea is the Lorentzian 
field. 

Then, from the Hamiltonian perspective adopted in this thesis the situation 
is as follows. Configuration variables hab and Aa remain real under the complex 
transformation, while their conjugate momenta 

become purely imaginary along with the Lagrange multipliers. 
Thus a matching is possible, though it isn't smooth. 



The conclusions of section |3.6| for spacetimes where there is an extrinsic curva- 
ture discontinuity across a timelike hypersurface apply equally well in this section 
where the discontinuity is across a spacelike surface. That is, the discontinuity 
corresponds to a thin shell of matter. In this case the stress-energy tensor repre- 
senting the matter is imaginary and since it is spacelike exists only instantaneously. 
This is unusual to say the very least, but then again the surface separates re- 
gions with different metric signature so perhaps it isn't surprising that something 
strange might occur at that surface. What is more of a concern however is that 

the presence of this strange matter at the borders of the instanton might shift the 
^In [H ^ Robert Mann and I took a slightly different view of this by letting Ua become 
imaginary over the instanton. Then the same extrinsic curvatures are induced on the surface. 
Here though I choose not to take this view. 
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action of the solution away from extremality. If that is the case then the instanton 
cannot be used to approximate the full path integral. This matter deserves further 
investigation, though I will not do that hereQ. 



In section O, where I select an appropriate action for evaluating the path 
integral, it will be seen that this situation of real and complex fields actually 
integrates quite nicely into the path integral formalism, but at first glance the 
discontinuities are a bit disturbing. It is clear however, that the complications 
have arisen from the inclusion of rotation. In earlier pair creation studies (such 



as 1^ ^ |7^, [Tj) there was no rotation which meant Kab = 

and the geometric matching was smooth. The discontinuity in the electric field 
remained, though it wasn't usually considered. 

Before moving on, I'll point out that by the traditional methods of instanton 



construction such as those used in ||73|, |6^, the situation would be even worse. 
The standard method would require that I analytically continue as many parameters 
of the metric as necessary to arrive at a real and Euclidean solution to the Einstein- 
Maxwell equations. For example, with the Kerr-Newman-deSitter solutions, which 
will soon be under consideration, the rotation and electric charge parameters would 
be made complex (a ia, Eq iE^) so that the Euclidean metric and electric 
field would be real. Although this approach avoids dealing with complex metrics, 
it incurs several serious problems of its own. Specifically, sending a —>■ ia and 
Eq iEq means that the surface metric hab itself is affected by the transformation. 
In detail, the polynomial Q, and functions Q and Ti (defined in and before equation 
(|6^ )) are all changed and so hrr ^ hrr, hoe ^ h^e, and h^^ ^ h^^. That this is 
not just a problem of coordinates is made clear most dramatically by the fact that 



^From the point of view adopted in |j, as discussed in the previous footnote, this problem 
doesn't arise because the extrinsic curvatures match exactly. 
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the change in Q ^ -f + (1 + ^a'^)r^ - 2Mr - Eq + Gq will certainly shift, 
and often entirely change its number of roots, which means that the "horizon" 
structure of the spatial surface will be different on each side of S^. Thus, by the 
traditional method the "Euclidean" solution emphatically would not match onto 
the Lorentzian solution. 

Given that the matching conditions are the only existing prescription that defini- 
tively link instantons with physical Lorentzian solutions, I choose to keep what 
matching conditions I can, abandon the requirement that the full spacetime metric 
be real, and proceed with the calculation. 

6.3.3 Putting the parts together 

With these general steps taken, I'm now ready to finish off the instantons. They will 
come in three classes: i) those creating spacetimes with two non-degenerate horizons 
bounding the primary Lorentzian sector (this case will create Nariai and lukewarm 
spacetimes), ii) those creating spacetimes with only a single non-degenerate horizon 
bounding the Lorentzian sector, (this case will create cold spacetimes and ultracold I 
spacetimes), and iii) and those creating zero horizon spacetimes (here, the ultracold 
II spacetime). 

Spacetimes with two nondegenerate horizons 

By the procedure described above I have found a complex solution that may be 
joined to the Lorentzian solution from which it was generated. However a subtlety 
arises in that the constant t spatial hypersurfaces of the nondegenerate KNdS and 
Nariai spacetimes both consist of two Lorentzian regions that are connected to 
each other across their corresponding horizons, while the constant t hypersurfaces 
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Figure 6.5: Construction of a two-horizon instanton. The radial/time sector is 
shown. The heavily dashed hues indicate that the solution continues in that direc- 
tion. 

of the complex solution consist of only one such region. The complex solution may 
be connected to both sections simultaneously by the following procedure (that is 
illustrated in figure |6.5| ). 

First, connect half of a full Lorentzian solution (the region bounded by the 
outer black hole and cosmological horizons) to each of the t = and i = ^ 
hypersurfaces of the "Euclidean" solution (as in figure |6.5| a). Next (figure |6.5| b) 



fold the construction over, and identify outer horizon to outer horizon, and inner 
horizon to inner horizon (figure |6.5| c). The t = constant hypersurfaces of the 
Lorentzian part of the construction now consist of two regions with opposite spin 
and charge, and are the complete t = constant hypersurfaces of the maximally 
extended but periodically identified KNdS solutions that I considered earher. 
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Next note that the metric at any point of the Riemannian part of the construc- 
tion is 

ds^ = N^dt^ + hijdx'dx^ (6.21) 

under the coordinate transformation that ehminates the shift at that point. At the 
horizons A^^ for these solutions. Therefore it is reasonable to identify the entire 
time coordinate along the horizons as a single time (figure |6.5| d). The instanton is 
nearly complete. The "Euclidean" part is smooth everywhere except at the points 
where I made the identification and probably introduced conical singularities. 

Now, for a given horizon at r = rh, I can find a period Pq such that lim^^rh '^"f''^ = 

y flrr 

2tt which in turn implies that the conical singularity has been eliminated. This is 
the same condition used in calculating the temperature of the horizons in section 
'6.2.3 , and so those results may be reused here. Hence the only double-horizon 
cases where the conical singularities at the two horizons may be simultaneously 
eliminated (figure |6.5| e) and so the only cases where the instanton will everywhere 
be a solution to the Einstein equations, are the lukewarm and Nariai instantons, 
for which 

= and Pr = T774^ (6-22) 



respectively. Q' = and r^h is the radius of the outer black hole horizon in the 
lukewarm solution. Then, the full construction of Lorentzian and Euclidean parts 
is smooth everywhere, except on the transition surface where there will be a 
mild jump discontinuity in the extrinsic curvatures. 

Next consider the single-horizon spacetimes. 



CHAPTER 6. QUANTUM CREATION OF BLACK HOLE PAIRS 



135 



Spacetimes with one non-degenerate horizon 

It is now fairly easy to build the single non-degenerate horizon instantons for the 
cold and ultracold I spacetimes (even though the cold spacetime has two horizons, 
the inner horizon is a degenerate, double horizon). For these spacetimes, attach 
half-copies of the Lorentzian spacetime at the t = and ^ = -y hypersurfaces of the 
complex Riemannian section (figure |6.6| a) . Then fold and identify the cosmological 
horizons to reconstruct the full Lorentzian t = constant hypersurfaces (figure |6.6| b 
and c). Next, identify the time coordinate along the cosmological horizon (figure 



5l6|d). Finally, with just one horizon choose 



pcoM ^ _^![^^^h+^ and P^^'' = 2n, (6.23) 

where Q'irch) = and Vch is the radius of the cosmological horizon. Then 

the instanton will have no conical singularities (figure |6 



No-horizon spacetimes 

This time the construction is less definite. With no identifications being made, 
and no horizons to define a period, the instanton has indefinite period creating two 
disjoint spacetimes (figure This corresponds to the ultracold II case. 



6.4 Choosing an appropriate action 

As was discussed in chapters ^ and^, if one chooses a Lagrangian action J, takes its 
first variation SI over a finite region M, and solves 61 = 0, then the solution includes 
not only field equations in the bulk, but also boundary conditions on the fields over 
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Figure 6.6: Construction of a one-horizon instanton. The radial/time sector is 
shown. The heavily dashed hnes indicate that the solution continues in that direc- 
tion. 
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Figure 6.7: Construction of a no-horizon instanton. The radial/time sector is 
shown. The heavily dashed lines indicate that the solution continues in that direc- 
tion. 

dM. Thus in choosing an action that is appropriate to a particular situation, one 
must keep in mind the implied boundary conditions that are attendant upon it. 

In particular, in the path integral formulation of gravity for finite regions of 
spacctime, the action choice also fixes boundary conditions that the possible paths 
must satisfy, and therefore restricts the allowed parameter space of those paths. 
There are two ways of approaching the choice of how the parameter space should 
be restricted. The first is from a geometrical/topological point of view. There 
one considers what geometrical properties the paths should have so that they will 
properly match onto the Lorentzian solutions. The second way is much more phys- 
ical and considers what physical restrictions should be placed on the paths so that 
they will produce the types of spacetime that one is interested in. That is, one 
demands that the created spacetime has certain physical characteristics such as 
horizons, temperatures of those horizons, and a particular angular momentum or 
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electric/magnetic charge and then enforces those same restrictions on the "paths" 
so that they will create the correct spacetimes. Happily, as will be seen below, these 
apparently disparate approaches complement each other and produce compatible 
lists of restrictions. 

First from a geometrical point of view, it is essential that the "paths" match, 
onto the Lorentzian solution along the interface surface S2. That is, they should all 
induce the correct surface metric hat and vector potential Aa on S2. Examining the 
Hamilton- Jacobi variation ( p. 201) (the orthogonal version is sufficient in this case) 
it is clear that the standard action functional has this property. Note however that 
the formalism does not guarantee that the conjugate momenta P"^ (or equivalently 
the extrinsic curvatures) and S"" will match as well. In an ideal world both would 
be fixed but since they are conjugate to each other this is not possible. Given this 
and the fact that the instanton work showed that for that solution the conjugate 
momenta don't in fact match across the transition surface, I'll fix the configuration 
variables and leave the other two free. 

Continuing with the geometry recall the conditions that were placed on the in- 
stantons. Namely I required that they have only one boundary (S2 that matches 
onto the Lorentzian solutions) and further that they be smooth and without conical 
singularities. That is I demanded that = (because the foliation of the space- 
time is othogonal to the boundary in this case, I'll drop the bar notation) at the 
coordinates values of r corresponding to non-degenerate horizons in the Lorentzian 
solution and further that 



where is the coordinate of the horizon. At first glance that second condition 
appears to be awkward and abstruse but in fact it is quite straightforward to show 



lim 



lim 



Pr^drN 



(6.24) 
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that if = at r/i then 

2 

hm Np = -n^daN, (6.25) 

where p is the pressure defined at the end of section p.4| . Now n"" = -^d^ and 
so to avoid conical singularities one must fix Np. To ensure this, it is more than 
sufficient to fix Ns""^. Turning again to the action variation ( p.20| ) the reader will 
note that N is already fixed while Ns""^ has been left free. Adding 

Alp = -I dt (fx^aNp, (6.26) 

to the action the situation is corrected and one can force the paths to be closed 
and smooth at the points corresponding to non-degenerate horizons. 

There is also a nice physical interpretation of these conditions. Namely fixing 
the lapse to be zero at Vh means that there will be an (apparent) horizon at that 
point while putting the restrictions on Np fixes the temperature of those horizons 
(see section |6.2.3D . Therefore enforcing these conditions at means that there will 
be a horizon of predetermined temperature there. If there are two non-degenerate 
horizons then each will have a temperature. By the no-conical-singularity require- 
ment of geometrical smoothness, they must have the same temperature and so 
geometrical smoothness is equivalent to the thermal equilibrium of the final state. 

Next consider what should be done at a degenerate horizon such as that found 
in a cold spacetime and what restrictions should be placed on the "paths" that 
might create it. To match onto the Lorentzian solution all paths must have the 
"tapered horn" shape characterized by A^ — > at the degenerate horizon. Since 
the horizon is an infinite proper distance from the rest of the spacetime, there is no 
need to worry about conical singularities, and therefore no need to fix the pressure. 
Instead leave aab fixed to ensure that the metric will have the correct asymptotic 
behaviour. Thus at degenerate horizons do not add Alp to the action. Note that 
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this geometrical behaviour is the basis of the claim |^ that an extreme black hole 



doesn't have a fixed temperature but instead can be in thermodynamic equilibrium 
with any background. 

Having set the boundary conditions to ensure that the spacetimes contain black 
holes it is natural to fix the angular momentum and electromagnetic charge of 
those black holes. Afterall the ultimate intention is to calculate the pair creation 
rates for pairs of black holes of specified mass, angular momentum, and charge so 
these quantities must be fixed in advance or else the path integral will calculate 
the creation rate for some other situation. At first it might seem natural to fix e as 
well so that one could specify the mass of the holes being created but as was noted 
above one can fix or £ but not both. A^ must be fixed so that the black holes 
can be guaranteed to exist, so e has to be left free. That said, to fix the angular 
momentum one must add 

AIj = - I dt I d^xy/^V^ja, (6.27) 

to the action. 

Next consider fixing the electromagnetic charges. First recall from chapter ^, 
that by choosing to work with 1^, F^fs, and A^ I have automatically excluded 
magnetically charged solutions from consideration. At the same time however. 



the electric charge has been left free (see the variation ( |4.26|) ). By the previous 
paragraph it should be fixed and so add 



A/ei = - [ dt [ d^xN^n^F'^^Ap, (6.28) 



to the action. There is a choice of whether to fix S"' or Aa on the boundary S2. 
As has already been noted Aa is the appropriate quantity to fix and that doesn't 
require an extra boundary term on E2. 
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By contrast if I want to consider magnetic black holes I use -kFap, and A* . 
Then electric charges are automatically eliminated from consideration and 

^Img = - I dt I d^xN^n^ (6.29) 

J JQt 

should be added to the action to fix the created magnetic charges. No additional 
boundary term is required to fix A* on S2. 

Switching to the thermodynamic interpretation of the path integrals it is imme- 
diate that what is being considered here is a canonical partition function. That is, 
extensive variables (angular momentum and electric/magnetic charge) are fixed ex- 
cept for the energy which is left free in favour of holding the temperature constant. 
This is then in accord with the standard approach to pair creation calculations 
which uses that partition function This choice then ensures that created 

spacetimes are in thermal equilibrium, that there is no discontinuity in physical 
properties such as electromagnetic charge and angular momenta at the juncture of 
the paths and the Lorentzian solution, and from a geometric point of view that the 
paths are smooth and match onto the Lorentzian solution. 

6.5 Evaluating the actions - pair creation rates 
and entropy 

As noted earlier, creation rates for these spacetimes are proportional to e"^^™"*, 
where Itnst is the numerical value of the action of the appropriate instanton. Now, 
as was laid out in section |6.1|, those rates are calculated only up to a normalization 



factor. Evaluating this normalization factor would involve fully evaluating another 
ill-defined path integral so I will side-step that issue by calculating the probability 
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of creation of these spacetimes relative to deSitter space. Then the normahzation 
factors cancel each other out and the relative probability of creation of the black 
holes in a deSitter background is 

P = exp(2/,5-2/), (6.30) 

where I is the action of the instanton, and Ids is the action of an instanton mediating 
the creation of deSitter space with the same cosmo logical constant. Conventionally, 
this probability is also interpreted as the probability that deSitter space will tunnel 
into a given black hole spacetime |T^. 

A further link to thermodynamics is found by the following argument. The 
spacelike hypersurfaces of the spacetimes that I have considered are all topologically 
closed and with finite volume. Then, the energy is trapped in the hypersurfaces and 
so they can be interpreted as having constant energy even though that condition has 
not been enforced by boundary conditions [[71] , |52| . By that reasoning the canonical 



partition function is equivalent to the microcanonical partition function in this case, 
and so as is standard for a microcanonical partition function, the entropy of the 
created spacetime is 

5 = In ^2 = -2/. (6.31) 

Thus, at least in the case where the created spacetime is quasi-static, there is a 
close connection between pair creation rates and the entropy of the spacetimes and 
in particular it is consistent with the idea that the entropy is the logarithm of the 
number of quantum states. With all of this in mind 1 evaluate the appropriate 
action for each spacetime. 

Momentarily leaving aside the matter terms, the appropriate action by all of 
the above considerations is 
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= ( terms that vanish for stationary solutions ) 

+Eah [ d^x^Ne - hlsH [ d'^x^Np, (6.32) 
Jb^ 2 Jb^ 

where the subscript SH indicates a sum over all single, non-degenerate horizons 
and AH means a sum over all horizons regardless of their degeneracy. Keeping in 
mind that = on all of the horizons, it is clear that the Ne terms are zero. 
Further, recall equation ( |6.25|) which says Np = {2/ K)n°'daN and equation ( |6.24|) 
which implies that n°'daN = (27r)/Po on a non-degenerate horizon. Then 

Ah 

hN,P,j) = -^SH-^, (6.33) 

where Ah is the surface area of the event horizons in the spatial surface E2. 

Next, consider the matter terms. For electric solutions it is a trivial use of 
Stokes's theorem to show that, 



L f d^x,/^{F^pF''^) + Ahi 
2fi; Jm 

= [ d^x^iApVaF""^)-- f d^x^/hE^Aa. (6.34) 

J M Jt, 

Of course the first term includes the constraint equation ( p. 121 ) and so is zero for 
solutions to the Maxwell equations. Thus all that is left is the second term, but it 
too is zero for the solutions in which I'm interested, and so the total electric term 
is also zero. Thus, the value of the action that keeps A^, p (if appropriate), ja, and 
Eq fixed (and Gq = 0) is 

I{N,p,j,Eo) = -^SH-^- (6.35) 

The same line of reasoning shows that the value of the action that keeps A^, p (if 
appropriate), ja, Gq fixed (and Eq = 0) is 

Ah 

^{N,p,j,Go) = — • (6.36) 
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Assuming that the spacetimes are at least quasi-static, then equation ( |6.31| ) 
says that the entropy of these spacetimes is equal to one-quarter of the sum of the 
areas of non-degenerate horizons bounding the Lorentzian region of the spacetime. 
Consistent with references ||9] and |7l], the degenerate horizon in the cold case 
does not contribute to the entropy of the cold spacetime. 

Using these general formulae for the pair creation rates and entropy of the 
spacetimes, I now consider each of the specific spacetimes separately. 

The lukewarm action: In this case, there are non-degenerate cosmological and 
outer black hole horizons. Therefore the numerical value of the action of the elec- 
tric/magnetic instantons is 



(6.37) 



8 

where Ac and Ah are respectively the areas of the cosmological and outer black 
hole horizons at and in the Lorentzian solution. 

The Nariai action: Again there are two non-degenerate horizons, this time at 
p = ±1. Therefore the total action of the electric/magnetic Nariai instantons is 

In = ^ , (b.38j 

where Ap=±i is the area of the horizon at p = ±1. Note that for the Nariai solutions 
Ap=i Ap=—i. 

The cold action: Here there is only one non-degenerate horizon, and so 



Ac 7r(rg + a2) 
Ic = -— = ^2 ' (6-39) 



where Ac is again the area of the cosmological horizon at Vc- 
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The ultracold I actions: Again there is only a single nondegenerate horizon, this 
time at i? = 0. The action of the magnetic ultracold I instanton is then 

Aq 7r(e2 + a^) 



lUCI 



(6.40) 



8 2x' 
where is the of the Rindler horizon as i? = 0. 

Ultracold II Actions: There are no horizons whatsoever for this case, and so 

lucii = 0, (6.41) 
irrespective of the chosen period Pq of the time coordinate, which is good since as 



was seen in section p.3|, that period is not specified by the formalism! 



In figure |6.8| , I plot these actions as a fraction of the action of the instanton cre- 
ating deSitter space with the same cosmological constant. The instantons/created 
spacetimes are parameterized by and ^M'^. For all cases I, Ids < and from 
the diagram it is clear that |/| < \Ids\- Then I^s — I < and so each of the space- 
times considered above is less likely to be created than pure deSitter space with the 
same cosmological constant. Note that the Nariai spacetime is the most likely to 
be created provided the parameter values are such that the instanton exists, while 
the cold spacetime is the least likely to be created. As might expected on physical 
grounds, smaller and more slowly rotating holes are more likely to be created than 
larger and more quickly rotating ones. As — and M 0, the creation rates 
approach those of deSitter space which again is physically reasonable. 



6.6 Reflections on the calculation 



The approach to pair creation taken here is a bit different from that taken in most 
of the literature and because rotation has been included new issues have arisen that 
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were not present in those papers. Thus in this section I will compare the methods 
and examine those issues a little further. 

First I will compare the way I calculated actions here with the way it was done 
in reference |j7l]| (which is representative of the more traditional calculations done 
for non-rotating black holes). There, the fact that the instantons are closed and 
smooth at the points corresponding to the non-degenerate horizons was taken to 
mean that no boundary terms need be considered there, implying that the basic 
action used for the lukewarm and Nariai instanton should be 



hid = [ d^x^ {n-2A-F^)-- I d'xVhK, 



(6.42) 



which from my point of view is the Lagrangian action ( ^.1| ) with the boundary term 

- ! d^Xyf^Q (6.43) 



Jb 

added on. It is easy to see that this term is equivalent to the pressure term 

'^^.t], (6.44) 

evaluated on the equivalent horizons. To see this, note that 6 = — k = 

— — i=9rlncr, and ^= at each horizon. So on those horizons G = — f and in 



[ d^xN^- = - [ d^xN^ 
JB 2 



the absence of rotation my approach is equivalent to that of . 

For the cold case, k still vanishes on the boundary and so the inclusion of the 
B term in is equivalent to the omission of the pressure term in my calcu- 
lation. Finally, in the ultracold cases = everywhere and so once more the 
omissions/inclusions are equivalent. 

For electric instantons in both calculations, electromagnetic boundary terms are 
added to the action to fix the electric charge for all paths considered in the path 
integral. Further, in both calculations for solutions to the Maxwell equations, these 
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boundary terms may be converted into the F"^ bulk term that was used in this 
work. For electric instantons that earlier work added a boundary term 

- / d^x^E^Aa (6.45) 

on the S2 boundary (the only boundary for those instantons). For solutions to the 
Maxwell equation this is then equivalent to adding a bulk term 

- f dS^gF^pF^" (6.46) 

to the action which then makes it numerically equivalent to the action + AI^i- 
For the magnetic case it was argued that nothing needed to be added since the 
magnetic charge was already fixed on the boundary. Note however that while that 
approach works out numerically it not quite right because, as was noted in chapter 
^, if one assumes that a single Aa covers M then no magnetic charge can exist. 

Although for non-rotating instantons the approach here is equivalent to earlier 
ones, differences arise when rotation is included. In earlier approaches [^, ^ 



39i [7T| , p^ , |53| , p2| , pTj P^ , |69| , |70| there was no provision made for fixing the 
angular momentum and so the action differs by the term AIj and its omission is 
tantamount to working with an incorrect thermodynamic ensemble. Evaluating the 
action of rotating instantons with (|6.42|) will not yield the preceding relationships 
linking surface areas, actions, and entropies. Indeed, using (|6.42|) the creation rate 
of rotating black holes is enhanced relative to that of non-rotating black holes and 
with an appropriate choice of physical parameters may be made arbitrarily large. 

Second, around the same time that this work was originally published, Wu 
Zhong Chao published a series of papers on the creation of a single black hole (see for 
example pGl E%) using a slightly different set of instantons to create spacetimes that 
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are not in any kind of equilibrium (thermodynamic or otherwise )0. He recognized 
that the angular momentum needed to be fixed but used an ad hoc approach to work 
out what the angular momentum fixing term should be. For the cases considered 
however, that term was equivalent to the one used here. 

However, despite the results being similar, his approach was quite different. In 
the first place he asserted that his approach could create a single black hole. From 
a physical point of view, this would violate conservation of angular momentum and 
electric/magnetic charge. Even apart from this, the instantons that he considered 
do not properly match to real Lorentzian solutions for two reasons. In the first 
place there are no periodic identifications of the universal covering space of the basic 
KNdS solution that can be made such that hypersurfaces of constant t will contain 
only a single black hole. The smallest number of black holes that may be contained 
are the two discussed here. Second, as argued earlier, an analytic continuation of a 
to ia and Eq to iE^ will mean in general that an instanton generated from a classical 
solution will not properly match onto that classical solution. In general, there will 
not even be the correct number of horizons available in the instanton to match onto 
the Lorentzian solution. In later papers (for example |55|, he considered the 



creation of pairs of black holes instead, but the other differences remain. 



^This type of instantons have also been proposed for use in cosmology by Hawking and Turok 
in 0. 
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Figure 6.8: The actions for the charged and rotating lukewarm, cold, and Nar- 
iai instantons plotted as a fraction of the action of deSitter space with the same 
cosmological constant. The Nariai instantons arc the meshed sheet, the lukewarm 
instantons are the lighter grey sheet, and the cold instantons are the darker grey 
sheet. The ultracold I instanton actions may be found at the "bottom" end of the 
cold sheet, while the ultracold II instanton actions are zero. 



Chapter 7 



Discussion 



In this thesis I have taken the quasilocal energy formahsm of Brown and York 
and generaUzcd it in several different directions. First, for the finite region of 
spacetime M, I dropped the requirement that the fohation hypersurfaces and 
spacehke boundaries Ei and E2 be orthogonal to the timelike boundaries B. Prom 
a theoretical perspective this obviously allows one to consider much more general 
regions M and further does not restrict the allowed variations of the metric when 
the variational principle is applied. The ensuing calculations then make it clear 
that the numerical value of the quasilocal Hamiltonian (and thus the quasilocal 
energies derived from it) is a function of the foliation of B and time evolution vector 
field T". It does not care about the foliation of the bulk M as a whole. As was 
repeatedly emphasized throughout the thesis, this is a very desirable characteristic 
for a quasilocal energy to have since the correspondence between foliations of the 
bulk and foliations of the boundary is many-to-one. Further from a practical, 
computational point of view, focusing only on the foliation of the boundary (as 
opposed to the bulk) makes it much easier to calculate the quasilocal energies seen 
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by moving observers. 

Second, I shifted the calculations from the usual Lagrangian framework into a 
pure Hamiltonian form. Of course the two are equivalent but this thesis was the 
first place where that was shown explicitly for a quasilocal region of spacetime. A 
side benefit of this shift of emphasis was that the variational calculations could 

easily be adapted to calculate rates of change of the quasilocal quantities and so 
give a slightly different outlook on conserved quantities than that discussed in the 
original Brown- York work. 

Third, I examined the reference terms that set the zero of the quasilocal energy. 
In the process of reviewing some of the extant proposals I showed that all had prob- 
lems dealing with moving observers in Minkowski space - namely such observers 
measure non-zero energies in flat space. To deal with this problem I proposed a 
new definition that embeds the two-surface of observers into a four-dimensional 
reference space. While this new reference term is by no means perfect at least 
it ensures that the action and quasilocal energy of flat space is zero in all cases. 
Unfortunately, the inclusion of this reference term (or any other reference term for 
that matter) complicates the Lorentz-like transformation laws derived for the ref- 
erenceless quasilocal energy. Specifically, the reference terms must be transformed 
with respect to a different boost velocity than the referenceless terms. 

Fourth, with the reference term discussed above, I showed that it is possible to 
recast the (generahzed) Brown- York QLE in an operational form. Roughly speaking 
I showed that the QLE contained by a closed two-surface is exactly equal to the 
total stress-energy of a particular thin shell of matter. That thin shell must be 
embedded in the reference space such that: i) it has the same intrinsic geometry as 
fit, and ii) outside of that surface the spacetime geometry is identical to that found 
outside of fit in the original spacetime. 
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Finally, I added Maxwell and dilaton matter fields into the mix. These have 
previously been considered for orthogonal foliations, but my work was the first 
to examine them in the non-orthogonal case. Their integration into the purely 
gravitational scheme of things proceeded smoothly resulting in small, but usually 
non-qualitative, changes. The only significant problem arose because the quasilocal 
formalism as constituted doesn't allow magnetically charged configurations of the 
Maxwell field. They can be accommodated by switching to a dual formalism but 
in doing so electric charges were excluded. Thus, the formalism was shown to allow 
either electric or magnetic charges but not both. 

Having developed this formalism, I applied it to a variety of situations, both 
classical and quantum, to investigate whether or not the QLE can reasonably be 
thought of as defining a physical energy. To get some orientation, I started by 
examining Schwarzschild and Reissner-Nordstrom spacetimes. For a spherical set 
of observers far from an RN source, I showed that the geometric QLE matches a 
Newtonian intuition of how energy should be distributed in the spacetime including 
the contributions from electromagnetic forces. I then showed that the total QLE 
included a contribution to the energy from the base Coulomb potential for the 
spacetime (by which I mean the Coulomb potential that would exist even if there 
was no charged matter in the spacetime). However, the KiUing- vector-adapted QLE 
didn't behave in such an intuitive way. 

From there I considered the quasilocal energy measured by boosted observers. I 
started with the easy case of radial boosts and showed that while the referenceless 
QLE increases in the expected Lorentzian way with the boost, the numerical value 
of the referenced QLE actually decreases! In particular 1 showed that for observers 
far from the hole, the boosted and referenced QLE equals M/7 where 7 is the 
usual Lorentz factor. Later on in the naked black hole section, I showed that this 
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is due to the competing relativistic effects of the motion of the observers and the 
gravitational field. I then examined the trickier case of z-boosted observers. The 
results were much more complicated, but far from the hole I found that the boosted 
and referenced Hamiltonian again measures M/j. 

Next, I examined the quasilocal energies measured by spherical observers who 
are either hovering around or falling into a class of naked black holes. Such observers 
respectively feel either negligible or Planck scale transverse tidal forces. In contrast, 
I demonstrated that the static observers measure a large geometric QLE while the 
infalling set measure a very small geometric QLE. This can be explained because 
the extremely strong tidal force corresponds to a massive Lorentz boost of the 
reference terms which in turn means that the relativistic effects of the motion 
completely overwhelm those of the gravity. Thus, even though the unreferenced 
QLE and reference terms are both hugely boosted, at the same time they converge 
towards the same value so the difference between them goes to zero. 

As a final classical example, 1 applied the formalism to investigate energy fiows 
that arise during gravitational tidal heating. I successfully used it to reproduce the 
standard Newtonian and pseudo-tensor result and explain their gauge ambiguities 
in terms of fiuctuations of the quasilocal two-surface. Thus I demonstrated the 
utihty of the formahsm in an astrophysical context which also helps to boost its 
claims to physical relevance. 

The thesis finishes up in the last, rather long chapter, by applying the quasilocal 
formalism to study pair production of rotating black holes in deSitter space. It was 
seen that the results for non-rotating black holes can be qualitatively extended to 
the rotating case. That is, created spacetimes can be classified as lukewarm (regu- 
lar KNdS solutions where the horizons are in thermal equilibrium), cold (extreme 
KNdS solutions), and Nariai (a hmiting case where the outer black hole horizon 
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approaches the cosmological horizon). The entropy of such spacetimes continues 
to be proportional to the surface areas of the non-degenerate horizons and the pair 
creation creation rates continue to be proportional to the negative exponential of 
those entropies and suppressed relative to the creation of a pure deSitter space. 

To obtain these results I was forced to make a choice between the real instan- 
tons that are usually employed to evaluate the path integrals and the standard 
Lorentzian solution/ "Euclidean" instanton matching conditions. Since the match- 
ing conditions are the only way that I know of to associate an instanton with a 
given Lorentzian solution, I opted for the matching conditions and allowed com- 
plex instantons. 

Using the quasilocal formalism to fix the ensemble of paths considered in the 
path integral, I showed that the standard Einstein-Hilbert action is not the ap- 
propriate action to use for rotating pair creation. In particular it does not fix the 
angular momentum of the ensemble and therefore docs not guarantee the creation 
of a black hole pair with a prespecified angular momentum. A careful application 
of the quasilocal formalism allowed me to identify the correct action and so obtain 
physically reasonable results. 

Possible future work 

I see two main directions in which to continue work started in this thesis. First, I 
have concentrated almost entirely on the Brown- York definition of QLE. However, 
there are many other Hamiltonian based QLE's and it would be of interest to 
examine them closely in the same way that I have dealt with the Brown- York QLE 
here. For example, it would be interesting to examine how the various definitions 
measure the energy flow in the tidal heating example. More generally one could 
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compare them in the hmit of weak gravitational fields where they could also be 
compared with the perturbative treatment of gravity as a spin-2 field propagating 
in a flat background. There is a gauge ambiguity in such a treatment and it seems 
possible that each measure of QLE might correspond to a different gauge choice. 
Even more generally, a close examination and comparison of their mathematical 
formalisms might help to shed light on each. 

A second project can be found in the pair creation calculation. There I noted 
that there was a difference between thermal and thermodynamic equilibrium in 
two horizon spacetimes. The first is defined by the temperatures of the horizons, 
which in turn are most easily calculated using Euclidean quantum gravity tech- 
niques. Thus such a spacetime is in thermal equilibrium if a regular instanton can 
be constructed from it. At the same time such spacetimes do not have to be in 
full thermodynamic equilibrium since angular momentum and/or electromagnetic 
charge could still be exchanged between the horizons. This dichotomy deserves a 
fuller investigation. To this end it would be profitable to investigate the evolution 
of black holes in deSitter space using the techniques of quantum field theory in 
curved spacetime to calculate rates of particle emission, and the mass, charge, and 
angular momentum carried off by those particles. Much work has been done in 
this area for asymptotically fiat spacetimes, but cosmological spacetimes are not 
as well studied. In particular no one has studied them when they are in thermal 
equilibrium. Apart from understanding the difference between the notions of equi- 
librium, this issue is quite topical with the recent interest in non-zero cosmological 
constant spacetimes, that has arisen from the astronomical measurements which 
indicate that our universe may have a positive cosmological constant and the string 
theory inspired AdS/CFT correspondence. 



Appendix A 



Hamiltonian Calculations 



This appendix presents the calculations behind the results of chapters |^ and 



A.l Foliating the gravitational action 



First I decompose the action (13^1 



/-/ = — / rf^xv^(7^ - 2A) + - / d^x^K-- [ d^xy/^Q 
2k Jm Jt. Jb 

+ — / d'^x\^ smh~^{ri), 
^ Jn 

as discussed in chapter ^ into three-fields and time derivatives of those fields defined 
on the foliation hypersurfaces and Qf To start, with the help of the Gauss- 
Codacci relations one can rewrite 

{TZ-2A) = R-2A- + K^pK'^^ - 2Va{Ku" + a°), (A.l) 

where = u^V = -^DaN is the acceleration of the foliation's unit normal vec- 
tor field along its length. Then using Stokes's theorem to move the total derivative 
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out to the boundary, it is trivial to show that 

[ c/^xv^(7^ - 2A) = [ d^x^{R-2A-K^ + K^pK^P) (A.2) 

-2 / d^xVhK-2 I £xy/^ [Kt] + n^a"") . 
Jy, Jb 



Next, referring back to the expressions for and given in equations ( |2.3| ) it 
is a simple matter to show that 



A 



(A.3) 



Then, these two results can be combined to rewrite the Lagrangian as 

/-/ = ^ f (fx^{R~2A- + K^pK'^^) (A.4) 



B 



(fxy/a sinh. '^{rj). 



The next step in the process combines (the matter-free versions of) the Einstein 
constraint equations (|2.39|) and (|2.40|) with the extrinsic curvature of in Ai 
written as Kajs = —\£uhai3 = {£Thaf3 — ^D^aVp)), to rewrite the integrand 
of the remaining bulk term of the Lagrangian as 



7^ 



(A.5) 



-9 



where P"^ = ^ [Kh'^'^ - fsT"^). Recalling that = NVh (eq. (|2l0|) ) and once 
again using Stokes's theorem, this time on the hypersurfaces to move the total 
divergence term out to the boundary surfaces Qt, I can rewrite the action as 

I-L = [ dS{P"^£TKp- NH- V'Ho) (A.6) 

J M 

-- [ dt [ d^x^ {Nk - ^[K^p - KKp\n^ - NXu^VaV) 
Jn 
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Up to this point I have been working with the fohation of M and therefore with 
the lapse A^, shift V"', and normal vectors u°' and n". On the term evaluated on 
B, I now switch to work with the foliation of B and therefore the boundary lapse 
N, the boundary shift V°', and normal vectors and fi". Then, 

Nk = ^Nk - r]Na''^VaUp, (A.7) 

and 

- V {K^p - KKp) vP = iVr/a"^VaM/3 - Nv^k + n^'Vl^VpUa + XV^Vpr]. (A.8) 

Writing the timelike vector T" in terms of the boundary quantities (eq. ( p.7[ )) it is 
easy to see that 

/ d'^x^smh-\r]) - dt d'^x^XNu'^V (A.9) 
= I dt I d^x{{£T^) sinh"^(r7) + v^AV^^V^r^) . 
Thus, the action takes its final form given in eq. ( |3.3| ). That is 

I-L = [ d^x{P''f^£Thaf3- NH-VH^) (A.IO) 
+ dt d^xP^{£TV^) - dt d^x^{Ne -V^ja) , 
where e = ^k = — -^cr^^Va^/? and = ^a^u"''V fsn^, and P^ = ^ sinh""^ r]. 

A. 2 Gravitational Hamiltonian variation 

Next, I calculate the variation of the Hamiltonian with respect to the Lagrange 
multipliers and V"", surface metric hat, and its conjugate momentum P""^. Be- 
cause ^/o^ and P ^ are functions of these other quantities they are automatically 
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varied as well. From eq. ( |3.4| ), the Hamiltonian is 

Ht= f £x{Nn + V''Ha)+ [ d?x^{Ne-V''3a) (A.ll^ 



where 



Ha = -2DtP' 



ai 



N = XN, 

= V - {V''nb)n\ 

1 2 . 
e = —k H ^vP"" naTib, and 

2 1 
Ja = jfaabP^^n^ a%{sm\i-^ri). 

The calculation is quite lengthy and so is tackled in parts. 



Variation of the bulk term 

I start with the bulk term 



Ebm = [ d'xVhiNH + V'na), (A.12) 

and calculate its variation with respect to each quantity. 

First, the variation with respect to the hypersurface momentum P'*'' is easily 
calculated as 

5p..Hbik = j^d?x{^ Pah-\phab + 2D(,H)) (A.13) 



-2 / £x^naVbSP''\ 

The boundary term arises from using Stokes's theorem to remove a total divergence 
to the boundary. 
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Next and more challenging is the variation with respect to the metric hat- To 
this end note that 



+VhDa {h-'^h'' [DMbd - DJh,]) . 



(A.14) 



This may be calculated from first principles, but the easiest way to do it is to 
simply adapt the variation of the four dimensional Ricci scalar with respect to the 
four metric Qais- Such a calculation may be found in any text book that deals 
with Lagrangian formulations of general relativity (for example Wald [^). Then, 
recalling from the previous section that the acceleration vector can be written in 
terms of the lapse as at = ^DbN, a not too lengthy computation obtains 



St 



2k 



■[i?-2A] 



(A.15) 



+ / d^x^^ {-n''h!"'DMab + n^h^'DaShab + a'n'^SKb - [a'^n^h'^'Shab] , 

where (^^G"'' = R"^ - \Rh''^. The hab variation of the rest of the NH is quickly 
found to be 



8h 

I >'n. 



d^x 



2k 



pab p 
1 If 



ab 2^ 



(A.16) 



d x—= 



pacp b }ipp<^b 

2 



6k 



ab- 



Slightly more difficult is the variation of the V^'Tia term. For that I use the relation 

1 



(A.17) 



where FJ;^ is the Levi-Cevita connection for hab- Then, keeping in mind that P'"^ is 
a tensor density (a relative tensor of weight one) and so DcP^'^ = dcP^'^ + T\^P'^'^, 
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it can be shown that 



St v/i V 



(A.18) 



-2 / d3x^fy'^PV5/ia6-^[l""rijP'^^5/i,fc 



Finally, it is trivial to calculate the variation of the bulk term with respect to 
the lapse and shift. To wit, 



(A.19) 



Focusing back on the boundary terms of (|A . 1 3| , |ArT5|JA . 1 8|) one can break up the 
variation of hab in terms of the variation of (Tab and n^, where 5aab = CTa'^b^hab- 
Further, is defined as the unit covariant vector normal to the surface Qt in 
and so it is defined up to a normalization constant without reference to the metric 
(crudely, if Qt is defined to be a surface of constant r, then na\\dr). Therefore 
6na = aua and = —an'^ + P'^ where a = ^n°-n'^6hab and = —a°'^n'^6hcd- Thus 
the total variation of Huh can be written as 



6H, 



blk 



d^X [mN + Ha5V'' - [P'^^dhab + [habh^P"') 



(A.20) 



-2/ d\^{naVb + V''n,nan,)5P^' 



2k, 
2k 



[a'^Dbf^c - k^'Saab + 2ak) 



+ / d^x^ {-2V''P'''nMab - SaW^'riaUb + W^'oyMab) , 
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where 



2k 
Nk 

and 



(^[P"^Pcd - ^P^lh""^ - 4[P'=('^P/) - ^PP^^]^ 



AkN 1 

[hab\T = -^[P-b - ^Phab] + 2D^M). (A.22) 



As the notation suggests (and is discussed in section |3.2.2| ) the above equations 
define time derivatives. 

Though this expression is quite a mess, it will be substantially improved once 
the variation of the boundary terms of Ht is added on. Thus, I now calculate that 
variation. 

Variation of the boundary term 

It is simplest to calculate the total variation of 

Hhnd= [ d^xV^{Ne-V''ja), (A.23) 

with respect to all of the variables simultaneously. Then, 

5Htnd= [ d^x{[Ne-V''Ja]6^/^+^/^[66N-Ja6V'' + N5e-V''5Ja]). (A.24) 

Tackling the individual terms one at a time, first note that 

S^=l^a-b^a,b, (A.25) 

just as Sy/h = lVhh''^6hab. The 6N and terms are left as they are while the 
6e and 6ja terms become, 

NSe = '^P''''nanbS[smh-^r]]-aNe + -Na''''DaPb (A.26) 



APPENDIX A. HAMILTONIAN CALCULATIONS 



163 



+3aV''ncP''^nanb + -^Vrijian^dP'^, 



and 



-y^'^da = --^P^VVd<j"''(^<Ja6+;^V^"P'W<Ja6 (A.27) 

Again the result is a bit of a mess. Luckily, however the unpleasant terms cancel 
each other out once this is combined with Huk- 



The complete Hamiltonian 

Putting the two variations together there is significant simplification. Apart from 
cancellations, the only other computational trick required for the recombination is 
to keep in mind that is a closed surface, so 

/ d'^x^/adaV" = 0, (A.28) 
Jcit 

for any smooth vector field f]"" G Tflf. Then, the total variation of is 



where 



5Ht = / d^x{n5N + naSV''-[P''']TShab+[hab]TSP''') (A.29) 
+ I d^x^ {e6N - jadV - {N/2)s'''Saab) 
+ [ d'x^{[^TSP^-[P^]^SV^), 



= \ (k"^ -[k- n'^odja'^*) - -^r]a>dP"^ (A.30) 
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and [P is an undetermined function over Q 



a 



ab 



( N^^P^^UaUb + N^k - -dbV^ 

XVh n n 



(A.31) 



A. 3 Foliating the matter action 

This section decomposes the matter action 

I'^-l = ^ [ c/^Xv^(7^-2A-2(V^0)(V» -e~2a^F«^^"/3^ 
2fi^ Jm 

+ - f dPx\/hK I d?x^/^Q + - I (i^XA/asinh"^(?7), 



from chapter ^?1| into three-fields and their time derivatives as defined on the hy- 
persurfaces of the fohations and Qf 

First, after breaking up the purely gravitational terms as before, the bulk term 
integrand is 



P''^£tK3- NH-V^H^ 



-9 



-9 



(A.32) 



Then bringing in the Einstein constraint equations ( p.39| ) and ( p.40| ) this may be 
rewritten as 



Nk 

P''^£tKp - NTT - V^H:: + ^ ^yajj^ 



(A.33) 



2Vh 



Next, from eq. (|2.45| ) it is not hard to rewrite Sa as 



-2a</i /./T _ /T 



(A.34) 
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Using this relation and the trivial £t4> = Np+ V^Dacj), the bulk integrand ( |A.32| ) 
may be written entirely with respect to fields on the hypersurface, time derivatives 
of those fields, constraints, and a total derivative. It becomes 

P'^^ItKp + p£t<P + S^'ItA - NTT - (A.35) 

Q, = —DaS°' is the constraint equation ( |2.26| ) for the electric field with no sources. 
Thus the action can be written as shown in equation (^4.4|). That is 



j dt j d^X^^P''^£Th^f3 + p£T(p + £'"£TAaj (A.36) 

+ j dtj d''x{P^{£TV^)} 



where 



- j dt J d'xiNH^' + VH^ + T'^AM} 

- J dtj^ rf2^7^{iV(£ + £™)-F"(j„+j:')}, 

\={np£^){\^ - r^i^n") = -l^{npB^)^ and (A.37) 



^ -^{npS^)A^ = -^{n^S^)A^, (A.38) 

and "I = —AaU°', S°' = —2/ VhF^l^up, and A^ = cr'^Ap. Here, I have once again 
used Stokes's theorem and therefore the assumption that there exists a single A^ 
defined over all of M. 



A. 4 Matter Hamiltonian variation 

Next I calculate the first variation of the matter Hamiltonian with respect to hah-, 
y^, Aa, and 0, their conjugate momenta P""^, B"" ■, and p, and the Lagrange 
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multipliers A^, V"", and From eq. ( ^.6| ) the Hamiltonian is 

= I d^x[N{l-r -^Q) + V''{l-C^ + AM)] (A.39) 



where Ti^ = 0, 7i™ = 0, and Q = are the constraint equations ( p.39| ),( p^40D , and 
(|2.26| ) respectively. 



Now, the variations of the purely geometric terms were calculated in section 
A.2|, so only those of the matter terms need to be considered separately here. This 



time, it is easiest to calculate the first variations of the full expression with respect 
to each quantity separately. First for the dilaton, it is trivial to show that 

6^Hr= [ d'xm^dp), (A.40) 

J St 



where 

Nk 



p+£v(i). (A.41) 



^^^^ - 2Vh' 

It is only a little more difficult to calculate the variation with respect to dilaton 



as 



f r 2N 

6^Hr = - / d'x{[p]^ 50) + / d'xV^— [<Pl 6<p, (A.42) 



where 



[p]^ ^ '^D\ND,<P) + D,{pV^) (A.43) 

Nk 

+a^{e-'''^B'Bh - e'^^^S^) and 
Changing to the four-dimensional perspective it is easy to show that [0]- = £n(p- 
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Variations with respect to the EM terms are a httle more difficult but still not 
too bad. A few lines of calculation are required to show that 

6eHr = [ d\ {[A,]T5t) , (A.45) 



where 



[Ah]T = ^^'""^^b + £vAt,- D,[Nct>l (A.46) 
and a few more give the variation with respect to $ and Aa as 

S^^H"^ = I d'x{-[S']^6S'-NQ6^) (A.47) 



+ / d^x^^(-\e-^'''^e'''''n,BMd + V^c^''SA 



where 

[S"] ^ = -e'''''D,[Ne^'"'>Bd\ + £vS'. (A.48) 

Switching again to a four dimensional perspective, it is only a little more involved 
to show that 

= alB, = ^B, + r^e^^'^e^^.n^^^ (A.49) 

which is a generalization of one of the standard Lorentz transform laws of electro- 
dynamics. Then the term in the brackets of the third integral of (|A.47|) is equal to 

d- 

Next consider the variation with respect to P°^ and the lapse N, shift y", 



BbU^e^^SA. 



and metric hab- The variation with respect to P"'' is unchanged from the pure 
gravitational case considered in section A. 2. At the same time the variation with 
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respect to the lapse N, shift F", and metric hab is fairly easily shown to be 

5gHr = J d^x {[H"^ - ^Q]5N + [TV: + QAaW - [P'''']^^5hab) (A.50) 



where 



Nk / 1 



4v^V 2' 
Thus, the total variation of is 

= J d^x(jH"'-^Q]SN + [n'^ + AaQ]SV''-NQS^^ (A.52) 

+ / d'x {[habhSP''' - [P''']^Shab) 

d'x ([<l>]TSp-[p]TS(l>+[Aa]TS£'' - [SV^a) 

+ f d'xV^ {[e + e^'jSN - [j„ + J-]5l^'^ - (iV/2)s«^5a„b) 

+ [ d'xV^{[V^^dP^-[P^]^d^), 
Jnt 



+ 



Appendix B 



Pair creation calculations 



B.l Reducing the generalized C-metric to KNdS 

As noted in section |6.2.1| , the general Plebanski-Demianski metric contains 
conical singularities that correspond to cosmic strings or rods that supply the energy 
necessary to give black holes their extra acceleration above or below the rate of the 
rest of the universe. In this section I show that one way of eliminating the conical 
singularities corresponding to those strings/rods reduces the Plebanski-Demianski 
metric to the Kerr-Newmann-deSitter metric. This serves to emphasize that the 



global KNdS metric contains at least two black holes (see section 6.2.1 for more on 
this point). 

The generalized C-metric takes the form 



169 
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with accompanying electromagnetic field defined by the vector potential 

eoq{dT + p^da) gop{da - q^dr) 

^~^ — 2"^ ^ ^~^ — 2"^ ' y^-^l 

1 + p^q^ 1 + p^q^ 

where p, q, r, and a are coordinate functions, 

P{P) = {— - gl + 7) + 2np - + 2mp^ + (-^^ - - 7)/, (B.3) 

and = P{q) + f (1 + 9"^)- A is the cosmological constant, 7 and e are constants 
connected in a non-trivial way with rotation and acceleration, eo and qq are linear 
multiples of electric and magnetic charge, and m and n are respectively the mass 
and NUT parameters (up to a hnear factor). This solution can be analytically 
extended across the coordinate singularity at p = g, so that on the other side of 
p — q there is a mirror image of the initial solution. Thus, it can be seen as 
describing a pair of black holes on opposite sides of that p — q hypersurface. 

To apply this metric to more specific physical situations, the coordinate func- 
tions are best converted to sphcrical-type spacetime coordinates as g ^, p 
Pa + acos^ for some constants a and Pa, cr <-> and t t. Now in general, a 
periodic identification of a will introduce conical singularities at the roots of P. To 
avoid such singularities, restrictions must be placed on the constants defining P. 
Defining pa, pp, a, and /3 so that the roots of P{p) are at Pa + a, Pa — a, pp + iP, 
and Pi3 — ip, one may write P as 

P{p) = -Ci\p - paf - a'){\p - pp? + f3'), (B.4) 

where C = — |- — Cq — 7. Specialize this by assuming that only p^ — a and pa + a 
are real roots, p^ — a < Pq, + a and p^, /3 G M, which means that there are only two 
real roots. Restricting p to lie between these two roots, I reparameterize by setting 
P — Pa + OiCOsO, where as usual 9 e [0, tt]. Then if pp — p^ (that is, P{p) has an 
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axis of symmetry along the line p = pa), potential conical singularities aX pa — a or 



Pa + Oi may be simultaneously eliminated by identifying a with period T 
where P' = 



4tt 



P'(Pc-a) 



Next, I make the following extended series of coordinate transformations and 
definitions: 

1 



Pa 

a 
a 

T 

n 
Q 

Q(r) 



1 + -c?. 



t — acf) 

1 + —0? cos^ 6*, 
+ {poL + « cos ( 



and 



(B.5) 

(B.6) 

(B.7) 

(B.8) 
(B.9) 
(B.IO) 

(B.ll) 

(B.12) 
(B.13) 
(B.14) 



Equating ( |B.3D and ( [B.4| ) obtains the following three equalities relating the two 
forms of P: 



m 



n 



2 I 2 

9o + eo 



2Cp^ 

Cpa{2pl- + P^), and 
Cil + [pl-aVa + P'])-i 



(B.15) 
(B.16) 
(B.17) 
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Then, after a significant amount of algebra, these transformations and equations 
modify the metric (|B.1|) so that it becomes 

.2 . G ,„2 



ds' = A 



g_^a2 I H sin^ e (~, 



adt + [r^ + a^jrf^) + ^dr'^ 



where 

Setting Co 
Q(r) 



A 



[£^ + 2p„cos( 



A/(3C) 



a sin^ 



;i - (A/3)/?2r[p« + acos^])2 



(B.18) 



(B.19) 



^EoP^, go = J l-Go/?^, and Pa = Q becomes. 



A 
3" 



{Ei + Gl){M^(3' - d'){l + ^M^f3^)(3' 



1 - {Ei + Gl)(3^ 



(B.20) 



A 



3 

-2Mr 



A 



2-M ( 1 + -{2M''p^ - a') ] + (1 + -{6M''p^ - a'))r 



A 



3' ' 7 '3 

E^ + GI + (a^ _ M^(3%1 + f M2/3^ 
1 + (Ei + GDP^ 



The term of the above is identified with the NUT parameter. To set this equal 
to zero while keeping the mass parameter M non-zero, one of /3 or 1 + — 



^.2\ 



a 



must be set to zero. Here I choose to take the limit as (3 (choosing 



l + f(M2/3^-52) = results in a metric that is similar to but not quite the KNdS 
metric - most notably it retains a leading conformal factor). Then, replacing a 
with the more traditional symbol a the metric becomes the standard KNdS metric, 
and similarly the vector potential A becomes a vector potential that generates the 
associated electromagnetic field. Thus, the KNdS metric describes two black holes 
in deSitter space that are accelerating away from each other due to the cosmological 
expansion of the universe. 

Note that there are other ways to eliminate the conical singularities in ( p.l 



Although most yield the KNDS metric, some will give rise to other spacetimes. 
They will not be considered in this thesis. 
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B.2 Range of KNdS spacetimes 

This section explores the allowed parameter range of KNdS spacetimes by analyzing 
the root structure of the polynomial Q. 

If Q has three positive real roots then they may be written in increasing order 
a.sd — 6,d + 6,e — 6, and e + e, where e and d are reals and e and 6 are non-negative 
reals. The absence of a cubic term in Q forces d = — e. Two additional conditions 

< e < e, and 

(B.21) 

e < 5 < 2e- e 
ensure that the roots are ordered as proposed. Then without loss of generality 

Q = -^{ir-er-s'){{r + er-5'). (B.22) 

Now, the requirement that Q has three positive real roots enforces restrictions 
on the allowed values of the physical parameters a, M, Eq, and Gq. Q is a quartic 
and therefore can be solved algebraically, so in principle it is possible to directly 
discover under what circumstances it has four real roots. In practice however, the 
exact solution to a quartic is too complicated to work with. Thus, I tackle the 
problem in reverse instead. First I determine the allowed ranges of the Q structure 
parameters e, 6, and e, and then use these to parameterize the allowed range of the 
physically meaningful parameters a, M, Eq, and Go- 

Matching ( |6.9| ) with ( |B.22| ) it is easy to obtain expressions for the physical 



parameters in terms of the structure parameters. Namely 

= j~5^-e''-2e^ (B.23) 

M = ^{6^-e^)e, and (B.24) 

El + Gl = ^{6'-e'){e'-e') + {6' + e' + 2e')-^. (B.25) 
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Requiring that each of these parameters be non-negative imposes further re- 
strictions (beyond the root ordering conditions (|B.21| )) on the allowed ranges of e, 
€, and 5. If > then 

j-6'-e^-2e'>0. (B.26) 

M will automatically be non-negative because of the root-ordering conditions while 
Eq+GI>0 implies that 

-(5^ - e2)(e2 - + (^2 ^ ^2 ^ 2e2) - f > 0. (B.27) 
3 -A. 

In order to disentangle these structure parameters, I rescale them as follows. A 

and e are non-zero so one can define A, E, and X such that 6 = Ae, e = Ee, and 

e = -^/s/AX. Then, the conditions (|B.26|) and (|B.27|) respectively become. 



1-{A'^ + E^ + 2)X^ > 0, and (B.28) 
(A^ - 1)(1 - E^)X^ + (A^ + E^ + 2)X^ - 1 > (B.29) 

The first of these provides an upper bound on the allowed range X for given values 
of A and E. > if and only if 

X <Xu= ^ (B.30) 

~ V2 + A2 + ^2 ^ ^ 



In the meantime, ( p.29|) is quadratic in and so may be easily solved. It 



turns out that over the allowed ranges of A and E, it has only one positive real 
root. Further, it is upward opening, and therefore the positive real root provides a 
lower bound for the allowed values of X. i?Q + Gq > if and only if 



^^^^ = \l 2(A2-1)(1-E2) • ^^-31) 

On plotting Xu and Xl one finds that for < £" < 1 and 1 < A < 2, < X,7 
and so there exists a non-zero range for X for all the possible values of E and A. 
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With this range of allowed values for X in hand, the possible KNdS solutions have 
been fully parameterized. This parameterization is given by the restrictions 

1 < A < 2, < ^ < 2 - A, and Xl < X < Xu. (B.32) 



These ranges are shown in figure [G^. In that figure the allowed parameter range 



of KNdS spacetimes is the region bounded by the five sheets defined by a = 0, 
M = 0, + Gl = 0, E = 0, and E = 2 - A. The last two conditions respectively 
correspond to an extreme (cold) black hole spacetime where the inner and outer 
black hole horizons coincide and a Nariai-type spacetime where the outer black hole 
horizon coincides with the cosmological horizon (though it will soon be seen that 
this apparent degeneracy of the metric is an artifact of the coordinate system and 
that the distance between the two horizons remains finite and non-zero throughout 
the limiting process). The intersection of the Nariai and cold sheets is referred 
to as the ultracold solution. This nomenclature is taken from the corresponding 



non-rotating instantons discussed in 



Having established the range of KNdS solutions allowed by the structure of the 
polynomial Q, it remains to be demonstrated that the extreme limits of the range 
are realizable as a set of well defined metrics. In particular the current coordinate 
representation of the metric breaks down in the Nariai {e ^ 0, 6 ^ 0) and ultracold 
(e — s> 0, 5 — s> 2e — £:) cases. The following three subsections show how these various 
limits may be achieved while the fourth provides some details of the lukewarm 
KNdS solution discussed in section |6.2.3| . 



B.2.1 The cold limit 



This limit can be taken within the current coordinate system. Therefore, the metric 
keeps the form ( |6.8| ) and the electromagnetic field and potential remain as ( |6.10 ) 
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and ( |6.11| ) respectively. The physical parameters are given by: 

= J - 2{3e^ - 2ee + e^) (B.33) 
4A 

M = —e^{e-e), and (B.34) 
o 

Eo+Gl = ^(3e-e)(e-e)2(e + e) + 2(3e2-2ee + £2) - |, (B.35) 
where the range of the parameters is limited by the relations 

< E < 1, and (B.36) 

-3 + 2E-E^ + 2V3^lET2W^ X ^ , ' (B.37) 



(3 - E){1 + E){1 - Ey - - ^2(^2 -2E + 3) 

As before, e = ■\J~^X, and e = Ee. 

In this spacetime, the double horizon of the black hole recedes to an infinite 
proper distance from all other parts of the spacetime (as measured in the hy- 
persurfaces) . Thus, the global structure of the spacetime changes - in particular, 
the region inside the black hole is cut off from the rest of the spacetime. Choosing 
the global structure so that the spacetime contains two (in this case extreme) black 



holes, this structure is shown in figure |6.3| . Note that in this case, the hypersurfaces 
of constant t consist of two extreme black holes, and so are not closed as they were 
in the lukewarm case (the horizons recede to infinite proper distance from all other 
points in the spacetime). 

Finally, note that for the cases where a = 0, this solution reduces to the cold 



solutions discussed in [71]. 



B.2.2 The Nariai limit 



The coordinate system breaks down in the e = limit. Specifically, for e = 0, r = e 
(becomes a constant), and Q = 0, so the coordinate system becomes degenerate. 
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and the metric ill-defined. These problems may easily be avoided however, if one 
makes the transformations 

r = e + ep, (B.38) 

(j) = ^+ " J , and (B.39) 

+ a"^ 

t = (£!±i?!M,. (B.40) 

€ 

Then, the e — limit may be taken without hindrance, and the metric becomes 

ds^ = -QGdr^ + idp^ + ^de^ + ( 2aepdT + ^^d^] ' , (B.41) 

Q T-C G \ J 

while the electromagnetic field becomes, 

—X F sin / + \ 

F = —^dp A rfr H i^:^d9 A 2aepdT H — d<^ . (B.42) 

y \ X J 

An electromagnetic potential generating this is 

^ (e^ - a^) , aEoesin^ 6 + Go{e'^ + a^) cos6' / , + \ 
A = -Eo \ ' pdT 2 2^ 2aept;r + —d^ . 

(B.43) 

In the above, Q = f (2e - - p2)(2e + 5), ^ = + cos^ ^, T = - cos^ ^, 
X = EqT + 2aGQe cos 6, and Y = GqT — 2aEoe cos 6. Note that the above potential 
is not the simply (|6.11| ) under the coordinate transformation as the A generated 



in that way diverges when e 0. The divergence is removed (and the above 
result obtained) if one makes the gauge transformation A ^ A — ^^dr before the 
coordinate transformation and limit. 

The physical parameters are given in terms of e and 5 as 

= ^-2e^-6^ (B.44) 

M = ^(5^6, and (B.45) 

E'o+Gl = ^(52-eV + (2e^ + ^^)-|, (B.46) 
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and the allowed ranges of e = \ \X and 5 = Ae are given by 



1 < A < 2 and (B.47) 



'-(A^ + 2) + AVA^ < X < (B.48) 
2(A2 - 1) - - v/2TA^ 

Note that the Nariai solution is no longer a black hole solution. Extending the 
metric through the horizons by the standard Kruskal techniques, its Penrose dia- 
gram appears as in figure |6.4| (for the (t,p) sector). There is no longer a singularity 
at finite distance beyond either of the horizons. In fact, the diagram is the same 
as that for two-dimensional deSitter space. If there were no rotation (a = 0), then 
this spacetime would just be the direct product of two-dimensional deSitter space, 
and a two-sphere of fixed radius. With rotation, of course the situation is not so 
simple. 

If a = 0, and one makes the coordinate transformation p = cosx, then this 



solution reduces to the non-rotating charged Nariai solution considered in |]7T|. 



B.2.3 The ultracold limits 

Finally consider the ultracold limits where both e ^ and 5 — > 2e — e. It turns 
out that there are two such limits which I label the ultracold I and II limits. In 
this subsection I only demonstrate how they may be reached from the Nariai limit. 
Similar coordinate transformations (which sometimes must be iterated two or three 
times) allow one to reach the same two limits both from the cold limit, and, taking 
6 ^ 2e ~ e and e — > simultaneously, straight from the non-extreme standard 
KNdS form of the metric. I deal with the two cases separately. 

Ultracold I: Making the transformations, 

p = ri-rik{2e-5)R, (B.49) 
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2 

m = $_2ri^^, and (B.50) 
e'^ + 

where 77 = ±1, and k — 8|^e, and taking the hmit as 5 — > 2e one obtains the metric, 

^ -QRdT'^ + ^(ii?2 + £(i^2 ^ ^ ^ / 2aeRdT + . (B.52) 

R n Q \ X J 

The electromagnetic field and potential become, 

~X ysin^ / p2 + \ 

F = -—dR AdT + —^d9 A 2aeRdT + , d$ , (B.53) 

^ V X J 



and. 



+ g{e^ + a^) \ X J 

(B.54) 

e (0,00), T G (—00,00), G [0,7r], and $ inherits a 27r periodicity from its 
predecessors. Q, 7i, x^, X, and F all retain their old definitions. Note that the EM 
potential and field have retained their Nariai form. 

The {R,T) sector of the spacetime is conformally the same as the Rindler space- 
time (which of course is actually a sector of two-dimensional Minkowski space). 
The Rindler horizon is at = and as this is the only horizon, the space does not 
contain black holes. Before giving the parameterization of this solution, consider 
the transformations leading to the ultracold II case. 

Ultracold II: Making the transformations. 



p = b + kV2e - 5R, (B.55) 
, aebx^r 

T 



= $ - 2^7;^, and (B.56) 
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where b ^ ±1, and k = 2 J ^{1 — b'^)e and taking the hmit as 6 ^ 2e, one obtains 



^d9^ + 11 sin2 9 ( 2aeRdT + 

H g \ X 



ds' = -gdT' + gdR' + ^dO' + —sin'e [2aeRdT + ^d^ ] . (B.58 



The electromagnetic field and potential again take the forms ( p.53|) and (p.54| ). 



R,T E (— oo, oo), 9 G [0, tt], and $ inherits a period of 27i from its predecessors, g, 
Ti, X, and Y again retain their meanings from the Nariai case. 

Clearly the {R,T) sector of this spacetime is conformally the same as two di- 
mensional Minkowski fiat space. There is no horizon structure and therefore no 
black hole. 

The physical parameters in both of these cases are given by 

= ;^-6e^ (B.59) 

M = 4-e^ (B.60) 
3 

E^^ + Gl = Ae^ + 6e2-|, (B.61) 
and the allowed range of e = \ j^X is given by, 



1 + -|< X <-L. ,B.62) 

Once more note that when a = these ultracold cases reduce to the two non- 



rotating ultra-cold solutions considered in [71|. However, neither of these spacetimes 



contains black holes. Perhaps one can make an argument for them decaying like 
the Nariai metric into black hole spacetimes, but in any case for completeness I 
shall include them in my considerations throughout the thesis. 
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B.2.4 The lukewarm solution 



As discussed in section |6.2.3| , the lukewarm solution is defined as a KNdS solution 
where the black hole and cosmological horizons are in thermal equilibrium. Their 
temperatures are given by equations ( |6.13| ) and a little algebra shows that they are 



equal (and not degenerate) when 26^ — 2a'^ — e"^ — 6"^ = 0. This relation can be used 
to eliminate 6 from the parameterizations of the physical parameters. Then 

= 4e^ - J, (B.63) 
M = 2e(l - ^(Se^ + ^2)), and (B.64) 
El + Gl = -^{7e' + e'Ke'-e')-2{e'-e') + ^. (B.65) 

The expression for the charge may also be written as E^ + Gn = ^ — a^x^- 
The range of the parameters is limited by the relations: 

< ^ < 1 (B.66) 
<X< ^ (B.67) 



V5 - 2E - E^ ~ \ E^ + 7 



1 2V2- E^ -1- E^ , 

where as earlier e = Ee and e = y^X. The second condition above is the 1 < A < 
2 — E inequality for this case, while the third is the > 0, E^ + Gq > condition. 



Plotting the two conditions over the allowed range of E one finds that ( p.67|) is 



redundant, and so the lukewarm range is given by the first and third conditions. 

These spacetimes are non-extreme KNdS spacetimes, and so have the global 
structure displayed in figure |6TI| . This spacetime was first discussed in Just as 



for the other special KNdS spacetimes that I considered, in the absence of rotation 
the lukewarm case reduces to its non- rotating counterpart discussed ||7T||. 
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